NASA-CR -203951 


HP i5g’-32?57 

■ ’ i ; , , V . « v * 

Computer Aided Grid Interface 


An interactive CFD Pre-Processor 
Final Report 

(performance period: 16 May 1991-15 May 1996) 

Technical Monitor 
Mr. Ted Benjamin 


Sponsoring Agency 
CFD Branch 
NASA Marshall Space 
Flight Center 
MSFC, AL 35807 

Principal Investigator: 

Dr. Bharat K. Soni 
Professor, Aerospace Engineering 


i 




Table of Contents 


Summary 

Introduction 

CAGI Development 

CAGD Techniques in Grid Generation 
A Ph.D Dissertation 

CAGI: User’s Guide 

Appendix 

A.l 

GENIE ++ : Brief Description 




Summary 

NASA maintains an applications oriented computational fluid dynamics (CFD) efforts complemen- 
tary to and in support of the aerodynamic-propulsion design and test activities. This is especially 
true at NASA/MSFC where the goal is to advance and optimize present and future liquid-fueled 
rocket engines. Numerical grid generation plays a significant role in the fluid flow simulations uti- 
lizing CFD. An overall goal of the current project was to develop a geometry-grid generation tool 
that will help engineers, scientists and CFD practitioners to analyze design problems involving com- 
plex geometries in a timely fashion. This goal is accomplished by developing the CAGI : Computer 
Aided Grid Interface system. The CAGI system is developed by integrating CAD/CAM (Computer 
Aided Design / Computer Aided Manufacturing) geometric system output and / or Initial Graphics 
Exchange Specification (IGES) files (including all the NASA-IGES entities), geometry manipula- 
tions and generations associated with grid constructions, and robust grid generation methodologies. 
This report describes the development process of the CAGI system. 


Introduction 

A multitude of techniques and computer codes have been developed for generating computational 
grids in arbitrary regions. However, in most of these codes and methodologies, the evaluation of 
the geometry input and realization of mapping between physical and computational space allowing 
appropriate zonal/block strategies are long, very laborious, and extremely time consuming. Geome- 
try-grid generation is considered as the most time and cost critical part in the typical CFD applica- 
tion. The development of CAGI was initiated in 1992 to fulfill simulation of complex geometries 
encountered in propulsion problems at MSFC in a timely fashion for engineering design. CAGI is 
developed in a modular fashion. A self explanatory pictorial view of different modules and their 
linkage is provided in Figure 1. The computer languages FORTRAN 77 and C in a Unix environ- 
ment are utilized for portability. An X-based interface is provided allowing SGI’s GL (Graphics 
Library) or X-Graphics access by command line argument. The dynamic allocation of memory 
and linked lists on a well-defined data structure is allowed in the CAGI development. The module, 
IGES transformer is developed to facilitate all critical geometrical entities. All NASA-IGES geo- 
metrical entities have been considered in this development. The Non-Uniform Rational B-Spline 
(NURBS) curve/surface/volume representation is selected for the geometric description. The 
NURBS offers a control point based parametric representation which is widely utilized for interac- 
tive design applications using CAD/CAM systems. An inverse formulation of NURBS is developed 
for evaluation of control points associated with the sculptured discretized specification of geometri- 
cal entities. The application of these developments in surface grid redistribution, adaptation, remap- 
ping and optimization are explored. 

The final version of the CAGI system offering an IGES transformer, conversion of the geometric 
entities into standard NURBS data structure and various geometry manipulation and generation was 
released to NASA/MSFC personnel in April 1996. The IGES transformer module was incorporated 
in the grid system GENIE++ and the National Grid System developed at the Mississippi State Uni- 
versity (MSU). The CAGI system has been applied to various grid generation problems based on 
the IGES supplied geometric information. 

The detailed development of the CAGI is published as a Ph.D Dissertation entitled, ”CAGD Tech- 
niques in Grid Generation” in December 1996. This dissertation is included under the section on 
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CAGI development. A user’s guide is also included in the same section. This user’s guide can be 
electronically accessed from the following address: http://www.erc.msstate.edu/thrusts/grid/cagi/ 
A brief description of the computer code GENIE++ is provided in Appendix A-l. 
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Figure 1. CAGI Modules 
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The objective of this study is to develop the algorithms for static and 
dynamic grid generation applicable to complex industrial configurations by 
utilizing Computer Aided Geometry Design (CAGD) techniques. The Non- 
Uniform Rational BSpline (NURBS) is used as a basis for the geometric defini- 
tion and in the development of grid generation schemes. The algorithms 
which bridge the gap between CAD/CAM systems and grid generation systems 
have been developed by utilizing the IGES (Initial Graphics Exchange Specifi- 
cation) file format, which is the output of the CAD/CAM systems. The IGES 
and NASA/IGES prescribed geometric entities have been supported by devel- 
oping algorithms which transfer all geometric entities into a common NURBS 
data structure. These robust transformation algorithms developed in this 
study provide enhancements and generalizations to the existing CAGD tech- 
niques pertinent to the grid generation process. The algorithms to interactive- 
ly construct NURBS curves, surfaces and volumes widely applicable to CFD 
configurations are presented. Projection and inversion techniques applicable 
to an interpolated sculptured discretized data set are developed and validated. 
The reparameterization algorithms to overcome the influence of NURBS geo- 



metric characteristics on the quality of the resulting grids have been devel- 
oped. These algorithms have been cast into a computer software CAGI (Com- 
puter Aided Grid Interface) to facilitate the geometric modeling and surface 
preparations associated with the grid generation process. The CAGI software 
provides a menu driven interactive environment based on the FORMS and GL 
graphics libraries. The interactive graphics capabilities of CAGI have been 
applied to initiate the surface grid generation associated with trimmed sur- 
face entity (IGES entity 144). The NURBS based griding algorithms have 
been applied to construct dynamic grids associated with complex solution 
adaptive and temporally deforming configurations. Computational examples 
demonstrating the success of these algorithms in allowing the treatment of 
CAD / IGES files, formulating very concise NURBS control polygon (control 
net or control volume) associated with industrial geometries, and addressing 
grid redistribution, solution adaptation and dynamically deforming grids are 
exercised. The algorithms developed in this study along with efficient comput- 
er memory requirement and fast evaluation . of NURBS entities make the 
CAGI program very attractive for addressing grid generation needs associated 
with complex industrial configurations efficiently and economically. 
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CHAPTER I 
INTRODUCTION 


With the advent of supercomputers and high powered workstations, 
computational field simulation has become a more common exercise in modem 
analysis, design and manufacturing. This simulation procedure can be di- 
vided into three steps: grid generation (pre-processing - generation of a dis- 
crete representation of a surface or volume for the solution domain), solution 
algorithms (processing - numerical solution of equations of fluid mechanics), 
and the scientific visualization ( the post-precessing - interpretation of simu- 
lated physical field characteristics). Numerical grid generation is usually the 
most labor intensive part of any Computational Field Simulation (CFS) ap- 
plication. In fact, at present it can take significantly more labor time to 
construct a typical CFS grid than it does to execute the flow simulation code 
on the grid or to analyze results. A multitude of techniques for grid generation 
with increasing capability has been developed over the past decade [Ref 57, 
59, 63 —66]. Among the more notable achievements are elliptic equation based 
smoothing procedures [Ref 71—75] and algebraic interpolation schemes [Ref 
57, 60-63] for the generation of surface and volume grids. However, 80%~90% 
of the grid generation labor time is spent on geometry processing. For a rou- 
tine application of CFS, in an industrial environment, the overall response 
time for CFS must be reduced considerably. As noted by Ives [Ref 36]: 

“ The industrial requirement is for reliable one hour grid generation 
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turnaround for one-time geometries when run by designers. The system must 
include CAD-to-grid links which resolve tolerance issues and produce grids 
with a quality good enough for the flow solver. The designer has to feel that 
the grid generation processes is under control and is predictable.” 

In many of today’s industrial applications, most of the geometrical con- 
figurations of interest to practical problems are designed using a CAD/CAM 
system. As pointed out by the NASA Steering Committee on Surface Modeling 
and Grid Generation [Ref 6,35,891, the linkage between the CAD/CAM sys- 
tems and the grid generation systems will significantly reduce the overall turn 
around time for CFS applications. Unfortunately, there are many different ge- 
ometry output formats which force the designer to spend a great deal of time 
manipulating geometrical entities in order to achieve a useful sculptured geo- 
metrical description for grid generation. In addition, there is a danger of los- 
ing the fidelity of the geometry in this process of data transfer between differ- 
ent I/O formats [Ref 6,881 . The other issue related to field simulation is the 
grid quality. It is well known that the quality of the grid affects the accuracy 
of the solution and the computation time [Ref 1, 73]. It may be necessary to 
reconstruct the grids for a more satisfactory result after obtaining the first 
solution. This reconstruction procedure involves a change of either resolution 
(the size of the grids) or the spacing (the distance between grid points) func- 
tions. However, this process, especially true for a complex grid, is tedious and 
very time consuming. This provides the basic motivation for this study. 

To bridge the gap between the CAD/CAM systems and grid generation 
systems, it is necessary to establish the communication paths so that the ge- 
ometries and grids defined within these two systems can be linked with each 
other. The Initial Graphics Exchange Specification (IGES) [Ref 35] is a widely 
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accepted standard for the geometry exchange. All CAD/CAM systems support 
the IGES format as an Input/Output of resulting geometries. Therefore, in 
this study, IGES format has been considered as a standard format for geome- 
try specification. The description of IGES format is provided in the following 
section. Regarding the grid quality and efficiency in geometry / grid manipu- 
lation, it is essential to develop grid generation schemes in conjunction with 
well defined analytical / semi-analytical sculptured (parametric) geometric 
representation. This parametric geometry definition must 'facilitate the gen- 
eration / manipulation / refinement of the high quality geometry (grids). The 
Non-Uniform Rational BSpline (NURBS) [Ref 18,43~47] has been selected as 
the standard for parametric geometric representation in this study. The at- 
tributes and characteristics of the NURBS are provided in the following sec- 
tion. Hence, in this study, the IGES file I/O format and the NURBS represen- 
tation are utilized as the standard for CAD/CAM system output and geometric 
definition respectively. 


IGES Format 

The IGES format was defined by researchers and engineers from indus- 
try, government organization and academics in the late 1970’s and early 
1980’s. This format is intended to describe every aspect of a geometric data set 
(such as the coordinates of the data set, what color to be plotted and even the 
line width to be plotted ...), and is designed to be the national standard of In- 
put/Output format for CAD/CAM systems. There are five sections in each 
IGES file, namely, the Start section, Global section, Directory entry section, 
Parameter data section and the Terminate section. All entities described in an 
IGES file are labeled with an integer number (for example, a circular arc is 
labeled with 100) and classified into four classes. These classes are “curve and 
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surface geometry entities”, “constructive solid geometry entities”, “annotation 
entities” and “structure entities”. Among these four classes, the widely uti- 
lized class for CFS engineers and designers is the “curve and surface geometry 
entities”. These entities are labeled from 100 ~ 199. An example of an IGES 
file describing the definition of a sphere is listed in Appendix A. In this IGES 
example, the Start section includes the human-readable prologue to the file. 
All records in this section shall have the letter “S” in column 73 and a se- 
quence number in columns 74 through 80. The Global section of this file con- 
tains the information describing the preprocessor and information needed by 
the postprocessor to handle the file. All records in this Global section shall 
contain the letter “G” in column 73 and a sequence number in columns 74 
through 80. Many important parameters are defined in this section, two of 
which are the characters of parameter delimiter and record delimiter. The de- 
fault characters are the “comma” and “semicolon” respectively. The Directory 
entry section has one directory entry for each entity in the file. This directory 
entry for each entity is fixed in size and contains twenty fields of eight charac- 
ters each, spread across two consecutive eighty character lines. All records in 
this section contains the letter “D” in column 73 followed by a 7 digit sequence 
number. For this example, the IGES file describes a sphere defined by a 
NURBS surface; hence, the entity number 128 is shown in the first and elev- 
enth fields of this directory entry. The Parameter data section contains the 
parameter data associated with each entity defined in the Directory entry sec- 
tion. For this example, the NURBS information, such as the control vertices, 
the knot vector and the weights, Eire defined in this section. The Terminate 
section is the last section. It contains a character “T” in column 73. There is 
only one fine of the Terminate section in an IGES file. This section describes 
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the number of lines used in the Start section, Global section, Directory section 
and Parameter data section in the IGES file. 

As aforementioned, the IGES file is an attempt to represent every as- 
pect of geometry data, and thus its documentation is long and difficult to com- 
prehend. Utilization of the IGES data file requires expert knowledge of the 
data format and the interpretation associated with the geometry description 
of respective entities. Hence, even though this format has been defined since 
1980*sdt has not been widely used by the grid generation community. To en- 
hance the utilization of an IGES file and to expedite the entire grid generation 
process, scientists at various NASA centers (NASA Ames, NASA Langley and 
NASA Lewis) have formed an IGES committee. This committee, referred to as 
the NASA IGES committee, selected several commonly used entities in CFS 
simulation as a group. This group of IGES entities is known as the NASA- 
IGES. NASA-IGES provides a common standard for geometric description 
between the CAD/CAM and the CFS community. All these efforts simply show 
the importance of the IGES format in creating quick and efficient results for 
the geometry transformation. 

In view of these efforts, IGES has been selected as the CAD/CAM sys- 
tem for geometry Input/Output in this study. The algorithms to interpret, 
analyze and manipulate the IGES geometric entities by utilizing the NURBS 
as a standard geometric representation have been developed in this work. 

NURBS 

NURBS stands for Non-Uniform Rational BSpline and is a parametric 
form consisting of the control polygons, knot vectors, orders and weights to 
represent the geometric entity accurately. If all the weights of a NURBS enti- 
ty are equal to one, then it is called a Non-Rational BSpline (or simply 
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BSpline). The BSpline was developed in the 1970’s and 80’s by Carl de Boor, 
M. G. Cox and J. Schoenberg. In 1973, R. Riesenfeld presented' the paper 
entitled “Applications of BSpline approximation to geometric problems of com- 
puter-aided design” which was the first application of the BSpline applied to 
CAD. The properties and applications of NURBS were initialized and exten- 
sively developed during the 1980s. The first commercial product which used 
the NURBS to represent the geometry came from the SDRC (Structural Dy- 
namics Research Corp.) in 1983. The Boeing company proposed the NURBS 
as an IGES standard in 1981, and now the NURBS curve and NURBS surface 
have been adapted as the IGES geometric entities 126 and 128. 

Recently, the IGES entities 126 and 128 have become increasingly pop- 
ular in geometric definition. This popularity is due to various characteristics 
of the NURBS, especially the shape preserving property. Many commonly 
used curve and surface definitions can be analytically represented by NURBS. 
These commonly used geometric entities include the circular arc, conic arc, 
line, cubic spline curve, cylinder, ruled surface, surface of revolution, paramet- 
ric surface ... and so on. In addition to this shape preserving property, NURBS 
has many other powerful features, such as the convex hull, local control, varia- 
tion diminishing and affine invariance. The convex hull property ensures that 
the NURBS curve (surface or volume) will lie entirely inside its associated 
control polygon (control net or control volume) [Ref 23-25,34]. The local con- 
trol property facilitates the local modification of the geometry without altering 
the entire geometric shape. This attractive property makes it possible to re- 
fine / reconstruct the desired portion of a geometry locally. The variation di- 
minishing property, taking a NURBS curve as an example, is defined such 
that if a straight line intersects the NURBS control polygon of a planar 



NURBS curve k times, then it can intersect the NURBS curve at most k times. 
This property ensures that the NURBS is always convex whenever its control 
polygon is convex, and the NURBS curve can be convex even when its control 
polygon is not convex. Furthermore, if a NURBS curve has a point of inflec- 
tion, then its associated control polygon must have at least one point of inflec- 
tion [Ref 23,341. This property is useful in predicting the shape of a NURBS 
representation. 

Also, the NURBS geometry tool kits which include the knot insertion, 
degree elevation and splitting algorithms, make NURBS useful in their own 
right. The knot insertion algorithm is used to increase the associated knot 
vector without changing the original NURBS shape. After increasing the 
number of knot values in a knot vector, the designer has more freedom to 
modify the associated NURBS curve (surface or volume). Similar to the knot 
insertion algorithm, the degree elevation provides the flexibility to modify the 
NURBS entity by increasing the associated degree without changing the origi- 
nal geometric shape. These two functions are important, . fundamental and 
unique to the NURBS representation. The splitting algorithm allows the de- 
signer to “cut” a NURBS entity into two different NURBS entities. 

All these geometric properties make the NURBS stable for representing 
complicated data sets and also very attractive to the designers and engineers 
in CFS and CAD/CAM. In fact, because of these attractive properties, many 
efforts have been proposed utilizing the NURBS as the standard and geomet- 
ric representation. For example, after presenting the the NASA— IGES for- 
mat, the NASA IGES committee further proposed the NINO (NASA IGES 
NURBS ONLY) standard [Ref 6,35]. The scientists and engineers from the 
Naval Surface Warface Center and Boeing company combined their efforts 
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and designed the software package “DT_NURBS” [Ref 21]. This software 
package utilizes NURBS definitions as a standard data structure for geomet- 
ric entities and provides evaluations and manipulations of NURBS entities. A 
general purpose grid generation package NGP ” [Ref 29,71] developed at Mis- 
sissippi State University / NSF Engineering Research Center also utilizes the 
NURBS as the standard database. Widely used CAD/CAM packages like “PA- 
TRAN”, “I-DEAS” and “ICEM CFD” all claim their products support the 
NURBS geometric database. Besides the IGES format, many national and in- 
ternational standards, such as the STEP and PHIGS, recognize NURBS as a 
standard for geometric design. Even in applications of computer graphics, 
NURBS is used to model the complicated objects. In addition, NURBS has 
even been implemented and supported as hardware functions. For example, 
the GL (Graphic Library) [Ref 39] in Silicon Graphics platform supports hard- 
ware functions such as “nurbscurve”, U nurbssurface ” , a nurbsc” and “nurbss ”, ... 
and so on. These functions account for the hardware implementation of 
NURBS curves and surfaces definitions. Similar NURBS functions such as 
the "gluGetNurbsProperty” y “ gluNurbsCallback ”, “ gluNurbsCurve ” and “glu- 
NurbsSurface” are supported by OpenGL [Ref 39] - smother popular graphics 
library. Even in virtual reality applications, NURBS has been used for model- 
ing scenes. In fact, the NURBS representation is becoming the “defacto” stan- 
dard for the geometry description in most modern CAD/CAM, grid generation 
systems and computer graphics. 

The application of NURBS to grid generation was first presented by 
Yoon [Ref 81] in 1991. However, in his study, most of the concentration was 
placed on the BSpline (by considering the case of all weights to be one). Fol- 
lowing Yoon’s work, Yang utilized NURBS for the application of grid adapta- 



tion [Ref 83] in 1993. Shih applied the NURBS for data interpolation in the 
turbomachinery [Ref 53] in 1994, Craft used it for re-splining the surface 
grids for a CFS analysis of a complete aircraft [Ref 11,67], and Boyalakuntla 
used it for simulation of the temporally deforming geometries [Ref 13]. These 
researchers have contributed in utilizing the NURBS interpolation and evalu- 
ation routines developed under this study in the area of numerical grid gen- 
eration. However, the NURBS has been used only as an interpolation tool by 
these researchers. As a matter of fact, those researchers only utilize the 
BSpline definition instead of NURBS. In this study, the detailed analysis and 
algorithms to utilize NURBS in grid generation are presented. 

As is mentioned, the communications between the CAD/CAM and grid 
generations are crucial. In order to communicate the geometry from grid gen- 
eration to CAD/CAM, the grid generation system must have the ability to read 
in / output to IGES file. The transform of those geometric entities, which are 
defined in the IGES standard, to discrete grid formats (such as PLOT3D for- 
mat) has been done in [Ref 86, 87]. As noted, these discrete grids formats have 
several drawbacks. They will lose the parametric definition, and the grid is 
difficult to re-evaluate or change the spacing. Also, this discrete grid format 
can not communicate to a CAD/CAM system anymore. However, many of the 
geometric representations can be transformed to NURBS analytically. And 
the NURBS can be used in grid generation to model various geometric defini- 
tions, it is possible to communicate with the CAD/CAM system through the 
IGES entities 126 and 128 which are the NURBS curve and NURBS surface. 
Hence, the NURBS representation provides the linkage between the grid gen- 
eration and CAD/CAM system. 
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This study concentrates on the development of generalized algorithms 
so that the NURBS can be fully utilized in grid generation system, and thus, 
eliminate the communication hurdle. The objectives and the organization of 
this study follows. 

Objective and Organization 

In this study, the objectives are to develop the algorithms to bridge the 
gap between CAD/CAM systems and numerical grid generation systems. The 
CAGD techniques have been utilized in CAD/CAM systems for decades. Many 
of them are well documented in the related literature. Several existing algo- 
rithms have been adopted in this research. These basic and fundamental al- 
gorithms include the (NURBS) knot insertion, degree elevation, splitting / 
joining algorithms [Ref 8,42~47, 76,77] and the evaluation of the BSpline ba- 
sis functions [Ref 18]. Other existing algorithms, such as the data reduction 
routines [Ref 38] and the FORM Library [Ref 39] which are available in the 
public domain, have also been utilized in this study. 

This work is carried out by first including the development of algo- 
rithms for transforming Non— NURBS entities encountered in standard IGES 
format to NURBS definitions and then following the development of algo- 
rithms for easy and efficient NURBS manipulation. The transformation algo- 
rithms presented in this study provide the enhancements and generalizations 
of those described in CAGD literature [Ref 7,42,43] . For example, the trans- 
formation algorithm associated with a circular arc to a NURBS curve is en- 
hanced to facilitate an arbitrary arc with different sector angle (the difference 
of ending angle and starting angle) without extra computation of using the 
knot insertion algorithm [Ref 45]. The transformation of the conic arc to a 
NURBS curve representation is generalized by solving an implicit equation 
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which provides more information necessary for grid generation process. This 
information also includes the types of the arc, the orientations, the sector 
angle and the semi-major and semi-minor axis. Also, the transformation al- 
gorithm developed for the surface of revolution is not limited to the full revolu- 
tion which is the only case described in CAGD literature. 

In addition to the transformation algorithms, the NURBS generation 
algorithms are also developed. In the CAD/CAM and CAGD area, the NURBS 
is only used for generating the curves and surfaces. However, the NURBS 
curve and surface algorithms have been extended to the volume generation in 
view of the application to grid generation. These new algorithms can be ap- 
plied for various 3D NURBS volume generation algorithms tailored for many 
CFS related configurations. Also, the modeling techniques for the superellip- 
tic arc and cascading surface by NURBS curve and surface presented in this 
study have not been discussed before in the literature. These algorithms can 
facilitate the use of CAD/CAM data for CFS analysis. The most commonly 
used NURBS evaluation algorithm discussed in related literature is the “de 
Boor” method [Ref 23, 34]. In this study, a new evaluation algorithm for 
NURBS representations (curve, surface and volume) is developed. This new 
algorithm provides a competitive alternative to the de Boor algorithm in terms 
of the memory used and the computational time. Since NURBS is a paramet- 
ric representation, the determination of proper parametric values for the de- 
sired (smooth) grids on NURBS entities is a difficult issue. The new repara- 
meterization algorithms developed in this study provide an efficient approach 
to solve this difficulty 

The grid adaptation is frequently used in many CFS applications. It 
provides the ability for obtaining more accurate solutions. Since the repara- 
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meterization algorithms can be used for precise grid distribution control on 
the NURBS entities, this reparameterization algorithm is then used to en- 
hance the adaptation algorithm developed by Yang [Ref 83 ~85] in this study. 

The organization of this dissertation is as follows. First, the back- 
ground and the motivation of this study are introduced in Chapter One. The 
NURBS descriptions along with the algorithms for transforming curves and 
surfaces to NURBS representations are described in Chapter Two. The algo- 
rithms for generating 3D surface and volume grids by various approaches uti- 
lizing NURBS control net and control volume are presented in Chapter Three. 
The general BSpline interpolation algorithm along with the projection and in- 
version methods follow. The NURBS re-parameterization algorithms are de- 
scribed in Chapter Four. These algorithms are designed to maintain the dis- 
continuity of the NURBS entities and accomplish a better grid distribution on 
the physical NURBS entities. A new algorithm developed for a fast and effi- 
cient NURBS evaluation is also included. The applications of NURBS in dy- 
namic grid generation are demonstrated in Chapter Five. These applications 
include the grid adaptation and temporally deforming geometry. In this chap- 
ter, the NURBS “local control” property and the re-parameterization algo- 
rithm are fully demonstrated. The grids are modeled with a concise NURBS 
control net (or control volume), and computational examples are included. 
The overview of the computer software package CAGI is presented in Chapter 
Six. The modules of CAGI and the associated functionalities are introduced. 
There are many packages which claim the capability of reading the IGES file 
for numerical grid generation. However, very few of them can handle two com- 
plicated and difficult entities — the bounded surface (entity 143) and the 
trimmed surface (entity 144) for “structured” grid topology. The approaches 
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and strategies of handling these two entities and the examples are shown in 
Chapter Seven. The summary and conclusions are listed in the last chapter. 



CHAPTER II 

NURBS AND THE TRANSFORMING ALGORITHMS 


NURBS Formulation 

The definitions of NURBS curve, surface and volume are presented as 
follows: 

A NURBS Curve of order k is defined as: 

f W. d. NHt) 

at) = ^ (2.D 

jr w t jv*«) 

1=0 

where the cfy i=0,..,n denotes the deBoor control polygon and the Wi are the 
weights associated with each control point. The Nf(t) is the normalized 
BSpline basis function of order k and is defined over a knot vector T = 3} i= 0, 
... , n+k by the recurrence relations as shown in equation (2.2). 
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Throughout this study, it is assumed that the knot vector has the form T = (0, 
... , 0, Tk , ... , T n , 1 , ... , 1} with the multiplicity k for the knot value 0 and 1 on 
both ends of the knot vector. If the knot vectors do not match this format, the 
knot insertion [Ref 8] technique must be used to achieve the multiplicity of k 
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on the ends of the knot vector, and if the end knot values are not 0 and 1, the 
knot vector must be normalized by the last knot value to match this format. 
Because shifting and scaling (normalizing) the knot value will not alter the 
shape of geometry, the basis function list in equation (2.2) is defined as “Nor- 
malization” basis function. 

The NURBS surface is the extension of the curve from ID to the 2D ten- 
sor product parametric surface and is shown as equation (2.3). 


S(s,t) = 
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Where dij denotes the 3D control net and Wy are the weights associated with 
each control point. The N\Hs), N^ 2 (t) denote the normalized BSpline basis 

r J 

functions of order kl and k2 over the two knot vector 2j=2J i=0,..,m+kl and 
22=2} j=0,..,n+k2 in the I and J directions, respectively. The definition of the 
BSpline basis functions of NURBS surface is exactly the same as for the curve 
shown in equation (2.2). 

The formula for 3D NURBS volume is defined analogous to NURBS 
surface and is a 3D tensor product form written as: 
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The di j i form the 3D control volume, and the Wi j k are weights 
associated with each control point. The N^Hs), Nj 2 (t) and N^ 3 (u) are the 
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normalized BSpline basis functions of order kl, k2 and k3 over the two knot 
vectors 7j=7J i=0,..,m+kl, T 2 =Tj j=0,..,n+k2 and 73=7} l=0,..,p+k3 in the/, J 
and L directions (i.e., the s,t,u directions), respectively. 

Transforming Procedures 

Transforming the Non-NURBS geometric curves and surfaces to 
NURBS definition is the main topic of this section. To model a NURBS entity, 
according to the equations (2.1) - (2.4), one should define the control polygons 
(or control net / volume), weights, knot vectors) and the order(s). 

CAD/CAM systems and grid generation systems often utilize different 
formats which makes transfer of information between systems difficult. The 
procedure presented in this section has been developed to facilitate this com- 
munication. A geometric entity defined in a CAD system can be represented 
in many ways (for example, cubic parametric spline surface, ruled surface, 
surface of revolution or extruded surface ...etc.), while still conforming to the 
IGES standard. It is well known that most geometric definition can be analyt- 
ically transformed to a NURBS representation. The transformations of most 
widely used entities to NURBS are described by Piegl [Ref 43,44,46]. Howev- 
er, there are many practical issues which are not covered in the transforming 
procedures published in those literature. For example, the IGES representa- 
tion of the implicit conic arc, an important entity, is not contained in those ref- 
erences. The transforming algorithm for a general circular arc (a circular arc 
with arbitrary starting and ending points) is also missing from those refer- 
ence. Another problem is the available literature may not provide sufficient 
detail. For example, procedures for transforming a surface of revolution into a 
NURBS are provided only for a 360° revolution, but many grid generation ap- 
plications require a specified range, such as 60°. The procedures for trans- 
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forming a ruled surface to a NURBS are only covered for sequence connection, 
but the IGES format defines two ways of connections [Ref 35]. Also, several 
transforming algorithms are never discussed in any of the literature. For ex- 
ample, the Transfinite Interpolation (TFI) for a NURBS volume and the mod- 
eling of the super-ellipse as the NURBS curve. Hence, the following sections 
provide enhancements and generalizations to existing transformations devel- 
oped to meet needs arising from the grid generation process for complex geom- 
etries defined in a CAD/CAM system. 

Transform the Curve Entities to the NURBS Representations 
The algorithms of transforming various curve definitions to NURBS 
curve representation follows. There are six different cases discussed in this 
section. 

Transform Straight Line (Entity 110) to NURBS Curve 

The first (and the easiest) entity to transform is the straight fine (entity 
110). In the IGES documentation, a straight line is defined as the connection 
between two 3D data points. Therefore, with regard to NURBS, the definition 
of the line uses the original two points as the two control points (d). Then, set 
order (k) equal to 2, set weights (W) all equal to 1, and set the knot vector as 
(0,0, 1,1). Since there are only two control points, the n is equal to 1. One may 
argue the need for using NURBS to represent a straight line, because the 
NURBS form needs extra storage than the traditional straight line definition. 
This is why even the NASA-IGES committee feels the straight line definition 
could remain unchanged. However, there are two reasons to transform this 
entity to NURBS form. First, transforming all entities to NURBS form keeps 
a simplicity of the database. Second, if this straight fine is one of the constitu- 
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ent entities of the composite curve (refer the “Transform composite curve to 
NURBS curve” for detail), then representing that composite curve with 
NURBS form requires that the straight line be represented in NURBS defini- 
tion so that one can perform knot insertion or degree elevation to join all con- 
stituent curve entities into one composite NURBS curve. 

Transform Circular Arc (Entity 100) to NURBS Curve 

A circular arc (entity 100) as defined in the IGES standard is repre- 
sented by a center point, starting point and ending point within a given 
constant Z plane. The two end points and the center point form an arbitrary 
sector angle which does not necessarily start from zero. It has been shown 
that any circular arc with sector angle less than or equal to 90° can be repre- 
sented by NURBS [Ref 42,45]. The basic control polygon for this NURBS rep- 
resentation is shown in Figure 2.1. 



Figure 2.1 The basic control triangle for a circular arc. 

In Figure 2.1, C is the center point, S, is the starting point and E is the 
ending point. The sector angle SCE (0) is less than or equal to 90°. The two 
tangent lines SD and ED intersect at D. The order of this control polygon is 
three, with the control points S, D, E (hence, the n is 2) and the weights are 1, 
cos(0/2) and 1 respectively. The associated knot vector is (0,0, 0,1, 1,1). A circu- 
lar arc with sector angle greater than 90° and less than or equal to 180° can be 
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represented by two arcs with one half of the original sector angle. For each of 
these two sections the previous procedure can be used to evaluate the corre- 
sponding control polygon. A 180° circular arc represented by two control poly- 
gons is illustrated in Figure 2.2. 



Figure 2.2 The NURBS control polygon for a semi-circle. 

These two control polygons can be combined and the common point M 
can be eliminated. The resulting NURBS information is setting the control 
polygon to SIMJE (hence, the n is 4), the knot vector to ( 0, 0, 0, 1/2, 1/2, 1, 1, 1) 
and the weights to (1., cos (0/n), 1., cos (0/n), 1.).' A similar procedure can be 
used for circular arcs between 180° and 270° (with n equal to 6) resulting in a 
final knot vector of (0, 0, 0, 1/3, 1/3, 2/3, 2/3, 1,1,1) and weights (1., cos (0/n), 1., 
cos(0/n), l.cos(0/n), 1.), and a knot vector of (0, 0, 0, .25, .25, .5, .5, .75 , .75 
,1,1,1) and weights (1., cos(0/n), 1., cos(0/n), 1., cos(0/n), l.cos (0/n), 1.) for arcs 
between 270° and 360° for n equal to 8. These four cases are shown in Figure 
2.3. 

This approach handles all possible circular arcs with no extra computa- 
tion (such as knot insertion) involved. Furthermore, the parameterization 
(distribution on the curves, see detail in Chapter four) is good for all cases. 




Figure 2.3 Arbitrary circular arcs with the NURBS control polygons. 

Transform Conic Arc (Entity 104) to NURBS Curve 

The transforming procedure for conic arc was discussed in [Ref 42,44], 
where they described the case of 3 given control points and changing the 
weight (conic shape factor) to produce a different family of conic arcs (elliptic, 
hyperbolic or parabolic arc). However, that case is complete different than the 
one defined in IGES format. The conic arc defined in IGES is represented by 
an implicit form Ax 2 + Bxy + Cy 2 + Dx + Ey + F= 0, with starting point S and 
ending point T supplied (counterclockwise). The transforming procedure for a 
basic conic arc is illustrated in Figure 2.4. In this figure m is the middle point 
of line TS. 



Figure 2.4 Basic NURBS control polygon for a conic arc. 

Since the two end points are known, the two slopes of the tangent lines 
at the end points can be obtained. The equations describing these two tangent 



21 

lines can be formed and the intersection point N can be determined. This is 
accomplished as follows: 

Differentiate the implicit form of the conic equation to obtain 2Ax + By 
+ Bxy + 2Cyy’ + D + Ey’=0. Solving this equation for the derivative yields: 

y’= (2Ax + By + D)/ f-2Cy-Bx-E) (2.5) 

Substitution of the coordinates of the two end points S and T into equa- 
tion (2.5) yields the two desired straight lines. The shoulder point h can then 
be obtained by solving for the intersection of the line Nm and the given implic- 
it equation. The control triangle is then defined by the polygon SNT (hence, 
the n is 2 for this case) with weights of (1, (jnh)/(hN), 1). The order can be set 
to 3 and knot vector is defined in a manner analogous to the circular arc. As 
long as this basic control triangle can be found, the procedure used for the cir- 
cular arc with the sector angle greater than 90° can be applied to conic arc by 
simply combining the different control triangles together to form the final con- 
trol polygon and by setting the proper knot vector. The definition of sector 
angle 0 for the conic arc is only applied to the elliptic arc, for the parabolic or 
hyperbolic arcs, three control points are sufficient to form the control polygon. 
Hence, for parabolic or hyperbolic arc, the knot vector is always (0., 0., 0., 1.0, 
1.0, 1.0) with n equal to 2. Figure 2.5 shows different conic arcs represented 
by the NURBS using this algorithm. From left to right, (I): Elliptic arc with 
equation 2x 2 +4xy+5y 2 -4x-22y+7=0, form by two NURBS control polygons. 
(II): Parabolic arc with equation 4x 2 - 4xy +y 2 -2x- 14y + 7 = 0, (III). Hyper- 
bolic arc with equation 2x 2 + 4/3xy - 2y 2 -l 6=0. 



Figure 2.5 NURBS control polygons represent different conic arc. 

Transform Parametric Curve (Entity 112) to NURBS_Curve 

The cubic parametric curve defined in IGES format is a sequence of 
parametric polynomial segments. More precisely, it is composed of N (2V> I) 
pieces of cubic parametric segments as illustrated in Figure 2.6. 



u=T(l) u=T(3) u=T(5) 

Figure 2.6 The definition of parametric curve in IGES format. 


In Figure 2.6, the T(i), i=l , .. , N+l are the breakpoints. For this case, 
N= 5, hence, there are 5 cube parametric segments constitute the final curve. 
For each parametric curve, it is defined as 

C(u) = a + bt + ct 2 + dt 3 T(i) ^ T(i+1) and t = u- T(i) (2.6) 

It has been proven [Ref 23,24,34] that the cubic Bezier curve is a special 
case of a BSpline curve with knots vector of (0, 0, 0, 0, 1, 1, 1, 1) (no interior 
knot value). Also, the BSpline curve is a special case of NURBS curve with all 
weights equal to 1. The mathematical transformation from parametric cubic 
spline curve in IGES definition to NURBS is accomplished as follows: 
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The matrix form of the each simple cubic parametric curve, according to 
equation ( 2 . 6 ), can be expressed as C(t) -[It t 2 t 3 J 14^4 [a be d] T where is 

the identity matrix and [a b c d] T \s the transposed matrix containing the coef- 
ficients of the cubic curve. The matrix form of the cubic Bezier curve is ex- 
pressed as C(t) = [1 1 t 2 t 3 ] B4y£ [bo bi l>2 bs] T . The £4^ is the cubic Bezier 
matrix and [bo bj b3J T is the transpose matrix containing the Bezier con- 

trol polygon. The strategy is to first transform the cubic parametric curve to 
Bezier form, since a Bezier curve can be treated as the special case of a 
NURBS curve. Each segment of parametric spline curve is transformed to a 
Bezier curve by finding the associated Bezier control polygon. This is done by 
setting the the two matrix equations to be equal 

'1 0 0 0] IV 

Bezier= [1 t t 2 t 3 ] ~ 3 3 0 0 

3 -6 3 0 b 2 

-1 3 -3 lj 6 3 

Tl 0 0 01 

0 10 0 6 ' 

= Cubic curve = [1 t t 2 t 3 ] 0 0 10 c- ( 2 . 7 ) 

0 0 0 lj d. 

Solving the equation ( 2 . 7 ) for the Bezier control polygon. Since the cubic para- 
metric spline defined in IGES is composed of N pieces of cubic curves, the 
range of parametric value t for each piece is not the same as for the Bezier 
curve. Hence, a re-parameterization of the cubic parametric curve is neces- 
sary. For each piece of cubic curve, the coefficients [04 ty q di] T can be ob- 
tained from the IGES file, therefore, the final equation to solve (for each seg- 
ment) is 
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( 2 . 8 ) 


where h = T(i+1) - T(i) and T(i) is the break value defined in the IGES file. 
After all the Bezier control polygons have been obtained, one can join them 
together and set the multiplicity of joint knot value equal to 3 to form the final 
Bspline curve. For example, if two cubic Bezier control polygon are obtained, 
the final knot vector will be set as (0, 0, 0, 0, 0.5, 0.5, 0.5, 1, 1, 1, 1) and the 
final curve would be C° continuous with order equal to 4 and all weights equal 
to 1. The knot removal algorithm described in [Ref 5,76] can be applied to 
eliminate the redundant knot and reduce the number of control points. Figure 
2.7 (not applying the knot removal algorithm) demonstrates this approach. 



Figure 2.7 BSpline control polygon for parametric curve with 2 segments. 

Transform Composite Curve (Entity 102) to NURBS Curve 

A composite curve(entity 102) is defined as a curve entity consisting of 
lists of constituent curves. The constituent curve can be any parameterization 
curve except another composite curve. And this entity is a directed curve, 
which means the direction of the composite curve is induced by the direction of 
the constituent curves in the following manner: The start point for the com- 
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posite curve is the start point of the first curve entity appearing in the defin- 
ing list, and the terminate point for the composite curve is the terminate point 
of the last constituent curve appearing in the defining list. Within the defin- 
ing list itself, the terminate point of each constituent curve entity has the 
same coordinates as the start point of the succeeding curve entity. It is quite 
difficult to represent the composite curve precisely without transforming all 
the constituent curves to the NURBS form. After transforming all curve enti- 
ties (like straight fines, circular arcs, conic arcs, parametric curves and ratio- 
nal BSpline curves), the “NURBS Joining” algorithm for all the constituent 
NURBS curves is performed to form the NURBS representation for the com- 
posite curve. The procedure is illustrated as follows: 

Suppose two constituent curves Ci and C 2 (already transformed to 
NURBS definition) form a composite curve. Then the first step is to perform 
the degree of elevation [Ref 14] of the lower degree curve so that the curves 
can have the same order. The second step is to adjust the knot vector of the 
second curve C 2 so that the first knot value of the second curve can have the 
same value as the last knot value of the first curve. Shifting the knot vector 
will not change the original NURBS curve because the basis function is a 
“normalized” basis function. The third step is to build up the final knot vector 
by joining the two knot vectors into one knot vector and set that knot value at 
the joint point to have the multiplicity equal to (< order — 1). For example, if the 
first knot vector is [0,0,0, 1,1,1] and the second knot vector is [2, 2, 2, 3, 3, 3], ad- 
just the second knot vector by shifting -1 to each value. Thus, the second knot 
vector becomes [1,1, 1,2, 2, 2]. Suppose the order of these two curves are 3, then, 
the final knot vector should be [0,0,0,1,1,2,2,21 (one may notice the interior 
knot 1 has multiplicity of {order —1) = 2). The fourth step is to match the 
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weights by timing the ratio of (the last weight of the first curve/ the first 
weight of the second curve) to all the weights of the second curve so that the 
weight at the joint point for the two curve are the same. The last step is to 
build up the final control polygon and weights by throwing the first control 
point and weight of the second curve away and jointing the others as one con- 
trol polygon and one weights vector. After these procedures have been applied 
to all the constituent curves, a composite NURBS curve should be formed. 
One more procedure that may apply to this final curve is to perform the knot 
removal to remove the redundant knot vector [Ref 5,76]. Figure 2.8 shows this 
algorithm for transforming the composite curve consisting of, from right to 
left, one straight line, one circular arc, one straight line, one ellipse arc and a 
straight line. 



Figure 2.8 NURBS control polygon for a composite curve. 

Transform Superellipse to NURBS Curve 

A superelliptic arc can be described as the equation (2.9) 

( | Y> + ( | V = 1 (2.9) 

where a is the semi-major and b is the semi-minor axes of the superellipse. 
Special cases of the equation (2.9) include a circle (with a = b, and rj = 2), an 
ellipse (with a =4= b, and rj = 2) find a rectangle (with a + b, and rj = oo). 
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The definition of superellipse is not included in IGES format, however, 
it is a commonly used geometric description. An example of this is the model- 
ing of a transition duct used for the test of a single-engine nozzle [Ref 41,48]. 
The transition duct was designed by using a sequence of constant-area, super- 
elliptic cross sections according to equation (2.9). In most of the literature, the 
process of obtaining the exponent of the superellipse rj was described as solv- 
ing the implicit function relating the quantities a, b and rj to Ac S of equation 
( 2 . 10 ). 

Acs = ^^(2/i/X 4o6) (2.10) 

In this study, the transforming of this superellipse to NURBS curve is 
presented in a more straightforward way as follows. 

This transforming approach is a combination of the circular arc and the conic 
arc algorithms. 

E=(0,h 

O=(0,0) 

Figure 2.9 Illustration of the NURBS control points for a superellipse. 

Consider a superellipse with semi-major a and semi-minor b in the 
first quadrant shown in Figure 2.9. This arc starts at the point (a,0) and ends 
at the point (0,b). Two tangent lines intersect at the point (a,b). Similar to the 
algorithm of circular arc, these three points can be used as the NURBS control 
polygon while setting the order to be 3 with knot vector (0., 0., 0., 1., 1., 1.). 
The weights at the starting and ending control polygon can be set to 1.0. The 
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only problem left is determining the weight at the middle point D of the control 
polygon. This is done similar to the algorithm of the conic arc. Construct the 
straight line OD and let this line intersect with the line SE and the superellip- 
tic arc at the points of m and h. The weight at point D is then set as the ratio 
of (hm/hD). 

This approach is self-explanatory. When the exponent of the superel- 
lipse r] increases, the arc is changing from a circular arc to a rectangular arc, 
this means that point h is approaching to the control point D. Also, the dis- 
tance of hD is decreasing and as a result the weight at point D is increased. 
This situation matches the NURBS theory — a NURBS curve is pulled to- 
wards to the control point when the weight of this control point increased. The 
mathematic verification can also be done by comparing the h point with the 
shoulder point evaluation from NURBS representation. Since the variables 
(a, b and rj) of superellipse are all given, the h can be solved from the intersec- 
tion of the line OD and the arc. On the other hand, after the entire NURBS 
representation is set up for this superellipse, the shoulder point h can also be 
evaluated with the parametric value t - 0.5. Comparing to the locations of 
these h’s, one can find out that the relative deviation is as small as 1.0e-9. 
Table 2.1 shows the selected values of the exponent rj of the superellipse and 
the corresponding weights values. 



Table 2.1 The relationship between exponent rj and weights. 


n 

weight 

2.000000 

0.7071067807 

2.076143 

0.7615055209 

2.184741 

0.8391550277 

2.310944 

0.9294727665 

2.446475 

1.0265482055 

2.588168 

1.1281144695 

2.736506 

1.2345144266 

2.894152 

1.3476587943 

3.064489 

1.4699782629 

3.250206 

1.6034070829 

3.453315 

1.7493976138 

3.676614 

1.9099667660 

3.924127 

2.0880154404 

4.201364 

2.2875017047 

4.515468 

2.5136151423 

4.875638 \ 

2.7729511992 

5.293192 

3.0736854139 

5.786112 

3.4287875496 

6.375087 | 

3.8531827169 

7.047038 

4.3374610450 

7.759080 

4.8507150955 

8.451551 

5.3499221183 \ 

9.061041 

5.7893464878 

9.533431 

6.1299460466 

9.836925 

6.3487773166 

9.975085 

6.4483977760 

9.999865 

6.4662654998 

10.00000 

6.4663630857 
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From Table 2.1, one can also notice that in the case of circular arc or the 
elliptic arc (when r\ = 2), the corresponding weight (for the sector angle equal to 
90°) is the same as cos(90°/2.0) which has been discussed in circular/elliptic 
arc section. 

Transform the surface entities to the NURBS representation 
The algorithms of transforming various surface definitions to NURBS 
surface representation follows. There are four different cases discussed in this 
section. 

Transform Cubic Parametric Spline Surface (Entity 114) to NURBS Surface 

The cubic parametric spline surface defined by IGES is composed of M 
by N cubic patches as illustrated in Figure 2.10. 


u=TU(l) 
v=TV(N+lh 


u=TU(M+l) 

v-TV(N+l) 



u=TU(l) 1 M / V 
v=TV(l) 


Figure 2.10 Illustration of cubic parametric surface defined in IGES. 


The definition of this surface is expressed as: 

S(u,v) = a + bs + cs 2 + ds 3 + 1 (e +fs + gs 2 +hs 3 ) + 

t 2 (k + Is + ms 2 + ns 3 ) + t 3 (o +ps + qs 2 + rs 3 ) 


( 2 . 11 ) 
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Two breakpoint vectors are TU(i),...,TU(M+l) and TV(i),...,TV(N+l) where 
TU(i) — u = TU(i+l) i= 1 , ... , M and s - u - TU(i) and TV(i) ^ u ^ TV(i+l) 
i= 1,.., N and t = v - TV(i). 

The strategy for transforming this entity to a BSpline tensor product 
surface is similar to the one for the cubic parametric spline. The matrix form 
for the parametric cubic spline surface, according to equation (2.11), can be ex- 
pressed in a matrix form as shown in equation (2.12). 


S(u ,v) 


[1 s s 2 s 3 ] 


a e k o' 
b f l p 
c g m q 
d h n r 


1 


t 

t 2 

t 3 


( 2 . 12 ) 


While the matrix form of the Bezier surface with Bezier control points fyj can 
be expressed as equation (2.13). 


S(u ,v) = [1km 2 k 3 1 Aij Bij Qj 


rn 


V 



(2.13) 


where 


■ 1 

0 

0 

O' 


'1 

-3 

3 

- 1 ' 

-3 

3 

0 

0 

it 

C? 

0 

3 

-6 

3 

3 

-6 

3 

0 

0 

0 

3 

-3 

-1 

3 

-3 

1 


1 

0 

0 

1 


The coefficients of this cubic parametric surface are given from an IGES 
file, therefore, the variables of equation (2.12) are all known, and the only un- 
known for equation (2.13) is matrix term of Bezier control points By. Hence, 
the Bezier control points for each bi-cubic patch are obtained by setting equa- 
tion (2.12) be equal to equation (2.13) and solving the matrix equation (2.14) 
with necessary re-parameterization. 
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where 


a eh 2 kh 2 oh\ 



bh x fh x h 2 lh x h\ ph x h\ 
ch\ gh\h 2 mh\h\ qh\h\ 


Q 



dh\ hh x h 2 nh\h\ rh\h\ 


L 


mi 

'3 0 0 O' 


'9 9 9 9* 

3 10 0 

n 1 

0 3 6 9 

3 2 10 

Q-q 

0 0 3 9 

3 3 3 3 


0 0 0 9 


(2.14) 


hi =TU(i+l)-TU(i) 
h 2 =TV(j+l)-TV(j) 


After all Bezier control patches B u are obtained, one can join each sub 
patch to form the final Bspline surface by setting the multiplicity of the knot 
value at the joint place to 3 in both directions (I and J). The advantage for 
using this algorithm is that there is no approximation or interpolation work 
involved. Therefore, the final Bspline surface represents the same geometry as 
the original parametric cubic spline surface defined in IGES the file. An IGES 
file generated by a CAD/CAM package is manipulated to demonstrate this al- 
gorithm. This IGES file, representing the nacelle of an engine, contains one 
cubic parametric surface with 34 by 24 sub patches. The final transformed 
BSpline surface ends up with 103 by 73 control points as shown in Figure 2.11. 
Similar to the parametric cubic spline curve, this resulting cubic parametric 
surface is only C° continuous (the composition algorithm is simply the inverse 
of the splitting algorithm, and in the splitting algorithm, knot insertion is re- 
peated until the multiplicity of the knot value is equal to the degree of the 
curve. The continuity is then defined as Qdegree - multiplicity ) According to this 
algorithm, the M x N bi-cubic patches will result in (3M+1) by (32V+1) control 
points of a single BSpline surface. One can perform the knot removal algo- 
rithm to reduce the redundant knots [Ref 5,38,76]. 



Figure 2.11 BSpline surface represents a nacelle of an engine. 


Transform Ruled Surface (Entity 118) to NURBS Surface 

As aforementioned, using NURBS to represent a ruled surface has been 
studied in many papers [Ref 42,43,44]. However, the definition for a ruled 
surface in IGES format is more general. The IGES defines the ruled surface 
(entity 118) as one surface formed by moving a line connecting points of either 
equal relative arc length or equal relative paranietric value (the issue of rela- 
tive arc length and relative parametric value is discussed in Chapter 4) on two 
curves and by defining a variable named “DIRFLG” as direction flag to deter- 
mine how the surface will be connected. If DIRFLG is 0, then the surface will 
connect the points on the same direction of the two curves. If the DIRFLG is 1, 
the direction of one of the curves will be reversed. The algorithm is described 
as follows: 

Suppose the two boundary curves have been converted to NURBS form 
according to the previous algorithms. Then, the first step is to make the two 
knot vectors have the same range by shifting and normalization of all the knot 
values without changing the multiplicity of any knot value. Next, use the de- 
gree raising algorithm to raise the low degree of the curve, which will yield a 
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new knot vector and new control polygon. If this new knot vector differs from 
the other knot vector, then perform knot insertion [Ref 8] by inserting the knot 
value which is not contained in both knot vectors to merge them into one final 
knot vector. For the J direction, set the order equal to 2 and the knot vector as 
[0,0,1,1]. If DIRFLG = 0, then the weights and the control net are set as W(iJ), 
d(ij) = ( weights and control points on the two boundary curve, i=0..m and 
j-0,1). If DIRFLG = 1, then the weights and the control net are set as W(i,0) , 
d(i,0) = {weights and control points of the first curve, i=0..m), and W(i,l) , 
d(i,l) = ( weights and control points of the second curve, i-m, m-1, ... , 1, 0). 
Figure 2.12 shows the NURBS rule surfaces with the same boundaries but 
with different direction flag (DIRFLG = 0 and DIRFLG= 1 , respectively). 


m 



DIRFLG=0 



DIRFLG=1 


Figure 2.12 NURBS surfaces represent the ruled surfaces. 


Transform Surface of Revolution (Entity 120) to NURBS Surface 

Again, the surface of revolution has been discussed in many places [Ref 
42,43,44,47]. However, the transforming procedures which can cover the more 
general case is presented. IGES defines the surface of revolution (entity 120) 
as the surface which is formed by rotating a boundary curve (called genera- 
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trix) with respect to a straight line (axis of revolution) from a starting angle 
(not necessarily zero) to an ending angle. This general algorithm can be 
stated as follows: 

First step is translating/rotating the axis of revolution so that it is coin- 
cident with the Z axis. It is assumed that the generatrix is defined as NURBS 
curve with the control polygon do to dm, order=k, weights wo to w m . Next, for 
each control point di (on the generatrix i = construct the surface control 

net dy j= 0 .. n at the i-th cross section according to the starting and ending 
angle by utilizing the circular arc algorithm. Based on the procedure de- 
scribed in the section of “Transform Circular Arc (entity 100) to NURBS 
curve”, n is determined by sector angle (equal to the difference between ending 
and starting angle). For example, if the angle is less than 90°, n is equal to 2, 
if the angle is in the range of 90° ~180°, n is equal to 4 ... etc. For the section 
angle 0 , the weights are set as Wy = wt, wi cos(Qln), W{, wicos(Q/n) ,....( repeat 
i vy WiCos(Q/n) with total n+1 terms). The knot vector in direction I (s) is the 
same as the one of the generatrix while the other one in direction J (t) is deter- 
mined according to the procedure described in the section of “Transform Circu- 
lar Arc (entity 100) to NURBS curve”. The control net and the weights are 
then transferred back to the original coordinates by reversing the translating 
/rotating operations. Figure 2.13 shows the construction of the associated con- 
trol polygon at each cross section for the case of n equal to 2 (section angle 
equal to 90°). 

The final NURBS definition for the constructed surface in Figure 2.13 
contains dy i=0,..,m, j=0,..,2 as the control net. The order and knot vector in 
the dir ection I are simply those of the generatrix, while the order in J direction 
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will be set as 3 and the knot vector is set as (0, 0, 0, 1, 1, 1). Weights are Wjj = 
(wi, WiCos(90/2), wi)i=0,..,m for all j. 



Figure 2.13 Illustration of surface of revolution by NURBS construction. 


Figure 2.14 illustrates an example for this algorithm. This figure dis- 
plays the “candle stand” NURBS control nets as well as the revolved surfaces 
for different starting and ending angles. Note the left figure also shows the 
generatrix. 
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Transform Tabulated Cylinder (Entity 122) to NURBS Surface 

IGES defines the Tabulated Cylinder (entity 122), or extruded surface, 
as a surface formed by moving a line segment (called generatrix) parallel to 
itself along a curve (called directrix). In other words, given a NURBS curve, 
one can generate another curve by extruding the given curve with a distance a 
along a vector V. Then the control polygon of this new curve is determined by 
Hi = di + aV. This NURBS tabulated cylinder is then defined as equation 
(2.15): 

m 1 

c, .. t-OJ-0 

S<s ’ l> (2.15) 

2 £ wtij)N^ s >mt) 

i = 0j = 0 

The knot vector in I direction is the same at the knot vector of the given curve, 
while the knot vector in J direction is set to (0., 0., 1.0, 1.0). The order kl is the 
same as the order of given curve, while k2 is set to 2. The weights W(iJ) are set 
to the weights of the two curves, and the control net is set to d(i,0)=di and 
d(i,l)=Hi for i=0 , ..., m. Figure 2.15 shows this algorithm for a tabulated cylin- 
der. 



Figure 2.15 NURBS surface represents the tabulated cylinder. 


Remark 
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In this chapter, the formulas of NURBS are first discussed. The trans- 
forming algorithms for converting different curve and surface entities to 
NURBS are presented. In the IGES standard, there are several other entities 
which are not discussed in this chapter. The reason of that is many of the enti- 
ties are not related to curve and surface definitions. For example, the entity 
106 (copious data entity), entity 116 (point entity) and entity 132 (connect 
point entity) will be treated as points, and the entity 124 (transformation ma- 
trix entity) is used for geometric transformation. And the other reason is the 
entities are already defined as NURBS, such as entity 126 (NURBS curve) 
and entity 128 (NURBS surface). Some entities will be discussed in later 
chapters. For example, the entity 130 (offset curve) and entity 140 (offset sur- 
face) will be discussed in Chapter Three. Two entities — entity 142 (curve on 
a parametric surface entity), entity 144 (trimmed parametric surface) will be 
discussed in Chapter Seven. 




CHAPTER III 

NURBS IN STATIC GRID GENERATION 


Utilizing the NURBS definitions for generating grids is the focus of this 
chapter. The organization of this chapter is as follows. First, the basic concept 
of surface grid generation is introduced. Utilizing the “parametric” property of 
NURBS entity to generate different topologies for the same surface structure 
is mentioned in the next section. The third section covers the surface genera- 
tion functions using NURBS for constructing the geometry which is commonly 
used in CFS analysis. The generating of volume grids by various NURBS op- 
tions will be presented in the last section. 

Basic Concept for Surface Grid Generation 
Structured surface grid generation involves the establishment of a one- 
to-one correspondence between the non-uniformly distributed physical sur- 
face and the uniformly distributed computational plane. Let r = ( xi(s,t ), 
X2(s,t), x 3 (s,t)) denote a parametric surface with Euclidean coordinates (xj, x 3 , 
X3) having parametric values (s,t). The surface grid associated with the com- 
putational domain { {g 1 , g 2 ) | will be denoted by r = ((xi)ij, 

( X2kj , (x 3 )ij), and Also there exists a one to one correspon- 

dence between the physical domain and the parametric distribution mesh. 
This one to one correspondence also exists between the surface distribution 
mesh and the computational domain. These relations are illustrated in Figure 
3 . 1 . 
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Physical space 


Distribution mesh Computation domain 


Figure 3.1 Physical space, Distribution mesh and Computation domain. 


Parametric Grid Generation 

From the NURBS formula listed in Chapter Two, it is easy to under- 
stand that the NURBS definition is presented in a parametric form. This 
means that for any parametric value t ( s, t for the surface or s, t, u for the 
volume) in the parametric domain, there exists' a unique point (excludes the 
singularity case) on the physical curve (or surface). This property actually can 
be applied to the other geometric definitions, such as the Bezier entity [Ref 
23,24,25], cubic splines [Ref 23,34] and so on. However, the NURBS represen- 
tation provides another unique property — the local control property which 
allows the designer to easily manipulate the geometry without altering the en- 
tire geometric definition. Utilizing these two properties, the NURBS can pro- 
vide great flexibility in grid generation applications. For example, the gener- 
ating of surface grids often requires different topologies — such as O type, C 
type and H type grid [Ref 72,73]. Generating a 3D surface from these different 
topologies for the same surface may require repeated applications of a routine. 
However, this goal can be achieved by using a 3D NURBS control polygon 
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(control net for the surface or control volume for volume case) and weights 
while generating the distribution mesh in 2D parametric space. This proce- 
dure is also referred as “griding in parametric space” [Ref 93]. After generat- 
ing the two dimensional distribution mesh tflz, 1 ,^ 2 )), one can evaluate 

the NURBS entity by utilizing these 2D grid points as the parametric values 
to achieve the final NURBS entity. The advantage of doing this is that it is 
always easier to generate the grid in 2D parametric space within the [0,1] 
square than as a complex 3D surface. Even though this approach may result 
in an unexpected grid lines in physical NURBS entity, yet this problem can be 
overcome by the “re-parameterization” algorithm which will be described in 
next chapter. Different topologies of NURBS surface grids generated by using 
the same NURBS control points and weights are demonstrated in Figure (3.2) 
~ (3.4). 




Figure 3.2 A O-type NURBS surface grid and the parametric values. 



Figure 3.3 A H—type NURBS surface grid and the parametric values. 
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Figure 3.4 A C— type NURBS surface grid and the parametric values. 


Here, the change in topologies of physical surface grid orientation is ac- 
complished by changing the topologies of the distribution mesh. The same 
strategy can also be applied to the unstructured and hybrid cases. Figure (3.5) 
demonstrates the unstructured approach, and Figure (3.6) shows two 3D sur- 
faces created from the hybrid approach. 



Figure 3.5 A unstructured NURBS surface grid and the parametric values. 



Figure 3.6 Two NURBS hybrid surface grids. 
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Surface Grid Generation bv NURBS Control Net 
The transforming procedures described in Chapter Two are used to 
transform an existing geometric definition from one form to another. For ex- 
ample, one can transform an existing cubic parametric surface to a NURBS 
definition. However, the generating procedures described in this chapter need 
to create the geometry definition from scratch. Two NURBS surface genera- 
tion functions which are frequently used for creating the grid for CFD analysis 
are presented. 


Transfinite Interpolation for NURBS Surface 

From the examples shown in Figures 3.2 to 3.5, one can understand 
that different topology grids can be generated by different types of parametric 
values. This goal can also be achieved by constructing different types of con- 
trol polygons and evaluating with regular parametric meshes. Taking the Fig- 
ure 3.3 as an example, in order to generate the H type grid, the intuitive way 
is to create the NURBS H type control net for the NURBS surface. In many of 
the numerical grid generation applications, the H type grid can be generated 
easily by the Jrans/inite interpolation algorithm (TFI) [Ref 15,73]. As a mat- 
ter of fact, the TFI algorithm is the most frequently used function for the grid 
generation. 

TFI, also referred to as Coons-Gordon patch, is a bivariate interpola- 
tion constructed from the superposition of a set of univariate interpolation 
schemes by the formation of the Boolean sum projector [Ref 73]. In other 
words, given a set of boundaries (or isoparametric curves), the TFI is a func- 
tion which constructs the interior surface grid bounded by the given bound- 
aries. The Boolean sum operator for a surface is defined in equation (3.1). 

PS = P K ®P n = P % (S) + P^S) - P^(S) (3.1) 
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where Pj-OS) interpolates the £ direction of boundaries (the given isoparametric 
curves) and P^(S) interpolates the ri direction of boundaries, while P^P^iS) 
captures the failures of P^(S) and Pr\(S). The final surface PS bi-directionally 
interpolates the given boundaries / isoparametric curves and forms interior 
surface grid points bounded with the given curves. There are many functions 
which can be applied to TFI. For example, one can use the linear, quadratic or 
even cubic interpolation for function P in equation (3.1). Taking the linear in- 
terpolation for a surface with the resolution N by M as an example, the equa- 
tion (3.1) can be re-written as equation (3.2). 

R t j = (1 - s-J) R tj + s t j R nj + (1 - t t J) R t l + t t j R i>M - (3.2) 

((1 - Sy)(l - t i J)R 11 + (1 - s i j)t i jR 1M + Sy(l - tjj)R N1 + 

Variables Rij in equation (3.2) are the control vertices which need to be deter- 
mined. For the NURBS case, the Rij could be dx, dy, dz (control points) and wij 
(the weights). 

This TFI function is a fundamental tool for generating grids in many 
grid applications. However, equations (3.1) and (3.2) can not be applied to 
NURBS TFI directly. The reason for this can be explained as follows: when 
four NURBS curves are given to generate a NURBS TFI surface, the interior 
control points can be created according to equation (3.2) (for the bi-linear in- 
terpolation) by supplying the control vertices of the boundaries without a 
problem. The problem comes when determining the interior weights. The 
addition and subtraction operations in equation (3.2) may lead to the interior 
weights being negative or zero values. Any negative weights will destroy the 
convex hull properly of a NURBS entity, while any zero weights will make the 
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control vertices lose their influence. This obstacle can be overcome by the mo- 
dified NURBS TFI [Ref 37]. The formula for this is shown in equation (3.3). 


WP(S) 


W ( wJfJ.S) + P,(S) - PjP„(S)) 

w 

ij 

W { W, 


Each term of the equation (3.3) (for the case of linear interpolation of P) is de- 
fined as follows: the P^(S) represents a NURBS “ruled surface” (described in 
previous chapter) with weights of formed in % direction (hence, the order is 
2, knot vector is [0, 0, 1, 11 in £ direction), P^(S) represents another NURBS 
“ruled surface” with weights of formed in rj direction ( order is 2, knot vector 

is [0, 0, 1, 1] in r| direction) and the P^P^(S) is a NURBS surface which is 
constructed by using the 4 comers points as the control net with orders 2 by 2 
in £ and rj directions. This is demonstrated in Figure 3.7 and 3.8. 



Figure 3.7 Four NURBS boundaries to form a TFI Surface. 


After creating the intermediate surfaces of P^(S), P^fS) and PfP^S), 
one has to perform the “knot insertion” [Ref 8] and “degree elevation” [Ref 14] 
algorithms to these three surfaces to ensure all of them have the same orders 
and same knot vectors in both the £ and T| directions. If the NURBS surfaces 
have the same orders and same knot vectors, then the dimension of control net 
would be the same also. Therefore, the control net of the final NURB TFI sur- 
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face can be obtained by adding the control net of P^(S), P^(S) and subtracting 
those of PtPx^S), while the weights are determined by multiplying and W^. 






Figure 3.8 A superposition diagram for a NURB TFI surface. 

Comparing the NURBS TFI with the traditional TFI shows that the 
NURBS TFI needs more computation because the weights need to be handled 
properly. In addition, the knot insertion and degree elevation algorithms need 
extra computation. However, this function is fundamental and useful when 
there is a need to create H type grids. Also, when generating the volume grids 
for a nozzle, this function is particularly useful to create the inlet and outlet 
surfaces. Figure 3.9 demonstrates this example. 



Figure 3.9 NURBS TFI creates the inlet / outlet surface for a nozzle. 
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Cascading Technique for NURBS Surface 
As it is discussed in the transforming procedures (described in previous 
chapter), the surface of revolution algorithm can be used to model the symmet- 
ric surfaces. In CFD applications, some of the geometries for analysis are 
symmetric objects. An example of this could be the simulation of the flow 
passing around a missile. Generally speaking, the surface of revolution algo- 
rithm can be used to model a “simple” symmetric surface, but for many of the 
CFD applications, the real geometric objects interact with other objects and 
can not be modelled by rotating a boundary curve to form a surface of revolu- 
tion. The example for this is shown in Figure 3.10 for a surface blade with the 
boundary intersected with a fin. Even though this surface is still symmetric, 
the surface of revolution ( SOR ) algorithm fails to model it. 



Figure 3.10 Example of symmetric surfaces couldn’t be modelled by SOR. 

This situation also happens to the “cascade” surface. A cascade surface 
is usually referred to the “blade-to-blade” surface in turbomachinery [Ref 56]. 
Even though most of the cascade surfaces are axis-symmetric, they can not be 
modeled by the NURBS surface revolution algorithm. Also, in the grid gen- 
eration area, creating the surface grid for the cascade is a challenge. The diffi- 
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culty of modeling the surface grids for a 3D cascade surface is that when the 
blade leading edge (or trailing edge) circle radius is too big, such as the ones in 
a turbine, or if the blade setting angle (the blade angle) is very low, it is hard to 
generate a well behaved H type grid. Particularly, there is often a grid cross- 
ing near the leading edge (or the trailing edge) for such a geometry. Tradition- 
ally, this kind of surface grid is generated by transforming the 3D surface of 
the (x, y, z) coordinate to 2D parametric (m\ a) space, griding in the parametric 
space and then transforming them back to 3D physical space according to the 
relation of the (m\ a). Detailed information can be referred to in [Ref 53,54]. 
In this research, the NURBS modeling approach is presented for modeling 
this type of geometry. 

The NURBS algorithm for modeling the cascade surface is described as 
follows: given a boundary curve of a cascade surface, transform it to a BSpline 
curve (as the curve A shown in Figure 3.11) by the interpolation technique (de- 
scribed in later section). Create a plane which bisects the surface sector angle 
(the 0, angle of aoib shown in Figure 3.11) and then vise the “Mirror” function 
[Ref 50] to reflect curve A with this plane to create the curve C. The curve C 
will have the same order , knot vector and number of control points as those of 
curve A. The next step is to create a straight line lying on the plane which 
contains the points oi,a and c. This is done by projecting the control polygons 
(described at the later section of this chapter) of curve A to this plane and set- 
ting the order , knot vector of this line to be the same as those of the curve A. 
After this line is created, perform the surface of revolution algorithm (refer to 
Chapter Two) to rotate this line with respect to the center of o^cg for a total 
sector angle of 0. A NURBS tabulated cylinder (refer to Chapter 2) with sector 
angle 0 will be generated after this step. However, this surface is not the de- 



50 

sired cascade surface. Hence, one has to perform the last step: replace the 
first and last iso-control polygons (in the axis direction) of this surface with 
the existing BSpline curves A and C. Because the tabulated cylinder is 
created by rotating a line which has the same order and knot vector as those of 
curve A, it is secure to replace the two control polygons of the surface with A 
and C without altering the entire shape of the surface. The control net, with 
curves A and C replacing the first and last iso-control polygons, is the final 
desired NURBS control net. A missile configuration, composed of the surface 
of revolution and cascade surface, is shown in Figure 3.12 to demonstrate this 
algorithm. Another example is shown in Figure 3.13 for a single rotation 
propfan modelled by this algorithm. 



Figure 3.11 Illustration of modeling cascade surface by the NURBS. 



Figure 3.12 A missile configuration modeled by the NURBS. 



Figure 3.13 A single rotation propfan modeled by the NURBS. 
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Volume Grid Generation by NURBS Control Volume 
Volume grid generation algorithms have been utilized in many CFD 
analysis procedures. A widely used technique to algebraically generate a 
three dimensional volume grid is to utilize the transfinite interpolation algo- 
rithm based on the bounding surface grids. However, the volume generation 
techniques are seldom applied to CAD/GAM applications. Even though 
NURBS representation has been widely used in many industry applications, 
the geometry modeled by the NURBS volume approaches are seldom dis- 
cussed in CAGD (Computer Aided Geometry Design) literature. In this chap- 
ter, using the NURBS volume to model the geometry for the volume grid is 
presented. Observing the curve and surface examples shown in previous 
chapter, one can realize that modeling the geometry with NURBS may only 
need very concise control polygons. Instead of storing the surface grid points, 
one can store the associated control polygon (or control net for the surface) 
with the associated weights to reduce the memory load. This is especially use- 
ful for volume grid generation. Even though the computer memory availabil- 
ity has been dramatically improved, a complicated geometry usually consumes 
a lot of computer memory for the volume grid. Storing the NURBS control net 
to reduce the size of entire volume grid is demonstrated in tha examples of this 
section. ; 

The ultimate objective of the present research effort is to explore vari- 
ous NURBS control volume options applicable to three dimensional grid gen- 
eration. In this chapter, the development of NURBS ruled volume, NURBS 
extruded volume, NURBS volume of revolution, NURBS composite volume 
and Transfinite Interpolation ( TFT) NURBS volume are discussed. 
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NURBS Control Volume for a Ruled Volume 

The easiest 3D NURBS volume to generate is the ruled NURBS vol- 
ume. The algorithm is described as follows: Given two NURBS surfaces, the 
first step to form a ruled volume is making the knot vectors of the surfaces be 
in the same range of [0~1]. Next, considering the I direction of both surfaces, 
use the degree raising technique to raise the low degree ( order - 1) of the sur- 
face. This procedure will yield a new knot vector and new control net. If the 
new knot vector differs from the other knot vector, then perform the knot in- 
sertion algorithm to merge them into one final knot vector. Then apply these 
steps to the knot vectors in J direction of both surfaces. After this step, the 
two NURBS surface will have the same orders and the same knot vectors in 
both I and J directions. This means the resolutions of control net of both sur- 
faces will be the same. Finally, connect the corresponding control point to- 
gether to form the 3D NURBS volume. The orders and knot vectors of the fi- 
nal volume in I and J directions will be the same as those of the surfaces after 
degree elevation and knot insertion, and the order in L direction will be set as 
two with knot vector set as (0,0, 1,1). Figure 3.14 shows a 3D “apple-like” 
NURBS volume and its volume grid while Figure 3.15 shows a missile configu- 
ration with the control volume formed by this algorithm. 

While defining the NURBS ruled surface, the IGES defines a variable 
named “DIRFLGT as a direction flag to control the direction of how the surface 
will be connected (refer to the previous chapter). If DIRFLG is 0, then the sur- 
face will connect the points in the same direction of the two curves, otherwise, 
the direction of one of the curves will be reversed. Similar to this definition, it 
is possible to set the two flags U DIRFLG_F and “DIRFLG_J” to control how 
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Figure 3.14 NURBS control volume and grids for ruled volume. 



Figure 3.15 NURBS control volume and grids for a missile configuration. 
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the directions of the two surfaces will be connected. The use of these flags in- 
creases the flexibility of generation options. 

NURBS Control Volume for an Extruded Volume 
The generation of a NURBS extruded volume is an extension from the 
extruded surface definition. The IGES defines the extruded surface as a sur- 
face formed by moving a line segment parallel to itself along a curve. In other 
words, given a NURBS curve, one can generate another curve by extruding 
the given curve with a distance a along a vector V. Similar to this definition, 
the NURBS extruded volume is defined as follows: Given a NURBS surface 
with control net dy and the associate weights, knot vectors and orders. The 
new surface dy can be generated by “extruding” the given surface with a dis- 
tance a along a vector V. Mathematically, this new NURBS extruded surface 
can be described as dy = dy + aV with the same orders, same weights and same 
knot vectors as those of the given surface. After this step, the algorithm of 
“ruled volume” (described in previous section, it simply connects these two 
surfaces and sets the L direction linearly) can be applied to these surfaces to 
form a final NURBS volume. Figure 3.16 shows an example of a 3D NURBS 
extruded volume. 



Figure 3.16 NURBS control volume and grids for extruded volume. 


NURBS Control Volume for a Volume of Revolution 
Another commonly used approach for generating volume grids is the 
“revolution” method. A revolution resulting in a surface is known as a “sur- 
face of revolution”, while a revolution resulting in a volume is then defined as 
a “volume of revolution”. The fact that this modelling technique can only be 
used for symmetric geometries is not limiting, as many objects in real world 
applications, such as turbomachineiy configurations, are symmetric. The ex- 
tension of a volume revolution modeled by NURBS is presented as follows: the 
definition of surface of revolution is a surface which is formed by rotating a 
given curve with respect to an arbitrary straight line from a starting angle to 
an ending angle (referred to Chapter 2). Likewise, the volume of revolution is 
defined as a volume form by rotating a given NURBS surface with respect to 
an arbitrary axis of revolution from any starting angle to an ending angle. 

The general algorithm is outlined as follows: the first step is translating 
/ rotating the axis of revolution by proper transformation matrix so that it is 
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coincident with the Z axis. This transformation matrix is also applied to the 
given NURBS surface so that the entire surface can be kept in the same posi- 
tion as the axis of rotation. It is assumed that the surface is defined (or trans- 
formed) as NURBS with the control net dy, order kl and k2, weights Wij and 
two knot vectors. The second step is to construct, for each control net djj (on 
the generatrix i= 0,..,n ), the control volume dyi 1=0 ... p at each j-th 

cross section through the starting and ending angle by utilizing the circular 
arc algorithm. In other words, this approach constructs the NURBS control 
net at each J constant plane by revolving the control polygon djj with respect 
to L direction and then “stack” them together to form a final NURBS volume. 
Figure 3.17 demonstrates this approach. The general procedure of generating 
the NURBS circular arc is described in Chapter 2. The p for the last dimen- 
sion of control volume is determined by the sector angle 0 (equal to the differ- 
ence between ending and starting angle). For example, if the angle is less 
than 90°, p is equal to 2. If the angle is in the range of 90° ~ 180°, p is equal to 
4, if in the range of 180° ~ 270°, p is 6, if it is greater than 270°, p should be 8. 
For the sector angle 0, the weights are set as (in each J constant plane, J=0,..n) 
Wijp = Wjj, Wjj cos(q/p), wuj, Wijcos(qlp ) ,... i=0,..,m ( repeat wu, wujcos(qlp) 
with total p + 1 terms). The knot vectors in directions of I (s) and J (t) are the 
same as the ones of the given surface while the knot vector in direction L ( u ) is 
determined according to the circular arc procedure. For example, when p is 2, 
the associated knot vector is set as (0, 0, 0, 1, 1, 1), for the case of p equal to 4, 
the knot vector is set as ( 0, 0, 0, 1/2, 1/2, 1, 1, 1), for the case of p equal to 6, the 
knot vector is set as (0, 0, 0, 1/3, 1/3, 2/3, 2/3, 1,1,1) and for the case of p equal 
to 8, the knot vector is (0, 0, 0, .25, .25, .5, .5, .75 , .75 ,1,1,1). Also, set the 
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orders in I and J be hi and k2 (as the orders of the original surface) while set 3 
as the order of L direction. 


rotation axis 


dmnl 

djnn2 


dy 

Figure 3.17 Illustration of constructing the NURBS volume. 

Because the NURBS has the translate and rotate invariant properties, 
the inverse transformation matrix can be applied to the control volume (with- 
out altering the weights and knot vectors) returning the volume to the original 
coordinates. Figure 3.18 shows a 3D volume grids and its control volume ac- 
cording to this algorithm. This example was developed by revolving the TFI 
surface from 0° to 180°. Because the NURBS surface TFI technique needs 
four boundary curves to define a surface, this will result in a “H” type surface 
grids. Revolving this “H” type TFI surface creates “H” type NURBS control 
volume and yields the “H” volume grids. This topology can be changed by re- 
volving another “O” type NURBS surface to form an “O” type volume grids as 
shown in Figure 3.19. Notice that the sizes of this control volume are only 







3x3x5 (for “H” type grids) and 9x9x5 (for “O” type grids), yet the resolution of 
the entire volume grid could be any number (for this case, 31 by 31 by 61). 



Figure 3.18 NURBS revolution volume for H type volume grid. 



Figure 3.19 NURBS revolution volume for G type volume grid. 


NURBS Control Volume for a Composite Volume 
A composite NURBS volume is defined as a volume consisting of lists of 
constituent volumes. The composing procedure is stated as follows: Suppose 
two constituent NURBS volumes V* and V 2 form a composite volume. Assume 
that Vi has control volume di[0:mi, 0:ni, 0:li ], weight Wi [0:mi, 0:ni, 0:li\, 
three knot vectors knot_i)i, knotj)i, knot_l)i and orders k_i)i, kj)i, k_Di 
while V 2 has control volume cfelO.vng, 0:ri2, O.'fe], weight W 2 [0:m2, Q'.ttq, O.-fe], 


60 

three knot vectors knotjfe, knotjfe, knotjfe and orders hj\, k j\, 

There are many possible combinations of the two volumes joined together. For 
example, one may join the volumes in I direction with the interface of J, L sur- 
face, or join in L direction with the interface of I, J surface, etc. Even though 
there are many cases, the procedure is similar. Take the case when joining in 
I direction as an example, the first step is to perform the degree elevation to Vi 
and V 2 so that these two volumes can have compatible degrees in I, J and L 
directions. If the two knot vectors in J direction for Vj and Vg are not the 
same, merge them together by setting the final knot vector as { knotj)i | 
knotJ )2 }, then apply the knot insertion to Vj and V 2 in J dimension. The 
same procedure should be applied to L direction if knot_l)i and knotjfe are not 
the same. After this step, Vj and will have the same degree in three direc- 
tions, and the number of control points and knot vectors in J and L directions 
will be the same. The second step is to adjust the knot vector knot_i )2 so that 
its first knot value can be the same as the last knot value knotjdi. Shifting 
the knot vector will not change the original NURBS because the basis function 
is a “normalized” basis function. The third step is to build up the final knot 
vector by joining the two knot vectors into one knot vector and setting that 
knot value at the joint point to have the multiplicity equal to (order -1). For 
example, if the knot vector knot_i)i is [0., 0., 0., 1., 1., 1.] and the knot vector 
knot_i )2 is [2., 2., 2., 3., 3., 3.], adjust the second knot vector by shifting —1 to 
each value. Thus, the knot_i)2 becomes [1,1, 1,2, 2, 2]. Suppose the final order 
of these two volumes in I direction is three, then, the final knot vector should 
be [0., 0., 0., 1., 1., 2., 2., 2.] (notice the interior knot 1 has multiplicity of (or- 
der- 1)=2). The fourth step is to match the weights at the interface surface by 
multiplying the ratio of Wi[m h j, l] / W 2 [0, j, l] to W 2 [i, j, l] for i=0, ... , m2, 
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j=0, ... , nl and 1=0, ... , 11. The last step is to construct the final control volume 
and weights by removing the d%\ 0:0 , 0:ri2, O.'fel and W 2 [0:0, 0:n2, 0:l2 ] and join- 
ing the others as one control volume and weights. 

Figure 3.20 demonstrates this algorithm. This 3D NURBS composite 
volume can be easily created by the following steps: First, use the algorithm of 
generating the ruled surface to form a NURBS surface by connecting two full 
NURBS circles, then extrude this surface with a proper distance a to form the 
first sub volume (the cylinder pipe) by the extruded algorithm. Second, ex- 
tract the last cross section surface and revolve it with 90° to form the second 
sub volume (the turning portion) by the revolving algorithm. The third step is 
exacting the last cross section surface of the second sub volume and revolving 
it 180° to form the third sub volume (the turning portion). The forth step is 
extracting the last cross section surface of the previous volume and extruding 
it with a proper distance a to form the last sub volume (the cylinder pipe) by 
the extruded algorithm, and the last step is applying the composite algorithm 
to form the final volume. The dimension of this NURBS control volume is only 
9 by 11 by 2, yet it can create volume grids as large as 121 by 91 by 31. The 
flow field solution obtained using the NPARC [Ref 16] flow simulation code for 
this configuration is also included. 
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Figure 3.20 A composite NURBS pipe volume with the NPARC solution. 
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Generally speaking, it is difficult to model a complicated geometry by a 
single NURBS control volume. However, one can construct the individual con- 
trol volume and then utilizing this composite algorithm to merge for a final 
volume. Figures 3.21 and 3.22 demonstrate the flexibility and the advantage 
of this approach. The NURBS control volumes are used to model the internal 
pipes. The griding of the turning portions of those circular pipes is a time con- 
suming and tedious task in grid generation. However, the turning portions 
can be constructed by “volume of revolution” without any difficulties. Assem- 
bling all the sub NURBS volumes makes the final single block NURBS control 
volume. 



Figure 3.21 Composite NURBS volume for a turning pipe (i). 


Figure 3.22 Composite NURBS volume for a turning pipe (ii). 

NURBS Control Volume for a TFI Volume 

Similar to the NURBS TFI surface algorithm, this approach is fre- 
quently used to generate an H type volume grid. Instead of providing 4 
NURBS curves, this algorithm requires six NURBS surfaces to generate a 
NURBS TFI volume. This algorithm is the extension from the surface to vol- 
ume, hence, the Boolean sum equation is defined as equation (3.4). 

PV = P^ = 

P{V + P n V + PtV- P^V- PyftV-PzPtV + PiPrftV (3.4) 

The P could be any interpolation function, such as the linear, quadratic 
hermit or the cubic interpolation and so on. The traditional definitions of 
each term in equation (3.4) can be found in [Ref 73]. However, as it is dis- 
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cussed in the NURBS TFI surface section, the traditional TFI approach [Ref 
61,73] can not be applied to the process of generating a NURBS TFI control 
volume because the addition and subtraction operations in equation (3.4) may 
lead to zero or negative weights in the interior control volume. Any zero 
weight will make the corresponding control point lose its influence and the 
negative weights will create undesirable grids, such as the un-bounded grids 
or crossing grids. Hence, when applying the equation (3.4) to a NURBS TFI 
volume, it is necessary to re-define the individual terms listed in equation 
(3.4). The procedure is shown as follows: suppose the six NURBS surfaces are 
all pre-defined, and the surfaces of Sj and Sfc are used for the \ direction. Sj 
and S2 have the same orders of k2, k3, and same number of control points of n 
by L (refer to Chapter Two), If the orders of these two surfaces don’t match, 
one should perform the degree raising algorithm to the low degree surface. If 
the resolutions of the control points of Si and S2 are different, then the knot 
insertion algorithm should be used to make them the same. The same proce- 
dures should be applied to the surface of S3, S4 (with the orders of kl, k3 and 
the resolutions of control net of m by L) and S5, S6 (with the orders of kl, k2 
and the resolutions of control net of m by n). After this step, the definition of 
each term for a linear NURBS TFI volume can be defined as follows: the P^V is 
a NURBS volume which is created by using the surfaces of Si and S2 with the 
algorithm of the “Ruled NURBS volume” (described in the previous section of 
this chapter). Hence, the three orders of P^V are 2, k2 and k3, while the reso- 
lutions of the control volume is 2 by n by L. The same procedures should be 
applied to P^V and P\Y. Therefore, the orders of P^V are kl, 2, k3 with the 
resolutions of the control volume of m by 2 by L, while the orders of P{V are 
kl, k2, 2 with the resolutions of the control volume of m by n by 2. P^P-^iV) is a 



66 

NURBS volume which is created by utilizing the boundaries (in £ direction) of 
Si , S2 and the comer points of S3, S4, S5, S6. In other words, it has orders of 2, 
2, k3 and the dimension of control volume of 2, 2, L. The P^P^V) and P^P^(V) 
are defined analogous — the orders of the P^P^(V) are hi, 2, 2 and the resolu- 
tion of the control volume is m, 2, 2, while the orders of the P^P^V) are 2, k2, 
2 and the resolution of the control volume is 2, n, 2. The last term of P^P^P^iV) 
is simply a NURBS control volume constructed by all the comer points of the 
six surfaces. Hence, the orders of this volume are 2, 2, 2 and the size of control 
volume is 2 by 2 by 2. These seven control volumes are illustrated in Figure 
3.23. 

After these seven intermediate control volumes are created, equation 
(3.5) should be used for the final linear NURBS TFI. This equation will avoid 
the creation of any undesired interior weights. 

In addition to the algorithm of NURBS TFI surface, one has to perform 
the “knot insertion” and “degree elevation” [Ref 8,14] algorithms to all of these 
seven intermediate control volumes to ensure all of them have the same orders 
and same knot vectors in all the Tj and C, directions, respectively. After this 
step is completed, the size of all the control volumes would be the same. 
Hence, the final control volume for NURB TFI can be obtained by adding the 
corresponding control points of P^(V), P^(V), P^(V), P^P^P^CV) and subtracting 
those of PfP^IV), P^PtfV) and P^P^CV), while the weights are determined by 
multiplication of W^, and W^. Figure 3.24 shows an H type nozzle gener- 

ated according to this approach. 
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Figure 3.23 Illustration of NURBS TFI volume. 
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Figure 3.24 A NURBS TFI nozzle with H type volume grid. 

General Interpolation Algorithm 

Another category of generating grids is by the interpolation methods. 
This technique is frequently used for modeling the sculptured discretized data 
points. Many interpolating algorithms, such as the “cubic spline interpola- 
tion”, “Lagrange interpolation”, “BSpline interpolation” and “Hermite inter- 
polation” are well known approaches. Among them, the BSpline interpolation 
technique has been frequently utilized in many of CFD researches. The exam- 
ples can be found in Yang’s General Purpose Adaptive Grid Generation System 
[Ref 83,84,85] and Shih’s TIGER (Turbomachinery Grid Generation System) 
[Ref 53,56]. Also, other grid generation softwares provide the interpolation 
option. These softwares include Genie ++ , EagleView and NGP. Many of the 
softwares use this technique only for certain order (most of them, order is set 
to four), In this study, a general BSpline interpolation algorithm which allows 
for an arbitrary order is presented. 
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The definition of interpolation for the curve data set is defined as fol- 
lows: given discrete data points Qj,j=0,..,n, construct the BSpline curve (set all 
the weights to 1 for the NURBS curve) such that the BSpline can pass through 
all the given data points Qj. Mathematically, this is described in equation 
(3.6). 


Q. = at) = C(Uj) = j = 0 , ../I (3.6) 

i = 0 

From equation (3.6), there are four unknowns: d (control polygon), k 
(the order), NHu.) (the basis function, which is related to the knot vector) and 

the parametric values Uj. The uj is used for evaluating the BSpline curve and 
can be obtained by the relative chord (arc) length parametrization or the “cen- 
tripetal” parametrization [Ref 23,34] from the given data points Qj. From the 
previous chapter, one knows that the normalized chord length parametriza- 
tion will lead uo = 0.0 and Tin-1 = 1.0. Then, the user has to decide what order 
(k) should be used (generally speaking, k=3 or 4 are sufficient). After the k 
and uj are known, the knot vector T is determined as T = [0, ... 0, tj , ... , t n _k~i , 
l,..l] with the knot value tj defined as equation (3.7). 

j+(.k — i)-i 

tj = 1 !) X U i J = 1 ’ -> n ~^ k ~ 1) (3 * 7) 

i=j 

Once the information of k (the order), knot vector and Uj are all set, solv- 
ing the equation (3.6) for the unknown d (the control points of Bspline) will 
result in a (n+ 1) by (n+1) linear equation. Notice that the knot values of 0 and 
1 have the multiplicity of k, hence, the interpolated BSpline curve will have 
the first control point identical to the first of the data set (do = Qo), and the 
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last control point identical to the last point of the data set ( dn = Q n ). This lin- 
ear equation is listed as equation (3.8). 


1 0 ... o' 
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(3.8) 


Based on the knowledge of NURBS, one knows that there are at most k 
non-zero basis functions for any particular Ui for order k. Hence, the system 
of equations in (3.8) are banded with width no larger than k. The LU decom- 
position or Gauss elimination [Ref 19] can be used to solve equation (3.8) for 
the control polygon dj. 

The interpolation of the BSpline surface for the discrete data set is 
analogous. It is defined as follows: given discrete data points Qij, i=0,..,m, 
j=0,..,n, construct a BSpline surface such that this surface can pass through 
all the given data points Qij. This can be presented as equation (3.9) 

m n 

% = S(u,v) = 2 X dfFt&iNjlty i = 0, -m J = 0, ..n (3.9) 

» = 0j=0 

The ui and Vj are the parametric values obtained by the relative chord 
(arc) length or the “centripetal” parametrization [Ref 23,34] from the given 
data points Qij. As is mentioned in previous section, the BSpline surface is 
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defined as a tensor product surface. Equation (3.9) can be rewritten to equa- 
tion (3.10) shown as follows. 


m 

Qij = Z 


i=0[j=0 




m 


|tf?0 u t ) ) = 'Zd'Nfa) 

i~0 


i-0, ..m 
j-0, ..n 


(3.10) 


Therefore, the interpolation of a BSpline surface with any arbitrary or- 
ders can be implemented by interpolating each row of the control polygon by 
using equation (3.6) to obtain the intermediate cross section control polygons. 
Then each column of these intermediate control polygons must be interpolated 
to obtain the final control net. 

A direct application of this interpolation is illustrated for the transform- 
ing algorithm for the IGES entity 130 (offset curve) and entity 140 (offset sur- 
face). 

Transform Offset Curve (Entity 130) to NURBS Curve 
An offset curve defined in IGES 5.1 is created by “offsetting” a existing 
parametric curve with a specified distance. The curve to be offset must be a 
planar curve, and must be slope continuous over the entire curve. Let this 
base curve be represented as r(t), and let V be a unit vector normal to the 
plane which contains r(t). Because r(t) is a slope continuous planar curve, its 
unit tangent vector can be obtained. Let this tangent vector be T(t), then the 
desired offset curve is defined as equation (3.11). 


OC(t ) = r(t) + f(s ) x (V x T(t)) t 1 <t<t 2 (3.11) 

Variables tj and t 2 are the starting and ending parametric values which must 
be chosen to be in the domain of the base curve r(t). The f(s) is a function con- 
trolled by a variable, FLAG, given in the IGES file. If the FLAG is equal to 1, 
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the equation (3.11) is a uniform offset curve. In other words, f(s) is set as a 
constant value, say Dj. The Di is used as the constant offset distance for the 
offset curve. If the FLAG is equal to 2, then the offset distance is set varying 
linearly from the distance of Di to D 2 by the equation (3.12). 

. . x N (s - TD.) 

f(s) -D 1 + ( D 2 - Dj) x (r£>2 _ rz>i) (3.12) 

s, TDi and TD 2 are either the accumulated arc length or the parametric values 
defined in the IGES document. 

After knowing how the offset curve OC(t) is created, it is easy to convert 
this offset curve to a BSpline representation. The first step is transforming 
the base curve to a NURBS representation according to the algorithms de- 
scribed in chapter 2. It is well known [Ref 3,4,31,33] that the OC(t) can not be 
generated by applying equation (3.11) to the NURBS control polygons of the 
base curve. The simplest reason for this is that most of the control polygons 
are not slope continuous. Also, the order of the curve is not preserved after the 
offset procedure. Hence, one should evaluate the base curve for the set of dis- 
crete data points with reasonable resolution. Furthermore, the tangent value 
should be stored at the vector T(t) at each data point. Then, equation (3.12) 
should be used obtain another set of points according to the variable FLAG (to 
have a uniform / non-uniform offset curve). The curve interpolation algorithm 
is then applied to these data points to obtain a BSpline curve with the desired 
order. Figure 3.25 shows the base curve (a NURBS curve with the control 
polygon) and a non-uniform offset curve obtained by this algorithm. 
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Figure 3.25 NURBS base curve and the non-uniform offset curve. 

Transform Offset Surface (Entity 140) to NURBS Surface 
The transforming procedure of this entity is analogous to the curve pro- 
cedure. Given a base surface S(u,v) which is defined by its specification pa- 
rameters and oriented by N(u,v), which is a differentiable field of unit normal 
vectors defined on the whole surface, the offset surface is then defined as 
equation (3.13), where D is a fixed nonzero real number used for the offset dis- 
tance. 


OS(u,v) = S(u,v ) 4- D x N(u,v ) u x < u < u 2 < v < v 2 (3.13) 

N(u,v) is obtained by the cross product of the two surface derivatives. In other 
words, N(u,v) - QS/9u x dS/dv) / ||3S/9u x 3S/9v||. Because this unit normal vec- 
tor could have two directions, the IGES documentation also provides an “offset 
indicator” [Ref 35] to avoid confusion connecting the orientation of the base 
surface. 

Transforming the offset surface to a BSpline surface is similar to the 
one of offset curve. The first step is to transform the base surface to a NURBS 
representation according to the algorithms described in Chapter Two. Again, 
it is not possible to apply the equation (3.13) to the control net of the base sur- 
face to form the offset surface. The correct procedure is to evaluate the base 
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surface for a set of data points with reasonable resolutions. Also, at each data 
point, the two surface derivatives should be calculated to obtain the union nor- 
mal value N(u,v). Then utilize the equation (3.13) to obtain the final set of 
points. The surface interpolation algorithm is then applied to these data 
points to obtain a BSpline surface (NURBS surface with all weights equal to 1) 
with the desired orders. 

In these transforming procedures, one may have a question: Since the 
equations (3.11) and (3.13) can not be applied directly to the control polygon 
(or control net) of the base curve (surface), why is there a need to transform 
the base curve (or base surface) to the NURBS representation? The answer to 
this question is related to the accuracy of the tangent vector obtained from the 
base curve (or the union normal vectors from the base surface). Equations 
(3.11) and (3.13) require T(t) and N(u,v), which are the tangent vectors of the 
curve and the union normal vectors of the surface. Once the base curve and 
the base surface are transformed to the NURBS representation analytically 
(based on the transforming algorithms described on chapter 2), the variables 
of T(t) and N(u,v) can be calculated exactly. Using the exact derivative can 
reduce the interpolation error, and thus, increase the accuracy of the final off- 
set curve / surface. 

Geometry Modeling by Interpolation Technique 

Besides the NURBS generation functions presented in previous sec- 
tions, the interpolation technique is also a widely used approach for creating 
the curves or surfaces. It has been commonly used for transforming a discrete 
data set to a BSpline (the special case of NURBS) definition. Specifically, it 
has been used to fit the ffee-form or sculptured curves and surfaces. Unlike 
the approximation method, which does not require the final curves or surfaces 
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pass through the given data set, the interpolation method does require the 
curve or surface which must satisfy the given data set precisely. 

From the equations (3.11) ~ (3.13), one can easily understand that there 
is a drawback for this technique — the disadvantage of this method is that the 
size of the control polygon of the interpolated BSpline curve will be the same 
as the size of the input data set. This destroys the advantage of saving the 
computer memory when using the NURBS definition for geometry modeling, 
when the size of data points are huge, this interpolated approach will result in 
a huge control polygon (or control net). Fortunately, the data reduction meth- 
ods [Ref 5,30], which have been discussed in many literatures, can be used to 
reduce the redundant control points and maintain the desired accuracy. Fig- 
ure 3.26 shows the BSpline curves which model an engine profile. The dis- 
crete data points are originally obtained from a digitizer. The curve interpola- 
tion method is then applied to these points to generate the BSpline control 
polygons. The associated control polygons have been reduced by the data re- 
duction algorithm described in [Ref 10,38]. After defining these curves with 
BSpline representation, the surfaces for this axial symmetric engine can be 
generated by a simple rotation with respect to the X axis. Figure 3.27 and 3.28 
show the associated surface grids and the NURBS control nets for this engine. 



Figure 3.26 The curve interpolation technique for a engine profile. 



76 



Figure 3.27 Surface grids for the multiple-duct engine. 



Figure 3.28 NURBS control nets for multiple-duct engine. 

Modeling the volume grid for this geometry requires the multi-block 
strategy. A valid domain decomposition is to break the entire surface control 
nets into several NURBS TFI sub-patches and then apply the volume of revo- 
lution technique to these sub-patches to construct the entire multi-block 
NURBS control volumes. Figures 3.29 and 3.30 show the NURBS control vol- 
ume and the volume grids. 
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Figure 3.29 Multi-blocks NURBS control volume for mock engine. 



Figure 3.30 Multi-blocks volume grids for mock engine. 
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Another geometry modeled by this technique is the wind tunnel design. 
The boundary of this wind tunnel is given with a set of discrete data points 
(121 points). The curve is interpolated by the equations (3.6) ~ (3.8), and the 
data reduction algorithm is used to reduce the redundant control polygon. 
Figure 3.31 shows the interpolated BSpline curve with the reduced control 
polygon. 



Figure 3.31 An interpolated BSpline curve for nozzle boundary. 

Because this curve is parallel to z axis, one can utilize the “surface of 
revolution” algorithm (described in chapter two) to construct a NURBS sur- 
face by revolving the curve with respect to the z axis. The surface will be sym- 
metric with respect to the z axis. One can apply the “scaling” algorithm to this 
NURBS surface with zero x, y scaling factor to generate another NURBS sur- 
face. This surface will be a degenerate surface which aligns to the z axis. This 
is demonstrated in Figure 3.32. After these two NURBS surfaces are ob- 
tained, one can perform the “ruled volume” algorithms (described in previous 
section) to construct the final NURBS 3D nozzle volume grid. Figure 3.33 
shows the final nozzle. 
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Figure 3.32 A nozzle geometry constructed by NURBS. 



Figure 3.33 A 3D NURBS nozzle. 
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Projection and Inversion 

After presenting the interpolation technique for modelling the free- 
form curve and surface, one may be interested to know how “well” the interpo- 
lated curve (or surface) describes the discrete data set. The algorithms to 
measure the deviation between the interpolated curve (or surface) and the 
sculptured data points are the “projection” and “inversion” algorithms. The 
definition of the projection problem can be stated as: Given an arbitrary point 
P = (x, y, z) in 3D space, find a point on a NURBS curve or surface such that 
the distance between this point and point P can be the shortest. The definition 
of the inversion problem is that given the point P, assuming this point is on the 
NURBS entity (curve or surface), find the associated parametric values t (for 
curve) or s, t (for surface) such that one can utilize this parametric value to 
obtain the point P. More precisely, finding t such that C(t) -P (for curve) or 
finding (s,t) such that S(s,t) = P (for surface). 

For the projection problem, if the point P is on the NURBS curve (or 
surface), then the shortest point would be the point P itself. Then, finding the 
parametric value for this P would become the inversion problem. Hence, the 
projection and inversion problems are the two sides of the same question. This 
question can then be combined with the two definitions and stated as: Given a 
point P (not necessarily lying on the NURBS entity), find the parametric value 
such that the parametric value can generate a point P*, which is on the 
NURBS entity, and such that the distance of PP* is also the shortest. If the 
point P is on the NURBS entity, then the P* would be identical to point P (or 
the distance of PP* would be zero). This is illustrated in Figure 3.34 for the 


curve case. 
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Figure 3.34 Illustration of projection and inversion for NURBS curve. 

The algorithm is described as follows (refer to Figure 3.34 for the curve 
case): Let P be the point in 3D space, C(t) be the point on the NURBS curve 
with parametric value t, and C’(t) be the tangent vector associated with para- 
metric value t. The value of t is obtained by solving the equation (5.1) formed 
by the inner product of the two vectors. 

fit) = C\t) • iCit) - P) = 0 (3.14) 

When the distance from point P to C(u) is minimized, these two vectors will be 
orthogonal to each other, and hence, the inner production function f(t) will be 
zero. The desired parametric value t which causes the shortest distance can 
be obtained by using the Newton Raphson iteration. The iterative equation is 
listed as equation (3.15). 

fiti) _ C%) • iCitf) - P) 

i+1 * f%) 1 • iCitf) - P) + \C%)\ 2 (3,15) 

The derivative of the NURBS curve, C’(tj) and C”(ti), can be found in 
[Ref 44,47]. The tolerance e is set to check the convergence, and the stopping 
criteria can be set by either | fy +1 -fy | e or | C(tj) - P | ^ e. Upon conver- 
gence, ti is the desired parametric value (for the inversion problem), and Cfc) 
is the projected point (for the projection problem). 
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The problems of projection and inversion for the NURBS surface are 
analogous. Because there are two independent variables (s, t) for a surface, it 
is necessary to create three vectors to form the equations. This yields the two 
inner production equations shown as (3.16) and (3.17). 

f(s,t) = GSM) - P) • S a (s,t) = 0 (3.16) 

g(s,t) = (S(s,t) - P ) • S t (s,t) = 0 (3.17) 

The S s (s,t) is the derivative of a NURBS surface with respect to I direc- 
tion evaluated at parametric value (s,t), while S t (s,t) is the derivative of a 
NURBS surface with respect to J direction evaluated at parametric value (s,t) 
[Ref 34,44]. When the distance between the point P to the point S(s,t) is the 
shortest, equations (3.16) and (3.17) will be satisfied. The Newton iteration 
for solving equations (3.16) and (3.17) is shown in equation (3.18). 


gs(.Si,tj) gtis^tj) 

Let vector R be R=S(s,t) - P. Then the elements of matrix on the left hand side 
in equation (3.18) can be rewritten as: 


fs(Si, tj)= I <S s | 2 + P(S; , tj) • S ss 

(3.19) 

ft(Si , tj) = S s -S t + P(s. , tj) • S st 

(3.20) 

gs(Si , tj) — S s ’ S t + R(s i , tj) • S st 

(3.21) 

gtei , ^ = |S ( | 2 + R(si , tj) ■ s u 

(3.22) 


Also A tj = tj+i - tj, A si = Si+ 1 - Si The convergence is achieved when the differ- 
ences of A tj and As* both are under the tolerance. The values of <* and s* are 
then the desired parametric values (for the inversion problem), and the point 
S(si ,ti) is the projected point. 



-AM,> 

- g(s if tj ) 


(3.18) 
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The algorithms of projection and inversion can be applied to many ap- 
plications. For example, one can use the projection algorithm to measure the 
deviation between the discrete data set and the BSpline curve (or surface) ob- 
tained by the interpolation technique. This is done by projecting the discrete 
pomts to the interpolated curve. If the corresponding projected points are 
wrthrn the tolerance with respect to the original discrete points, then the in- 
terpolated BSpline curve is good enough to model the data set. The other ex- 
ample, as shown in Figure 3.35, projects an arbitrary curve on a NURBS sur- 
face. In this Figure, two points are selected on the NURBS surface. The 
inversion algorithm is then used to obtain the associated parametric values in 

the domain of the surface. A line joining these respective parametric values is 
evaluated to form the curve on this surface. 



Figure 3.35 The projection of a 


curve to a NURBS surface. 




CHAPTER IV 

GRID REDISTRIBUTION AND EVALUATION 

NURBS has been used in the CAD/CAM system for decades. It has 
been used for free— form surface representation and modeling for a long time. 
And from the generation and manipulation abilities described in previous 
chapters, one can realize that NURBS has become a CAD/CAM industry stan- 
dard due to its versatile properties. However, there are several difficulties 
which inhibit people from fully utilizing this representation. The software 
tools available to designers for creating and reshaping such geometry are 
often inefficient for industrial parts, thereby consuming a large amount of de- 
signer time. 


Obstacles of Using NURBS 

The obstacles of using NURBS are discussed as follows: 

Geometry Fidelity Maintenance 

The NURBS representation is a parametric representation, and it is difficult 
to obtain the desired parametric value to maintain the critical geometric loca- 
tion. An example of this would be a NURBS representation which contains a 
sharp comer, which is defined as a slope discontinuity. How can the proper 
parametric value be obtained to maintain the original geometric information? 
Figure 4.1 shows the control polygon and its grid points on a NURBS compos- 
ite curve. The curve does not maintain the sharp comer at the joint portion 
of the straight line and circular arc. Hence, the geometric fidelity is lost. The 
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challenge of this problem is to obtain the proper parametric value t such that 
the final curve would be the one shown in Figure 4.2. 



Figure 4.1 Improper parametric values lose the geometry fidelity. 



Figure 4.2 Proper parametric values keep the geometry fidelity. 

Distribution Control of the Grid Points on a 3D Physical Space 
The NURBS entity (curve, surface, or volume) is presented as a para- 
metric format, and the grid point on a NURBS entity is generated by evaluat- 
ing the parametric value t (or s, t for surface, s, t, u for the volume). However, 
the designer desires the distribution of the grid points on the physical NURBS 
entity, not the distribution of the parametric values. For example, evenly dis- 
tributed parametric values t may not result in a sequence of evenly distrib- 
uted grid points of C(t) on the physical NURBS curve. Finding the correct 
parametric values to obtain the desired distribution on 3D physical space has 
presented a problem to engineers for a long time. Figure 4.3 shows the control 
polygon and grid points on a NURBS curve. Even though the parametric val- 
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ues used for evaluation of this curve are distributed evenly in parametric do- 
main, one can observe that the grid points are actually packed toward the 
middle. The location of the control polygon, the weights and even the knot vec- 
tor are all possible factors in controlling the final NURBS entity. The chal- 
lenge of this problem is, without altering the NURBS definition (control poly- 
gon, weights and knot vector), calculating the proper parametric values to 
obtain the desired distribution on the physical NURBS entity as shown in Fig- 
ure 4.4. 



Figure 4.3 Illustration of undesirable distribution on NURBS curve. 



Figure 4.4 Illustration of desirable distribution on NURBS curve. 

Bad Parameterization 

It has been previously mentioned that, setting the communication be- 
tween CAD/CAM and grid generation system is a ideal for shortening the en- 
tire CFS process, because most the geometry of interest is created from the 
CAD/CAM package. However, when importing the geometry from a CAD/ 
CAM package, the engineer has to handle this geometry with all information 
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pre-defined. Unfortunately, that predefined information may lead to unex- 
pected surface grids. Taking the Figure 4.5 as an example, this surface is ob- 
tained from the package CATIA (a product of the IBM CAD system). The sur- 
face is presented as a NURBS surface. From the Figure 4.5, one can see that 
the control polygons are clustered in the middle of the surface in the J direc- 
tion as well as the lower index of the I direction of the surface. The knot vec- 
tors are set evenly in both I and J directions. These NURBS definitions (both 
the location of control net and the setting of the knot vectors) make the surface 
grid fines packing towards the clustered control nets. This unexpected and 
undesired packing of the surface grid fines will produce a poor quality grid 
which might result in a divergence of CFS simulation. The situation is re- 
ferred to as “bad parameterization”. This is the first penalty one has to pay for 
the geometry communication between systems due to an unexperienced CAD/ 
CAM designer. Also, the improper geometry definitions happen frequently 
during this communication. Since the engineers from grid generation did not 
participate in designing the geometry when it was in CAD system, the chal- 
lenge becomes finding a solution to obtaining the desired distribution without 
altering the geometry shape as shown in Figure 4.6. 

Another penalty is related to geometry manipulation. One of the 
NURBS properties is the local control property which states that altering the 
location of the control points or the values of weight will only modify the shape 
of geometry locally. This local control property is very attractive to the design- 
er because the modification will not change the entire geometry. However, im- 
proper manipulation of the NURBS control points or weights may lead to poor 
grid points (fines) packing, and this situation also results in “bad parameter- 
ization”. Maintaining the desired shape of geometry (after changing the con- 
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trol points or weights) and the proper geometric parameterization is challeng- 
ing. Figure 4.7 shows the control net and the surface grid lines of a NURBS 
surface after increasing the weight of one of the control points. One can see 
that, the surface grid lines are packed toward the control points with in- 
creased weight. The challenge of this problem is finding a better parameter- 
ization without altering the desired shape of the geometry as shown in Figure 
4.8. 



Figure 4.5 Bad locations of a NURBS control net. 



Figure 4.6 A uniform distribution NURBS surface. 
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Figure 4.7 Improper manipulation leading to bad parameterization. 



Figure 4.8 Improved parameterization surface after re-parameterization. 

Computation Intensive for NURBS Evaluation 
The other complaint, which is not related to geometry communication, 
for using the NURBS modeling is the computationally intense evaluation. 
Given the control polygon d and weights W from equations (2.1)~(2.4), the 
challenge for efficient evaluation is finding the basis function Nf(t) at particu- 
lar value t as quickly as possible, since this basis function is defined recursive- 
ly. Of the many approaches to evaluating the grid points on a NURBS entity, 


90 

the “de Boor” [Ref 18] method is the most commonly used. Even though many 
computer languages, such as C and Pascal, allow the recursive call, evaluating 
this basis function with an internal recursive call needs substantial memory 
and computation time. Also, when the order(s) or the dimension of the 
NURBS control points is high, especially for the case of a 3D NURBS control 
volume, the evaluation will take significant computation time. Avoiding the 
reclusive call and quickly evaluating the grid points on a NURBS entity is an 
important issue. 

Strategies for Overcoming NURBS Difficulties 
The strategies for overcoming the problems associated with using 
NURBS are presented in this section. 

As one can understand from equations (2.1) ~ (2.4), the NURBS is de- 
fined in a parametric format. Thus, any parametric value will create a physi- 
cal grid point on a NURBS entity. Taking the NURBS curve as an example, 
any value t in parametric space, which must be inside the domain of the knot 
vector, will result in a point C(t) in the physical space. In other words, the 
evaluation domain for any parametric value t must be located in the range of 
the knot vector of T, t €E [3\_j , ... , T n+ i ] . The first problem — maintaining the 
NURBS geometric fidelity, is related to the knot vector. 

Maintaining the Sharp Corner with NURBS 
In order to know how to maintain the fidelity of the geometry, it is nec- 
essary to know the relationship between the NURBS entity and the knot vec- 
tors. This is illustrated as follows: a NURBS curve has order equal to 3 and 
control points a, b, c, d and e, as shown in Figure 4.9. Its associated knot vec- 
tor is [0., 0., 0., 0.3333, 0.6667, 1., 1., 1.]. The points ml and m2 are the points 



evaluated with the parametric values t equal to 0.3333 and 0.6667 (notice that 
these two values are exactly the same as the interior knot values of the knot 
vector). The arc ami is created from those parametric values located in the 
first knot span [0., 0.3333], and the arc mlm2 is created from the parametric 
values located in the second knot span [0.3333, 0.6667], while the last arc m2e 
is created from the parametric values located in the last knot span [0.6667, 
1.0]. In each knot span, the associated NURBS arc is a smooth curve. If the 
knot value ti = 0.3333 increases, and the knot value #2 = 0.6667 decreases (but 
ti then the arc mlm2 is “shrunk”. Before the tj and #2 values collide, the 
segments are all C 2 continuously differentiable due to the order equal to 3. In 
other words, in each arc segment, there is no discontinuous point. However, as 
one can imagine, when the and #2 value collide, the arc mlm2 becomes a 
singularity point, and the discontinuous problem may arise. It is well known 
[Ref 23,34,44] that the continuity of a NURBS entity at knot value t is equal to 
( k-l-m ) where k is the order and m is the multiplicity of the knot value. 
Hence, if t± and #2 collide, then the multiplicity of this knot value becomes 2, 
and the continuity is = C°, which defines the sharp comer (discontinu- 

ous point). This situation occurs when the two NURBS entities are joined to- 
gether by the “composite” algorithm described in chapter 2. In that algorithm, 
the multiplicity of the knot value has been set to ( order-1 ) in order to main- 
tain the discontinuous portion. 

Knowing the reason for the creation of the sharp comer, the strategy for 
maintaining the discontinuous point is as follows: when the sequence of para- 
metric values ti are generated, one should first check the domain of the knot 
vector (for this example, the domain of knot vector [0.0, 0.3333, 0.6667, 1.0] ). 
If any knot value T\ has multiplicity ni which gives C k-l-m ) equal to 0, then 
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one should chose the parametric value tj with the absolute minimum differ- 
ence of (tj - Ti ) and assign 2J to tj. For example, if the knot vector for the curve 
shown in Figure 4.1 is set [0., 0., 0., 0.25, 0.52, 0.52, 0.75* 0.75, 1., 1., 1.0], and 
the user would like to distribute 11 points on the curve with even spacing in 
parametric space (ti = O.li i=0, ..., 10), since T 5 = T 6 =0.52 (with multiplicity 
equal to 2), then #5 = 0.5 must be adjusted to the value of T 5 - Tq =0.52. The 
same procedure should apply to the parametric value #7 = 0.7 (or ts = 0.8, since 
| ty -Tj\ = | ts -T 7 1 ), updating it to the value of T7 = Tg =0.75fe After this pro- 
cess, the parametric sequence is ensured to maintain all possible discontinu- 
ous points. 



Figure 4.9 The relationship between knot value and the NURBS curve. 

It is necessary to add a “physical” checking algorithm to find out the 
“real” sharp comers after the previous procedures. In many cases, the knot 
value which has multiplicity greater than 1 does not intend to represent the 
discontinuous point. For example, the circular arc shown in Figure 2.2 has the 
knot vector set [0., 0., 0., 0.5, 0.5, 1., 1., 1.]. Even though the knot value 0.5 
has the multiplicity 2, according to the previous discussion, it should have a 
sharp comer because Ct 3-1 '® = C°. However, one knows there are no discon- 
tinuous points in a semi-circle. There is no harm in setting any parametric 
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values to the knot value 0.5, but to avoid taking the point M as the sharp cor- 
ner, one should add a distinguishing algorithm described as follows: If ti is the 
suspect parametric value which might create the discontinuous point, then 
setting three additional parametric values ti+j and ti+ 2 , with t^i slightly 
smaller than ti, and ti+j, ti +2 slightly greater than ti (“slightly” can be defined 
as the distance of ti and %.i to be 1.0e-7) will discourage discontinuity. Eva- 
luating with these four parametric values will yield four points — C 
C(ti), C (ti+i) and C (ti+ 2 ) in the physical curve entity. These four points will 
form two angles, 0i and 02, as shown in Figure 4.10. If the absolute difference 
of these two angles is greater them the tolerance (for example, 5°), then this 
parametric ti is indeed a value which can create a real physical sharp corner. 



Re-Parameterization Algorithm 

The second and third barriers can be overcome by the re-parameteriza- 
tion algorithm. This algorithm allows the user to determine the parametric 
values which yield the desired distribution on physical NURBS entities. Two 
approaches are introduced: one is performed with an iterative method, and the 
other is implemented by linear interpolation. 
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Before discussing the algorithms, it is necessary to define several nota- 
tions. For a NURBS curve with resolution ni, there are several ID arrays 
which need to be defined: 

( 1 ) csi (i), i=l, ..., ni be the parametric values associated with the desired dis- 
tribution of the curve in physical space; 

(2) cs 2 (i), i= 1 , ...» ni be the normalized chord length of the curve with desired 
distribution (for example, cs 2 (i) = (i-1) / (ni- 1 ) for even distribution); 

(3) C$ 3 (i), i= 1, ..., ni be the normalized chord length of the curve evaluated at 
parametric values csjd), i=l , ..., ni. 

The explanation of these ID arrays are referred to Figure 4.3 and 4.4. 
Consider the example shown in Figure 4.3. If the designer would like to have 
the final curve as shown in Figure 4.4, then cs 2 (i) will be a ID array which 
contains the distribution packing evenly, and csi(i) are the parametric values 
which are to be determined, such that the cs 3 (i), the normalized chord length 
of final curve, would be the same as cs 2 (i) (or | cs 2 (i) - cs 3 (i) | be minimized for 
all i=l,.., ni). 

For the NURBS tensor product surface with resolution ni by nj, the 2D 
arrays are defined as follows: 

( 1 ) (ssi(ij), sti(ij)), i=l, ni,j= 1, ..., nj be the parametric distribution mesh 
associated with the desired distribution of the surface in physical space; 

(2) (ss 2 (ij), st 2 (ij)), i= 1 , ..., ni,j= 1 , ..., nj be the normalized chord length of the 
surface with desired distribution in direction I and J respectively; 

(3) (ss 3 (ij), st 3 (ij)), i= 1, ..., ni,j=l , ..., nj be the normalized chord length of the 
surface evaluated at parametric values (ssidj), stidj)) i=l, ..., ni,j= 1 , ..., nj. 

Consider Figure 4.7 and 4.8 as examples. If the designer would like to 
have the final surface, as shown in Figure 4.7, then (ss 2 dj), st2dj)) would be a 
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2D array which contains the even distribution, and (ssi(ij), stiGJ)) would be 
the parametric values which are to be determined such that (SS3GJ), st3(ij)), 
the normalized chord length of fined surface, would be the same as (ss2(ij), 
st2(ij)) or within certain tolerance. 

For the NURBS tensor product volume with resolution ni, nj and nl, the 
3D arrays are defined as follows: 

(1) (vsi(ij,l), vti(ij,l), vui(ij,l)), z=l, ..., ni, j= 1, ..., nj and 1=1, ..., nl be the 
parametric distribution volume associated with the desired distribution of the 
volume in physical space; 

(2) (vs2(ij,l), vt2(ij,l), vu2(ij,l)), i=l, ..., ni, j= 1, ..., nj and 1=1, ..., nl be the 
normalized chord length of the volume with desired distribution in direction I, 
J and L respectively; 

(3) (vs3(ij,l), vt3(ij,l), vu3(ij,l)), i= 1, ..., ni, j= 1, ..., nj and 1 = 1 , ..., nl be the 
normalized chord length of the volume evaluated at parametric values 
(vsi(ij,l), vti(ij,l), vui(ij,l)) i= 1, ..., ni,j= 1 , ..., nj and 1=1 , ..., nl. 

The illustration of these 3D arrays are analogous to curve and surface 

cases. 

1> Iteration for Re-Parameterization Algorithm. 

The procedure for the reparameterization algorithm for a NURBS curve 
is stated as follows: The first step is to initialize the values of csi(i) by setting 
csi(i) = cs2(i). Once CS2G) are determined, the second step is to evaluate the 
NURBS curve using the parametric values of csi(i) and calculate the normal- 
ized chord length, CS3, of this curve. Generally, CS3G) will not be the same as 
cs2(i). If any value of CS3G), say CS3Q), is greater (or smaller) than cs2(I) by the 
tolerance, because CS3Q) is obtained by evaluating the NURBS curve with the 
parametric value of csi(I, then the value of csj(l ) should be decreased (or in- 
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creased if it is smaller). This process needs to be checked with all the values 
of cs2(I) and cs3(I) iteratively until all the differences of | cs2(i) - cs3(i) | are all 
within the tolerance for all i=l,...,ni. More precisely, the iterative algorithm 
can be described in the computer pseudocode shown as Algorithm 4 . 1 . 

Algorithm 4.1 

for (j=0; j<max_iteration ; j++) { 

check whether \cs2(i)-cs3(i)\ < tolerance, 0<i<ni 
if (yea) break the for loop and program stop; 

else 

{ 

for (i=l;i<ni;i++) { 
if(cs2(i) < cs3(i)) 

csi(i) += (cs2(i) - CS3 (i)) *( csi (V-csj (i- 1 )); 

else 

csi(i) += (cs 2 (i) - cs3(i))*(cs 1 (i+l)-cs 1 (i)); 

} 

j Repeat Step 2. 

) 

The maxjterotion and tolerance are set according to different request. If the 
tolerance is set too small, the program needs more iterations. For the exam- 
ples shown in Figure 4.3 and 4.4, the tolerance is set as 1.0e-6, and it needs 51 
iterations. 

The reparameterization algorithm for the NURBS surface is analogous 
to the curve algorithm. The first step is to initialize the values of (ssi(ij), 
sti(ij)) as the desired distribution mesh (ss2(ij),st2(ij)). The second step is to 
evaluate the surface with these parametric values to obtain (ss3(ij), st3(ij)). 
Then the normalized chord length of the entire surface is calculated for 
(ss3(ij), st3(ij)). Again, (ss3(ij), st3(ij)) will be different than the desired dis- 
tribution mesh of (ss2(ij), st2(ij)). The idea of adjusting the parametric mesh 
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is simply the extension of ID arrays of the curve to 2D arrays and is presented 
in the computer pseudocode shown as Algorithm 4.2. 

Algorithm 4.2 

for ( k-0 k<max_iteration; k++) { 

check whether | ss 2 dj)-ss 3 dj) j < tolerance for 0<i<ni, 0<j<nj 
if (no) { 

for (j=OJ<njJ++) 
for (i=l;i<ni-l;i++) { 
if(ss 2 (ij) < ss 3 (ij)) 

ssjdJ) += (ss 2 (ij)-ss 3 (ij)) *( ssj (ij)-ssi (i-lj)); 

else 

ssidj) += (ss 2 (ij)-ss 3 (ij))*(ssi (i+lJ)-ssj (ij)); 

} 

} 

check whether \st 2 (ij)-st 3 (ij)\ < tolerance for 0<i<ni, 0<j<nj 
if (no) { 

for ( i=0;i<ni;i++ ) 
for (j=lj<nj-lj++) { 

if (st 2 (ij) < st 3 (ij)) 

stidj) += (st 2 dJ)-st3dJ))*(st 1 dJ)-st 1 dJ-V); 

else 

stxdj) += ( st 2 dj)-st 3 dj)) *(stt dJ+D-stj dj)); 

> 1 

i/(both checking return yes) 

break the loop and program stop. 
j else Repeat Step 2. 

From this pseudocode, one can realize that adjusting the parametric 
value in the I direction for ssidj) will effect the calculating of stjdJ), and ad- 
justing the value for stjdJ) will effect the calculating of ssidj) too. This situa- 
tion is easy to understand because the surface is defined as a tensor product 
surface. However, due to this mutual effect, the convergence of the entire pro- 
cedure needs more computation time that that of a curve. For the examples 
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shown in Figures 4.7 and 4.8, the tolerance is set as 1.0e-5, and it needs 216 
iterations for the 31 by 31 surface resolution. If the tolerance is set to 1.0e-6, 
the number of iterations increases to 1500 iterations for the same resolution. 

The volume reparameterization algorithm for NURBS can be extended 
from 2D surface arrays to 3D volume arrays. However, since the NURBS vol- 
ume is also defined as a tensor product format, adjusting the parametric val- 
ues of any one of the directions (/, J or K) will automatically effect the calcula- 
tion of the other two directions. This situation makes it difficult for the 
algorithm to reach the convergence of a small tolerance, and it also greatly in- 
creases the computation time for one iteration of volume evaluation. Due to 
this drawback, the iterative algorithm is seldom utilized for the NURBS vol- 
ume case. An alternative way for this algorithm is the linear interpolation 
approach described in the next section. 

2> Linear Interpolation for Re-Parameterization Algorithm. 

For the curve case, in order to obtain cs 3 (i), the initial values of csj(i) 
must be initialized. One can set csi(i) as evenly distributed initially, evaluate 
the curve with these parametric values to obtain the final curve, then normal- 
ize the entire chord length to obtain cs 3 (i). Generally, the CS 3 W will be differ- 
ent than cs 2 (i). How does one adjust csi(i) to have new cs 3 (i) which is close to 
cs 2 (ip. The entire procedure is described in the Algorithm 4.3. 

Algorithm 4.3 

for (i=l; i<(ni-l); i++) { 

find the location of j such that CS 3 Q) cs 2 (i) < csq(j+ 1 ). 

let a = cs 2(0 - css 3 ( 7 ); P = cs^fj+l) - cs 2 (i); 

} cs(i) = (cs 1 (j+l)*a + csj (/)*p) / (a + p); 
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The theory behind this pseudocode is the simple “linear-interpolation”. 
The final cs(i) should replace csj(i) as the final desired parametric distribu- 
tion. Using this distribution to evaluate the NURBS curve shown in Figure 
4.3 will yield the result shown in Figure 4.4 — which has the even distribution 
(point packing) on the physical curve. The linear interpolation may not be 
very accurate, especially when the resolution is very low. Therefore, if the fi- 
nal curve does not achieve the desired distribution, one should repeat this pro- 
cedure again. According to research experience, the result of this algorithm is 
very satisfactory after 3 iterations. 

For the surface case, the procedure is similar to the NURBS curve. One 
can set (ssj (ij), stidj)) as evenly distributed initially, then evaluate the sur- 
face with these parametric values to obtain (ss 3 (ij), st 3 (ij)). Most likely 
(ss 3 (ij), st 3 (ij)) will be different than (ss 2 dj), st 2 dj))- The procedure to adjust 
(ssidj), stidj)) so that (sssdj), stsdj)) can be as close as to (ss 2 dj), st 2 dj)) is 
presented in Algorithm 4.4 as a computer pseudocode: 

Algorithm, 4.4 

for (j=l; j<(nj—l); j++) 

for (i=l; i<(ni- 1); i++) { 

search the index of //Jsuch that (ss 2 dj), st 2 dj)) is located within the 
cell of (ss 3 ff, J),st 3 (I,J)), (ss 3 ff+l,J),st 3 ff+l,J)), (ss 3 (l J+l),st 3 (I, J+1)) 
and (SS 3 (J+ 1 , J+ 1 ) ,st 3 (7+1, J+1)) 

let (ss 2 dj)>st 2 dj)) = (l-aXl-PXssafl, J),st 3 ff, J)) + (l-a)(3fss3(7, J+1), 
st 3 (I f J +l))+a( 1— p)(ss 3 (1+1, J), sf 3 (I+l, J))+afi(ss 3 (1+1, J+1), st 3 (I+l, 
J+ 1 )) and solve for a and |5 ; 

new (ssdj),stdj)) = (l-a)(l-$)(ssi(I, J), stiff, J)) + (l-a)fXss/(7, J+1), 
stiff, J+ 1 )) + <x(l— §)(ssiff+l, J), stiff+1, J)) + a$(ssiff+l, J+1), stiff + 1 , 

} J+1) >> 

One can verify that the bi-linear interpolation is used in this algorithm 
for finding the new parametric values. Even this bi— linear interpolation is 
more accurate than linear interpolation, but when the resolutions are not high 
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enough, two or three iterations may be needed to improve the quality. After 
this algorithm, the final ((ss(ij), st(ij)) should update ((ssjGJ), sti(ij)) and 
then ((ssi(ij), stiGJ)) is the desired distribution mesh in which to generate the 
desired surface grids shown in Figure 4.8. In this algorithm, the Newton 
method [Ref 191 is used for solving the non-linear equation for multi-vari- 
ables a and p. 

After knowing how the algorithms work for curve and surface, the algo- 
rithm for the volume is analogous. The only thing different is that the trili- 
near interpolation should be used for solving the equation (4.1) by Newton it- 
eration for variables a, p and y. 

(vs 2 GJ,l), vt 2 vu 2 (ij,l))=( l-a)( l-p)(. l-y)( vs 3 (I,J,L), vt 3 (I,J,L), vu 3 (I,J,L)) 

+a (l-fi)(l-y)(vs 3 (I+1,J,L), vt 3 (I+l,J,L),vu 3 (I+1,J,L)) 
+(l-a)f}(l-y)(vs 3 (I,J+l,L), vt 3 (I,J+1,L), vu 3 (I t J +1,L)) 

+(l-a) ( 1 -fi)y(vs 3 (I,J,L+1), vt 3 (I,J,L+1), vu 3 (I,J,L+l)) 
+af}(l-y)(vs 3 (I+l,J+l,L),vt 3 (I+l,J+l,L),vu 3 (I+l,J+l,L)) (4-1) 

+a(l -fi)y(vs 3 (I+1,J,L+1), vt 3 (I+l,J,L+l), vu 3 (I+1,J,L+1)) 

+(l-a){fy(\ vs 3 (I,J+1,L+1), vt 3 (I,J+1,L+1), vu 3 (I,J+1,L+1)) 

+a/fy(vs 3 (I+1,J +1,L+1), vt 3 (I+1,J +1,L+1), vu 3 (I+l,J +1,L+1)) 

After a, P and y are obtained, the new parametric values are determined as 
shown in equation (4.2). 

(vs(ij,l), vt(ij,l), vu(ij,l))=( 1 -a)( l-fi)(l-y)(vsi (I,J,L), vtj (I,J,L), vui (I,J,L)) 

+a(l-fi)( l-y)(vsi (I+1,J,L), vti (I+1,J,L), vui (I+1,J,L)) 

+(l-a)f}(l-y)(vsi (I,J+l,L),vti (I,J+l,L),vui (I,J+1,L)) 

+(l-a) ( 1 -fi)y( vs i (I,J,L+1), vti (I,J,L+1), vuj_ (I,J,L+1)) ^ ^ 

+a(3(l -y)(vsi (I+1,J+1,L), vti (I+1,J+1,L), vui (I+1,J+1,L)) 

+a(l -fi)y( vsj (I+1,J,L+1), vti (I+1,J,L+1), vuj(I +1,J,L+1)) 

+(l-a)/3y(vsi (I,J+1,L+1), vti(I,J+l,L+l), vui (I,J+1,L+1)) 

+a(fy(vsi (I+l,J+l,L+l),vti (I+l,J+l,L+l),vui(I+l,J+l,L+l)) 



101 

This approach, compared to the volume case of the iterative algorithm, is 
much less computationally expensive. 

Figure 4.11 and 4.12 show the results for a volume grid re-parameter- 
ization algorithm. 



Figure 4.11 A volume grid with undesirable packing. 



Figure 4.12 A volume grid with desirable packing after re-parameterization. 
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Singularity Control 

As one examines these re-parameterization algorithms, one can easily 
find out that there is a serious flaw which may lead these algorithms to fail. 
This flaw is related to the “singularity” problem. Taking a NURBS curve as an 
example, if all the points on this curve collapse to one point, then this NURBS 
curve is a singular curve. The same definition can be applied to NURBS sur- 
face and volume. If any of the iso-parametric line on a NURBS surface col- 
lapses to one point, then this line is defined as singular lines. For a volume, if 
any of the iso-plane on a NURBS volume collapses to a line, then this plane is 
defined as singular plane. When these singularity problems occur, the total 
arc length from calculating it according to the re-parameterization algorithms 
will be zero. Since the normalized arc lengths are obtained by dividing the 
total arc length to the individual ones, this will lead to the operation of 0/0 
which is mathematically undefined and can not be implemented by any com- 
puter language. These singularity problems happen in many cases. For ex- 
ample, the surface grid which represents the canopy of a aircraft has a singu- 
lar line at the nose position; a surface grid which represents the missile has a 
singular line in the nose position; and the volume grid of a cylinder (or any 
cylindrical pipe) has a singular plan at the axial direction. From these exam- 
ples, one notices that the singularity problems occured because the control 
points collapse to one point (for the surface case) or one line (for the volume 
case). While evaluating the NURBS entities with certain parametric values 
by utilizing these collapsed control points, the singularity problem arises. 
Hence, it is necessary to enhance the algorithms to handle this problems. 

The strategy of the enhancement algorithm is related to the machine 
accuracy (or called machine precision). The machine accuracy, commonly rep- 
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resented as symbol e, is defined as the smallest positive real number such that 
1 + e > e. On the Silicon Graphics Personal Iris, this number is equal to 1.0~ 16 
(double precision). In many numerical simulations, this number is needed to 
represent the finite precision arithmetic of the computer architecture. For ex- 
ample, the stopping criteria of an iteration scheme is dependent upon the ma- 
chine accuracy. A variable expected to be zero in numerical representation 
may not be reached due to the finite precision of the computer memory repre- 
sentation. Hence, in many numerical applications, the exact zero is replaced 
by a value related to e, say if a variable is less than Js, then this variable is 
assumed to be equal to zero. This concept is also utilized to avoid the singular- 
ity problems. Taking the example of the NURBS surface with a singular line 
shown in Figure 4.13, the grid line evaluated with the parametric values of 
(ssi(0J), sti(0J) j=0, ..., nj will shrink to one singular point due to the control 
vertices doj collapsing to one point. However, if one perturbes these paramet- 
ric values by a small values, say Je, and then re-evaluates the surface, the 
returned grid line will not be the same as the singular one since these para- 
metric values are no longer exact zero. Instead, it will return a grid line with 
a small but recognizable total arc length. Even though the total arc length is 
small, the normalization process will make the values of (ss 3 ( 0 J), st 3 ( 0 J) 
j=0,...,nj to 0.0 ~ 1.0 and will avoid the uncertain situation of 0/0. This proce- 
dure is shown as the Algorithm 4.5. 

Algorithm 4.5 

for (j=0; j<nj; j++) { 

ssi(0j) += e; ssi(ni-lj) -= e; 

} 

for (i=0; i<ni; i++) { 

stitt, 0) += e; stitt, nj-1) -= e; 


} 
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Figure 4.14 illustrates the result of this algorithm applied to a NURBS surface 
with a singular line. 



Figure 4.13 A NURBS surface with a singular line. 



Figure 4.14 Reparameterization for a NURBS surface with singularity. 

The associated approach can be applied to 3D NURBS volume. These 
“offset” procedures for the 3D parametric values can be summarized in Algo- 
rithm 4.6. Figure 4.15 shows a 3D NURBS cylindric pipe evaluated with even 
parametric values. Notice that in its L direction, the surface degenerates to a 
singer line. The result of reparameterization for this volume is shown in Fig- 
ure 4.16. 


Algorithm 4.6 


for (k=0; k<nk; k++) 
for (j=0; j<nj; j++) { 

vsi(OJ,k) += e; vsi(ni-lj,k) -= e; 

} 

for (k=0; k<nk; k++) 
for (i=0; i<ni; i++) { 

vti(i,0,k) += e; cti(i,nj-l,k) -= e; 

} 

for (j=0; j<nj; j++) 
for (i=0; i<ni; i++) { 

vui(ij,0) += e; vui(ij,nj-l) -= e; 

} 



! Sul 


Figure 4.15 A NURBS volume grid with a singular plane. 
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Efficient Evaluation for NURBS Entity 
As it was described in the beginning of this chapter, using the intrinsic 
computer recursive call for evaluating the NURBS entity is the worst choice. 
Currently, the “de Boor” algorithm [Ref 18,231, which is based on the geomet- 
ric construction, is widely utilized. The illustration of this algorithm is dem- 
onstrated with a BSpline curve, which has the same curve configuration as 
shown in the Figure 4.17. 


cU 

Figure 4.17 Illustration of the de Boor algorithm fora NURBS curve. 

Before discussing the evaluation of a NURBS curve, the BSpline curve 
should be introduced first. The definition of the BSpline curve is similar to the 
NURBS curve except the weights are all equal to 1. According to the property 
of the basis function, if all the weights are the same, the denominator term of 
the equation (2.1) will sum to 1. Hence, the BSpline curve is presented as 
equation (4.3). 

m 

C(s) = X d i N f (s) (4.3) 

i = 0 

If one would like to evaluate t = 0.5 for this curve, the first step is locat- 
ing the proper knot span such that Tj ^ t< Tj + i for the knot vector T. For this 
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case, the I would be 3 since the knot vector is set as [0, 0, 0, 0.3333, 0.6667, 1, 
1, 1]. The second step is mapping the ratios of T 2 : T 3 , T 3 : t, t : T 4 (let these 
ratios be a, (3 and y) and the ratios of T 3 : t, t : T 4 , T 4 : T 5 (p, y and a) to the 
proper control polygon shown in Figure 4.17. This mapping will yield two in- 
termediate points <2q, d\. The third step is mapping the ratio of T 3 : t and t : 
T 4 to the line of dgdj. The intermediate point created after this step is the 
point on the curve at parametric t. These procedures are demonstrated on 
Figure 4.17. 

The algorithm for obtaining all the intermediate points and the final 
point on a BSpline curve (with order equal to k) for a particular parametric 
value t can be summarized as Algorithm. 4.7. 

Algorithm 4.7 

find out the index I such that T[I] ^ t < T[I+ 1 ]. 
for (r=l; r «S(&-1); r++) 

for (i= I - (£-1); i ^ (I-r); i++) 

Jr — jjO ~~ t) rJ r ~l 4- — — ZKi — H r ~^ 

i (T[i + k]~ T[i + rl) * (T[i + k] - T[i + r]) i+1 

The evaluation of a NURBS curve by the “de Boor” method is related to 
“homogenization”. The 4D control vertices d* = [ d^w, dyW, dzW, w ] are ob- 
tained by multiplying the weights to the corresponding 3D control polygon [c^, 
dy, dg 1, and these new control vertices d* (including the weights) are called the 
4D homogenous coordinates. Performing the “de Boor” algorithm to these 4D 
control polygons for the particulate parametric value t will yield a new 4D 
point C*(t) = [ dfx Wi, diy W{, d^ wi, Wi ]. The corresponding point of C*( t ) on 
3D space is then obtained by “inhomogenizing” the 4D point and is accom- 



108 

plished by dividing the C*( t ) with the last component, W{. Hence, the 3D 
point C(t)=[ dix(t), di y (t), diz(t)]. 

The definition of the BSpline surface is similar to NURBS definition. 
The only different is the weights of the BSpline surface are all same. This 
leads the denominator term of equation (2.3) vanish [Ref 34,441. Hence, the 
formula is represented as equation (4.4). 

n m 

S(s,t) = 2 IX *?■<•> N T {t) “ (4.4) 

7 = 0 i=0 U-. 

Since the BSpline surface is also defined as a tensor product surface, the equa- 
tion (4.4) can be rewritten as equation (4.5). 


I 

0 


i=o 


m n 

X d tj Nf(s) Nf(t) = ^d* N fit) 


(4.5) 


> 0 


This equation shows that the evaluation of a BSpline surface can be ac- 
complished by applying the BSpline curve algorithm to the isoparametric 
curves in I direction to obtain the intermediate control points, and then treat 
them as the isoparametric control polygon in J direction, applying the curve 
algorithm to these control polygons to obtain the final BSpline surface. In oth- 
er words, the surface algorithm is implemented by applying two of the curve 
algorithms. 

The detailed implementation of the NURBS surface can be found in 
Farian’s [Ref 23,34,44]. Basically, one should apply the “homogenization” to 


obtain the 4D control vertices dy - [dxyWy, dy^Wy, dzyWy, Wy\, then treat 
these new control vertices as the control points of a BSpline surface and per- 
form the BSpline evaluation algorithm in 4D space for the particular paramet- 
ric values. Finally, project the values back to 3D space by “inhomogenizing”. 
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The evaluation of NURBS volumes by “de Boor” algorithm can be ex- 
tended analogously. Compared to the method using the recursive call imple- 
mented in many computer languages, this approach avoids the recursive call 
and imbeds the evaluation of the basis functions to the entire process. The 
evaluation process is much faster. Also, the geometric construction process 
helps the user to know how the points are generated. Many software packages 
adopt this algorithm. For example, the NGP (National Grid Project) [Ref 
29,71] used this approach to evaluate the NURBS entities. 

Even though the computation of the “homogenizing” and “inhomogeniz- 
ing” of the “de Boor” algorithm is not very expensive, it is an extra cost for the 
entire process. Also, the algorithm requires the iterative process repeated for 
evaluating the x, y, z and w components. Examine the equations (2.1) ~ (2.4) 
thoughtfully. One should know that after the basis functions are evaluated, 
the control points regarding the x,y, z components can share these basis func- 
tions without repeated calculation. The detail approach is described as fol- 
lows: 

Consider all the basis functions of the NURBS curve shown in Figure 
4.9. According to equation (2.2), those basis functions can be plotted as shown 
in Figure 4.18. 
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Figure 4.18 Non-zero quadratic basis functions. 

From Figure 4.18, the first basis function exists in the knot span [0.0, 
.3333], and second basis function exists in the knot span of [0.0, .6667], and so 
forth. Also, notice that in any of the knot span, the non-zero basis functions 
sxim to 1. For example, for parametric value t - .3333, the basis functions 

Nf(t) + Np) + Np) = 0.5 + 0.5 + 0.0 = 1.0. This situation also indicates 
that the number of non-zero basis functions in any of the particular knot span 
is equal to the order of the curve (3 for this case). (Even though Np) is zero at 
t = 0.3333, but as t slightly increases, this basis function becomes non-zero.) 
Due to the local control property, one can also understand that not all the con- 
trol points are involved for a particular t, only 3 (equal to the order of this 
NURBS curve) control points are used for the evaluation. For example, if the 
parametric value t is located in the knot span [0.0, .3333], the control points of 
do, di and d% are used, if it is in the knot span [.3333, .6667], the control points 
of di , cfe and da are used, if it is in the knot span [.6667, 1.0], then the control 
points of c? 2 , do and d* are used. Generally, if the parametric value t is located 
in the knot span of [Tj, Tj+jl, then the control points of dj_( or der-l)i 
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di-(orcler-l)+l, •••> dj are involved for the evaluation. This statement can also be 
applied to the NURBS surface and volume. Take the surface as an example: if 
the parametric values (s, t) are within the criteria of Till] ^ s < TtfI+1] and 
T 2 U] ^ t < T 2 &+I], then the control net involved for the evaluation is in the 
region of di+^j+j where i=I- (order _i - 1), ..., I j = J- (order J - 1), ..., J. The 
volume can be extended analogously. 

With this knowledge, the evaluation for the NURBS entity is fairly 
straightforward. If the user would like to evaluate a NURBS curve for ni 
points with parametric values ti i=0, ... ni-1, then the evaluation procedure for 
this NURBS curve (with order k and the control polygon di i=0, ..., n) can be 
presented by Algorithm 4.8. 

Algorithm 4.8 

for (i=0; i<ni; i++) { 

I = locate the proper knot span for t[i]; 
find all non-zero basis functions, B[0:k-1]; 
vari - var 2 = var 3 - var 4 =0.0; 
for <j=0; j<&; j++) { 

varj += w(J+I—(k—l )] x B[j]; 
var 2 += w(j+I-(k-l) ] x dx{j+I-(k-l)]xB{j]; 
var 3 += w[j+I-(k-l)] x dyU+I-(k-l)]yB!j]; 
var 4 += w[]+I-(k-l)] x dz[j+I-(k-l)]xB[j]; 

^ curx[i]=var 2 /vari; cury[i]=var 3 / vari; curz[i]=var 4 / vari ; 

From. Algorithm 4.8, one can understand that when the non-zero basis 
functions are obtained for the parametric value ti in the particular knot span, 
then only k {order of the curve) operations are needed to evaluate the point of 
C( ti ). Compare to the “de Boor” geometric construction, this algorithm is 
more intuitive and easy to understand based on the equation (2.1). 
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The other advantage of this approach is that it is fairly easy to extend 
from the ID curve case to higher dimensions. Since the NURBS surface and 
volume are all defined by the tensor product form, the corresponding algo- 
rithms for the surface and volume are shown as Algorithm 4.9 and 4.10, re- 
spectively. 

The NURBS information ( the order (a), the control vertices, the weights, 
and the knot vectors)) is all given in the Algorithm 4.8 ~ 4.10, the only chal- 
lenge left for a fast computation is obtaining all the non-zero- basis functions 
efficiently in the particular knot span. Indeed, the evaluation of the recursive 
defined basis functions is the key point oi Algorithms 4.8 ~ 4.10. 


Algorithm 4.9 


for (j=0; j<nj; j++) 
for (i=0; i<ni; i++) { 

/ = locate the proper knot span for s[i][j]; 

J = locate the proper knot span for t[i][j]; 

find all non— zero basis functions, Bi[0:ki~l], Bj[0:k2~l]; 

vari = var 2 = var 3 - var 4 =0.0; 

for (jj=0; jj<& 2 ; jj++) 

for (ii=0; ii <ki; ii++) { 

vari+=w[ii+I^i-l)][Jj+J-(k 2 -l)]xBi[ii]xBjljj]; 
var 2 +=w[ii+I-(ki -1) ][jj+J-(k 2 -l) ]xdx[ii+I-(ki -l)][jj+J-(k. 2 -l) ]x 
BdiiJxBjUj]; 

var 3 +=w[ii+I-(ki~l)][jj+J-(k 2 -l)]xdy[ii+I-(ki-l)][jj+J-(k 2 -l)]x 
Bi[ii]xBj[jj]; 

var4+=w[ii+I-(ki~l)][jj+J-(k2-l)]xdz[ii+I-(ki-l)][jj+J-(k2-l)]'X 


} 


} BrfiiJxBjljj]; 

surx[i][j]=var 2 / vary, sury[i][j]=var 3 / varj; surz[i][j]=var 4 / vary, 
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Algorithm 4.10 


for (1=0; 1 <nl; 1++) 
for (j=0; j<nj; j++) 
for (i=0; i <ni; i++) { 

locate the knot span I, J, L for s[i][j][l], t[i][j][l], u[i][j][l]; 

find all non-zero basis functions, Bi[0:ki-1], Bj[0:k2~l], Bi[ 0 :k 3 ~l ]; 

vari - var 2 = var 3 = var 4 =0.0; 

for (11=0; ll<fc 3 ; 11++) 

for (jj=0; jj<& 2 ; jj++) 

for (ii=0; ii <kj; ii++) { 

varj +=w[ii+I-(ki-l)]{jj+J-(k2-l)][ll+K-(k3-l)]xBi[ii]xBj[jj]xBi[ll]; 
var2+=w[ii+I-(ki-l)][jj+J-(k2-l)]xBi[ii]xBj[jrj]xBi[ll]x 

dx[ii+I-(ki-l)]ljj^^2-V][ll+L-(k3-V]; 

var3+=w[ii+I-(ki-l)][jj+J-(k2-l)]xBi[ii]xBj[jj]xBi[ll]'x 

dy[ii+I-(k 1 -l)][jj+J-(k 2 -l)][ll+L-(k 3 -l)]; 


} 


var 4 +=w[ii+I-(ki-l)]ljj+J^(k 2 -l)IxBi[ii]xBjljj]xBi[lU'x 

dz[ii+I-(k 1 -l)][jj+J-(k 2 -l)Mll+Lr-(k 3 -l)]; 

} 

surx[i][j]=var 2 / vari; sury[i][j]=var 3 / vari; surz[i][j]=var 4 / vari ; 


An efficient evaluation of the basis functions was first presented in de 
Boor’s [Ref 18]. This approach was implemented by an iterative (non-recur- 
sive) approach. This evaluation process is briefly discussed as follows: 

Based on the previous observation, one knows that for a NURBS curve 
with order equal to As, if the parametric value t is located in the knot span I 
such that Ti t < Ti+j, then there are at most k non-zero basis functions. 
Also, according to equation (2.2), one knows a basis function of k order is de- 
fined by the combination of two basis functions with k—1 order. These two 
facts conclude to the relationship, shown as Figure 4.19, of all the non-zero 
basis functions with all possible orders for the parametric value t in the knot 
span I. 
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Figure 4.19 Nonzero BSpline basis functions on [Tj, Tj+j] 


In each column of Figure 4.19, there are exactly k non-zero basis func- 
tions, and any basis function in the jth column is obtained by the combination 
of two other basis functions in the j-1 th column. Also, notice that when com- 
puted the basis functions N J . and 1 in thej'th column, the fact that one of 

their neighbors to the j-lth column is zero is utilized. For example, when cal- 
culating Nf, the Nj_ 1 is zero, and when calculating Nf_ v the Nj +1 is zero. 
In addition to this, there exists a simple relationship between all the non-zero 
basis functions in the jth. column. Take j equals 3 ( order equal to 3) as an ex- 
ample, according to equation (2.2), these non-zero basis functions can be 
listed as follows: 


JVf- 2 <0 = (t T _ T ‘ -* } JV?_ 2 (fl + r T ‘ + l t ) tf-i® 

2 i 2 i- 2 2 i + l 2 i - 1 


1 ¥?_,«> + ‘' Jtj 

2 i + 1 2 i- 1 2 i + 2 ~ 2 i 


W i+ 2 ~ * > 


Nf(t) 


(4.6) 


(4.7) 
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Nf(t) = if '‘jl N?tt) + (r i*3 * l 

■*1+2 ^1 + 3 ^i + l 


(4.8) 


The relationship is that the 


-Nf_ 1 which appears in the second 


T i+ 1 ~T i _ 1 

term of equation (4.6) also appears in the first term of equation (4.7). Similar 
results also can be applied to equations (4.7) and (4.8). Furthermore, two 
qualities are introduced and shown as equation (4.9) and (4.10). 


5Jeft(j) = t - T i+1 _j (4.9) 

h_right(j) - 7}+j - 1 (4.10) 


This simple relationship can be presented in a clearing fashion by applying 
equations (4.9) ~ (4.10) to equations (4.6) ~ (4.8). The re-formed equations 
are shown as equations (4.11) ~ (4.13). 



djeft(3) ^ djightiX)) Ar2 

d_right(0) + d_left(3) i ~ 2 d_right(l) + djeft(2) * -1 


(4.11) 


N z = 
iV i- 1 


Nf = 


djeft( 2) 


:N ? , + 


d_right(2) 


N ? 


d_right( 1) + d_left(2) i_1 d_right( 2) 4- djefti 1) * 


djeftd) 


at2 , d_right{3) ^ 

• 1 ./A\ . C 1 


d_right(2) + djeftil) 1 d_right(3 ) + djeft(0) 


i + 1 


(4.12) 

(4.13) 


Equations (4.11) ~ (4.13) provide the intrinsic understanding for an effi- 
cient evaluation of all non— zero basis functions for a parametric value t in the 
knot span of [7J, TJ+j]. This efficient evaluation process, which utilizes an it- 
erative approach, for these non-zero basis functions is listed as Algorithm 
4.11. 
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Algorithm 4.11 

basis [0] = 1.0; 

for (j=0;j< order; j++) { 

delta_l(/] = t - knot[J+l-fl; 

delta_r[/] = knot[/+/] - 1; 

save = 0.0; 

for {i-O; i<j; i++) { 

term = basis [i] / ( delta_r [i+1] + delta_l[/-i]; 
basis [i] = save + delta_r [i+1] * term; 
save = delta_l \j-i] * term; 

1 

basis (/] = save; 

} 

Currently, CAGI (Computer Aided Grid interface) applies Algorithm 
4.11 to Algorithms 4.8 ~ 4.10 for evaluating the NURBS curves, surfaces and 
volumes. 

Compared to the different evaluation approaches discussed in this sec- 
tion, the recursive algorithm needs the most computational time. The geomet- 
ric construction of “de Boor” method is currently the one used most for many 
packages [Ref 21,28,44]. However, the one adopted in CAGI is the most effi- 
cient algorithm — both in memory and computation time. The reasons for 
this is that it avoids the processes of “homogenization”, “inhomogenization”, 
and also when the non-zero basis functions are obtained, these non-zero func- 
tions can be applied to any coordinates (x, y and z) without repeating the same 
iterations. Besides these advantages, it is the easiest to extend to any higher 
dimensions. In order to compare the computation speed for different evaluat- 
ing approaches, a C program is implemented for the test. This program is 
executed in a personal IRIS machine which has “R2000A/R3000” processor 
chip, 56 Mbytes main memory and 64 Kbytes of cache. Let the algorithms 
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adopted in CAGI be denoted as algorithm A, and the “de Boor” be denoted as 
algorithm B, Tables 4.1 and 4.2 show the CPU time required for algorithms A 
and B for different NURBS curves and surfaces. 


Table 4.1 Comparison time for evaluation of NURBS curves. 


order 

(k) 

resolution 

(np) 

control point 
(n) 

CPU time for A 
(microseconds) 

CPU time for B 
(microseconds) 

4 

1000 

5 

49333 

128000 

4 

5000 

5 

248000 

650000 

6 

1000 

18 

93000 

293000 

6 

5000 

18 

464500 

1478250 


Table 4.2 Comparison time for evaluation of NURBS surfaces. 


BSSSi 

resolutions 
(ni, nj) 

control net 
(m, n) 

UFUtirnTforA 

(seconds) 

CPU time for B 
(seconds) 

4,3 

100, 100 

5,3 

0.4650 

0. 5710 

4,3 

500, 500 

5,3 

11.719 

18.509 

12, 10 

100,100 

79, 230 

3.7690 

46.398 

12, 10 

500, 500 

79, 230 

92.739 

562.068 


The command “getrusage” was used to measure the CPU time needed. 
The CPU time was obtained by repeating the same case for 3 times and then 
averaging the total time spent. 

From the results shown in Table 4.1 and 4.2, one can observe that the 
two methods are very completive when the orders, resolution and control 
points are low. However, the CPU time needed for algorithm A is much small- 
er than that of algorithm B while the orders , resolutions and control net are 
highly increased. 
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As is aforementioned, the NURBS representation is becomeing the 
standard for geometry description in CAD/CAM design systems. In addition, 
because of its properties, NURBS has been used for modeling objects in many 
areas, such as the manufacturing industry, entertainment industries, art and 
even modeling the scenes for virtual reality applications. In terms of numeri- 
cal grid generation and CFS applications, the difficulties and barriers of using 
NURBS must be solved to encourage the CFS community to utilize this model- 
ing technique. These problems discussed in this chapter are important. If the 
volume grids generated by NURBS can not maintain the geometry fidelity, 
then the CFS algorithm simply solves the different (or wrong) problems. The 
distribution control on grid points and the bad parameterization problems are 
highly related to grid quality — which has a direct influence on the flow con- 
vergence rate and solution precision. A grid with bad quality may even lead to 
the divergence of the CFS simulation results. Even though the computational 
ability has been improved for many different types of computers, a fast, robust 
and efficient evaluation for NURBS is required to facilitate the grid genera- 
tion process. This is especially true when handling high resolutions moving 
grids (or moving boundary problems) which require fast evaluation in each 
time step. 



CHAPTER V 

NURBS IN DYNAMIC GRID GENERATION 

After discussing the difficulties of utilizing NURBS and the associated 
overcoming strategies, one can fully apply the NURBS to many CFD simula- 
tions. As is discussed in previous chapters, the reparameterization algorithm 
allows the user to control the distribution of the grid points (lines) in physical 
space. This is useful for the surface (volume) grid refinement, particularly for 
the application of grid adaptation. The approach used here is different than 
that in Yang’s [Ref 83,84]. In Yang’s work, the BSpline interpolation is used 
for the iso-curve interpolation, while in this application, the entire distribu- 
tion (adaptive) mesh / volume is utilized for the new adaptive grids. The inter- 
polation involves the accumulation of numerical errors which may lead to the 
losing of the original geometry definition, especially for the viscous turbulent 
flow simulations. In addition, this interpolation is utilized for the entire sur- 
face grid by iso-curve interpolation, it will be a computational intensive work 
if the resolution of the surface grid is high. The approach used in this applica- 
tion will avoid these challenges. The computational examples are presented 
in this chapter. 

An efficient grid generation with a good quality grid is important for 
most CFS simulations. This is particularly true for the unsteady deforming 
geometry type of simulation. For this type of problem, there is a need to gen- 
erate the grid at each time step. If the grid generation is not efficient, a sig- 
nificant part of the computation time will be used for generating the grid. Also 
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if the grid quality is not satisfactory, the final results may not be accurate, or it 
may take more time for convergence. Hence, the grid generation procedure 
plays an important role in this type of simulation. NURBS is a good candidate 
for generating grids for this type of problem. The reasons for this are the 
parametric formula of NURBS can easily rebuild the geometry, and the robust 
evaluation algorithm described in the previous chapter makes the generating 
grid efficient. Many NURBS geometries are constructed by very concise con- 
trol polygons (control net or control volume). This means that comparing the 
grid sizes, the size of the NURBS control net is much smaller than that of the 
entire grids. And the NURBS local control property allows the geometry to be 
locally modified without altering the entire geometric shape. If there is a need 
to create a local deforming geometry, it is ideal to perform the deformation to 
the NURBS control points. Since the resolution of the NURBS control net is 
smaller, it does not take much computational time to apply perturbations to 
the geometry. In addition, the reparameterization algorithm provides a de- 
sired distribution on the deformed geometries to maintain the grid quality at 
different time steps. Also, the grid smoothness can be achieved by the com- 
bination of the proper NURBS orders and knot vectors. This is because the 
continuity of the NURBS entity is determined by the associated order and 
knot vector and is equal to where k is the order and m is the multiplic- 

ity of the evaluated knot value. These applications can be found in the moving 
grid example and the temporally deforming geometries. These issues are dis- 
cussed as follows: 

Grid Adaptation 

Since the re-parameterization algorithm can control the grid packing 
on the geometry in physical space, its application can then be used for the grid 
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adaptation. This is demonstrated by Figures 5.1 ~ 5.5 which show a 2D 
NURBS surface grid modeled for simulating a supersonic flow over an oval- 
type object. A NURBS control net is generated to model the 2D flow field con- 
figuration. As it was discussed in chapter 2, a circular arc can be represented 
by NURBS with very few control polygons. For a semi-circular arc with a ra- 
dius r, as shown in Figure 5.1, the 2D control points are (counterclockwise) 
(r,0), (r,r), (0,r), (-r,r) and (~r,0), while the corresponding weights are 1, 
cos(45°), 1, cos(45°) and 1. The knot values are (0, 0, 0, 0.5, 0.5, 1, 1, 1) with 
the order set to 3. In order to simulate the flow passing through the oval-type 
object, one can apply the algorithm described in Piegl’s [Ref 42,43,441 for the 
semi-ellipse. However, the easiest and simplest way to model this conic arc 
(the ellipse) is by utilizing the “scaling” operation: since this r is not fixed, 
another semi-circular arc can be created with radius r’ such that r’ < r while 
keeping the weights , knot values and the order unchanged, then “scaling” the 
entire control polygon by scaling the x andy components of each control point 
by (a/r, b/r) — yielding a new control polygon (a,0), ( a, 6 ), (0,6), (-0,6) and 
(-a,0) of the semi-ellipse, hence, the semi major is a and semi minor is 6. Af- 
ter these two NURBS curves are created, one can perform the ^NURBS Ruled” 
surface algorithm (described in Chapter Two) to generate the surface grid for 
the flow field simulation. The NURBS information of this surface is listed in 
Appendix B for reference. Supersonic flow at Mach number 3.0 is applied to 
this geometry. After the initial grid is generated, the NPARC flow simulation 
code [Ref 16] is used to create the initial solution. These examples are demon- 
strated in Figures 5.2 and 5.3. 




Figure 5.2 NURBS control net with the initial surface grid. 




Figure 5.3 Initial solution (Mach number plot). 


After obtaining the initial grid and solution, the adaptive program 
“GAGE2D” by Yang [Ref 83,85] is used for calculation of the distribution mesh 
and the weight functions. GAGE2D was modified, however, with this new ap- 
proach for full NURBS evaluation. The idea is that after the adaptive code 
calculates the adaptive mesh, instead of interpolating iso-curve line by line, 
the re— parameterization algorithm is utilized to generate a new surface grid 
which can reflect the adaptive mesh to the physical space. Figures 5.4 abd 5.5 
show the adaptive grids (with NURBS control net) and the corresponding 
solutions after the adaptive iterations. 

Similar to this 2D case, a 3D NURBS control volume is used to model a 
missile for flow simulation. The 3D NURBS control volume and the initial 
grid for this generic missile configuration are shown in Figure 5.6. This con- 
trol volume is generated by creating two NURBS surfaces - one for the missile 
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Figure 5.5 Adaptive solution (first iteration, Mach number plot). 
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configuration and the other for the outer free stream surface. The “ruled vol- 
ume” algorithm is then applied to these two NURBS surfaces for constructing 
the final control volume. Supersonic flow at Mach = 2.5 and 14 degree angle of 
attack is applied to this geometry. Because of the simulation of the angle of 
attack, there is no need to generate a symmetric grid around the body. Hence, 
the NURBS control volume should not be symmetric in the circumferential 
direction. The outer free stream NURBS surface was originall|created from a 
surface (if revolution algorithm. In order to model a non-symmetric surface, 
half of the control net must be reduced by a scaling factor. This, however, re- 
sults in the initial grid with the uneven distribution in circumferential direc- 
tion C K direction) even though the even parametric values are used. One 
should be able to observe this “bad parameterization” situation from Figure 

5.6. 

Figure 5.7 presents the flow solution obtained from NPARC [Ref 16] 
flow solver with the Baldwin-Lomax algebraic turbulence model option. This 
Figure shows a streamwise cut of the initial grid and solution. The grid size 

for this simulation is set as 121 x 89 x 89. 

Similar to the 2D adaptation case, the adaptive volume (the new para- 
metric values) is calculated according to the the initial solution. Unlike the 
approach of line by line interpolation used for obtaining the adaptive grid in 
Yang’s [Ref 83,84] and Thornburg’s [Ref 67], the 3D re-parameterization algo- 
rithm is utilized, along with this new adaptive volume, to obtain the final 
adaptive volume grid. Figure 5.8 shows this adaptive volume grid in the per- 
spective view. This re-parameterization technique provides the capability of 
precise grid distribution control and also maintains the geometry fidelity. One 
can notice that the “bad parametrization” shown in Figure 5.6 has been fixed 






127 

after the adaptive procedure. The adaptive grid and solution in streamwise 
cut are shown in Figure 5.9. Figure 5.10 presents the adaptive solution with 
the NURBS control volume in the perspective view. 


Figure 5.8 Adaptive grid in perspective view. 




Figure 5.9 Adaptive grid and solution (Mach number plot). 


Figure 5.10 Adaptive solution in perspective view. 
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Temporally Deforming Geometry Model bv NURBS 
Unlike the other geometric representation, NURBS allows the engineer 
to modify the geometry locally by changing the associated control points or the 
associated weights. This property is called the “local control” property. The 
other geometric representations, such as the Bezier and cubic spline entities, 
do not have the similar characteristics. Hence, if one changes the control 
points of the Bezier curve or the coefficients of a cubic spline surface, the en- 
tire geometry will change. This NURBS local control property attracts many 
designers because it provides the flexibility to model the geometric shape easi- 
ly. Figure 5.11 demonstrates this property. This figure shows the original 
NURBS surface and the NURBS surface with the perturbed control point. 
One can observe that after the control point is perturbed, only the local surface 
is effected. 



Figure 5.11 NURBS local control property. 

NURBS can model some physical phenomenons by utilizing this local 
control property. For instance, by properly changing the location of the control 
polygons, the NURBS can be used to model the wave propagation. As is dem- 
onstrated in previous chapters, many of the geometries modelled by NURBS 
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have a smaller size of control polygon (net or volume) than that of the 
associated grid points. Manipulating the control polygon (control net or vol- 
ume) instead of the grid points on curve (surface or volume) will save the com- 
putational time. It is reasonable to work on the control polygon and then re- 
evaluate the entire geometry with the efficient NURBS evaluation algorithm 
(described in chapter four). This strategy has been applied to CFD simula- 
tions, such as the moving boundaries problem or the temporally deforming ge- 
ometry [Ref 13]. For the unsteady problems, the grids needed to be generated 
efficiently in each time step to shorten the simulation procedure. Manipulat- 
ing the control polygon and then evaluating the grid will meet this require- 
ment. The following examples show the deforming geometry problems with 
this approach. 

NURBS Models the One Dimensional Wave Movement 
A straight line can be easily modelled by NURBS with two control 
points and set the order equal to two. However, to model the wave propaga- 
tion along a one dimensional curve, the order of the NURBS curve must be set 
to at least three (or greater than three) to get a smooth wave movement. Also, 
the approach presented in this study is applying the manipulation to the con- 
trol polygon, hence, one has to set enough numbers of control points to catch 
enough simulating information. In this one dimensional case, a NURBS curve 
with nine control points is used. The order of this curve is three. In order to 
simulate the wave propagation along this curve, the one dimensional wave 
equation [Ref 1] is solved and applied to the associated NURBS control poly- 
gon at different time steps. The one dimensional wave equation can be de- 
scribed in equation (5.1) shown as follows. 
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u t + au x = 0 (5.1) 

where u is the displacement iny axis direction, and a is the propagation speed. 
This simple wave equation can be solved by many available finite difference 
schemes [Ref 1]. However, in this study, instead of solving this equation to all 
the grid points (121 points) on the curve, the wave equation is applying to the 
NURBS control polygon. Originally, the control polygon (consisting of nine 
control points) lines in the x axis. At the first time step, the initial condition is 
set so that the y coordinate of the first control point is perturbed. The bound- 
ary condition is set to the left hand side of the curve. The y coordinate of the 
first control point is oscillating with certain damping functions at each time 
step. The explicit upwind scheme [Ref 1] is used to solve the equation. At 
each time step, when the new control points are obtained according to equa- 
tion (5.1), the NURBS curve evaluation algorithm (described in Chapter four) 
is utilized to generate the new curve. Solving equation (5.1) will result in a 
propagating wave based on the initial condition and boundary condition. 
Hence, the final curve evaluated with different control polygons will result in 
the new curve which simulates the wave propagation. This simulation is dem- 
onstrated in Figure 5.12 at different time steps. 

One may doubt why this approach is needed for this simulation, be- 
cause many other approaches can reach the same result without using 
NURBS representation. The reason is that the proof of the versatility of using 
NURBS for different applications. And for the complicated cases, such as the 
surface or the volume with large grid points, this approach can reach the goal 
with less computational time. 
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time = 0 sec 


time = 1 sec 



time = 2 sec 



time = 3.5 sec 



time = 5 sec 


NURBS control polygon. NURBS curve. 

Figure 5.12 NURBS curve simulates the wave propagation. 
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NURBS Models the Two Dimensional Wave Movement 
Similar to the one dimensional case, the NURBS surface can be used to 
simulate the two dimensional wave propagation. This is done by applying the 
two dimensional wave equation to the control net of the NURBS surface and 
then evaluating the surface grid with the new control nets at each time step. 
The resulting surface grid will be smooth if proper orders and knot vectors are 
selected for the NURBS evaluation. The 2D wave equation is described in 
equation (5.2) as follows. 

u t + au x + bu y = 0 (5.2) 

A NURBS surface with orders 3 by 3 is used to model this wave propagation. 
The knot vectors are set so that the multiplicities of knot values in the evalua- 
tion domain are all equal to one. Hence, the surface grid will be C 1 continu- 
ously. Initially, the control net (consists of 19 by 19 control points) sits on the 
xy plane. At the time step equals to 0, the initial condition is set so that none 
of the control points is perturbed. The boundary condition is set to the center 
of the control net (the indices of the (9, 9) control point). The z coordinate of 
this control point is oscillating at each time step with a designed damping 
function. The function will create different z attitudes for the central control 
point at different time steps. The z displacement verses time steps of control 
points (9,9) is shown in Figure 5.13. 
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Figure 5.13 The oscillating function of the central control point. 

An explicit finite difference scheme is used to solve equation (5.2) for 
propagating the z displacement. At each time step, when the new control net 
is obtained, the efficient evaluation algorithm (described in chapter four) is 
then used to generate the corresponding surface. In this simulation, the wave 
speeds (variables a, b in equation (5.2)) are set to be the same, and the time 
step is set to 0.1. This simulating process is demonstrated in the Figure 5.14 
for different time steps. 

The NURBS is composed of the control polygon (net or volume), orders , 
knot vectors and the weights. A previous example shows the deforming geom- 
etry obtained by manipulating the locations of the control points. Another al- 
ternative for modeling the motion grid is manipulating the weights of the 
NURBS representation. The next example of deforming geometry is demon- 
strated by manipulating not only the control points but also the associated 
weights. As is discussed in chapter two, the NURBS modeling algorithms for 
the circular arcs, elliptic arcs and the superelliptic arcs are analogous. For an 
arc (circular arc, elliptic arc or the superelliptic arc) with sector angle less 
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than 90°, three control points with proper weights will be enough to model the 
arc. The differences between these three entities are the determination of the 
weights and the values of the radii (or semi-major and semi-minor for the el- 
lipse). Utilizing these three NURBS (circular arc, elliptic arc and superelliptic 
arc) representations, one is able to model a deforming grid — from a circular 
arc to elliptic arc (by changing the radii) and then to a superelliptic arc (by 
changing the weights). This example is demonstrated by a 3D transition duct. 
This 3D transition duct was designed to connect a typical circular engine ex- 
haust to a high aspect ratio rectangular supersonic nozzle. The corresponding 
CFS simulation and the experimental results were first presented by the sci- 
entists in NASA Langley Research Center to investigate the effect of internal 
transition duct length on nozzle performance [Ref 41,48]. Many of the re- 
searches were concentrated on the CFS solution, however, the grids of this 
transition duct can be modelled easily by NURBS with a very concise control 
net (or control volume). Since the transitional cross sections were represented 
by the superellipse, the algorithm described in chapter two can be used to 
model each cross section with NURBS representation. Utilizing the “ruled 
volume” algorithm, one can generate the associated NURBS control volume 
for this circular-to-rectangular duct. Figure 5.15 shows the 3D grid as well as 
the associated NURBS control volume. 



Figure 5.15 Circular-to-Rectangle nozzle (NURBS control volume & grids). 

Based on the knowledge described in chapter two, one knows the fact 
that for a NURBS circular arc and elliptic arc with the same sector angle, the 
associated weights for these two NURBS representation are the same. The 
only difference is the radius of the circular arc changes to a pair of semi— major 
and semi-minor axes of the ellipse. Also, for a NURBS elliptic arc and a su- 
perelliptic arc with the same sector angle, the associated semi-major axes and 
semi-minor for these two NURBS representation are the same, however, the 
weight of the middle control point differs. Hence, one can apply these deform- 
ing relations to model this transition duct. The first procedure is to construct 
a NURBS cylinder with proper radius. A cylinder can be modelled by only 9 
by 2 control points, however, in order to model the final circular-to-rectangu- 
lar duct, more control points are needed in J direction so that it can catch the 
changing shape of the throat portion of the duct. Since only the radii change 
from the shape of circular cross sections to elliptic cross sections, the second 
step is then manipulating the associated radii of the control net at different 
cross sections to proper semi— major and semi— minor values without altering 
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the weights. During these processes, the shape of the cylinder will deform to 
the shape of circular— to-ellipse. The last step is then to increase the weights 
at the middle control point at each cross section of control net according to 
Table 2.1 for the superelliptic configurations. Figure 5.16 shows the deform- 
ing processes at different time steps for this configuration. 

This grid generating process may be tedious and time consuming, which 
makes the simulation of deforming geometry problem difficult. In previous 
cases, the NURBS representation has been utilized for modeling the ID and 
2D wave equations. One may neglect the advantage of this NURBS modeling 
approach since the resolutions used for those two cases are not extremely 
high. However, for the volume case, the CPU time required for evaluating all 
the grid points is much larger than those of the curve and surface cases. Be- 
sides, the fact that the volume grids must be re— generated at each time step 
makes it a computational intensive problem. Utilizing the NURBS for this 
deforming geometry problem will make a substantial difference. 
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After demonstrating the capability of using NURBS to model the de- 
forming geometry, an unsteady simulation using this approach is presented. 
This case simulates the flow passing through a 3D pipe which deforms at each 
time step. This 3D pipe is generated first by creating two NURBS cylinders 
with radii of r and 0 (a cylinder with radius 0 indicate a singular surface), the 
“ruled volume” algorithm is then applying to these two NURBS surfaces to 
create the final control volume. The size of control volume of this pipe is only 
2x2x9. However, in order to catch the deformation, more control points are 
needed in flow direction. The knot insertion algorithm is then used to achieve 
this goal. Figure 5.17 shows the control volume and the O type volume grid. 



Figure 5.17 NURBS control volume for a deforming pipe. 

Figures 5.18 and 5.19 show the deforming control volumes and the 
associated volume grids. The efficient NURBS evaluation algorithm (de- 
scribed in chapter 4) is utilized to generate the deforming grid at each time 
step. The unsteady solutions (courtesy Boyalakuntla [Ref 13]) obtained from 
UBI [Ref 2] flow solver are shown in Figure 5.20 ~ 5.22. 
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Figure 5.22 Grid and solution of a deforming pipe (III). 



CHAPTER VI 

CAGI: COMPUTER AIDED GRID INTERFACE 


As is discussed in previous chapters, NURBS representation has many 
important properties which are useful in many research fields such as com- 
puter graphics, CAD/CAM, and the grid generation. It has also become a stan- 
dard for geometric modelling in the industry. The transforming procedures 
and generation algorithms described in previous chapters show the versatility 
of using NURBS. The issue of geometry communication between the different 
packages is supported by the use of the standard IGES data format. In order 
to integrate all the NURBS properties and enhance the portability of geome- 
try/grid, a computer software is developed to provide the efficient tools for CFS 
analysis. The structure and the current status of this software, CAGI — Com- 
puter Aided Grid Interface, is briefly introduced in this chapter. 

CAGI Overview 

CAGI is sponsored under the grant from NASA Marshall Space Flight 
Center. The motivation for developing this software is to try to build a bridge 
between the CAD/CAM system and grid generation codes. This code is aimed 
at developing a software system which integrates CAD/CAM geometric system 
output and IGES files for widely utilized grid generation systems and expect- 
ing it to allow for fast, efficient and economical responses to the pre-proces- 
sing required in CFD applications. 
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The development of CAGI is carried out in a modular fashion by using 
the FORMS library [Ref 391 for designing the friendly Graphical User Inter- 
face (GUI). The main modules include:(i) Geometry transformer — developing 
a process of transferring widely utilized geometry entities defined in the IGES 
to NURBS format; (ii) Geometry manipulation - including NURBS geometric 
toolkits for easy manipulation of those sculptured geometries. The techniques 
of the “obstacles of using NURBS and the overcoming strategies” described in 
chapter 4 are also included to redistribute candidate geometric entities (curve, 
surface and volume) with desired distributions; (iii) Geometry generation - 
providing easy inputs for generating the NURBS sculpted geometry (curve, 
surface and volume). The overall objective of this work can be illustrated as 
Figure 6.1. 



Figure 6.1 The modules of CAGI and their links. 

The modules of “Flow Visualization” and “Flow Code” are not fully im- 
plemented in current CAGI. 
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Each module implemented in CAGI may be performed independently 
for different tasks. All of the modules are linked together with a common 
graphical user interface (GUI). This intuitive graphical user interface pro- 
vides a user friendly environment for easy and efficient manipulation and gen- 
eration of the desired sculptured geometry. The graphic visualization is devel- 
oped utilizing the GL (Graphic Library) available on the SGI machines. 
Figures 6.3 and 6.4 show the graphical user interface with the entire screen 
layout of CAGI. 

From Figures 6.3 and 6.4, one can realize that CAGI is constructed by 
creating a main panel window, an image display window, and several sub- 
panels. The main panel is the one which holds the commonly used function 
buttons or those buttons which can access the other sub-panels. The graphi- 
cal image display window located in the middle of the main panel is the win- 
dow that allows the user to interactively view the grids and the sculptured ge- 
ometry. A message panel, located below the lower-left comer of the image 
display window, is designed to report any error message or provide some use- 
ful messages for the users. Next to this message panel are those function but- 
tons which control the transformation. The user can adjust the proper sensi- 
tivity and utilize either those buttons or the mouse to control the translating, 
rotating and scaling of the displayed objects. Some other commonly used func- 
tions, such as refreshing the screen or changing the background color of the 
display window, are provided for the easy access on the global panel. 

Transformer Module 

As is discussed in the previous chapter, grid generation system and 
CAD/CAM design tools often utilized different formats which makes the trans- 
ferrial of information between systems difficult and tedious. The Geometry 
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Transformer module designed in CAGI is developed to facilitate such commu- 
nication. The theory behind this transformation is to use the NURBS as the 
interface. This is expressed in Figure 6.2. 


Discrete grid points 




IGES file 


Figure 6.2 NURBS - interface between grid system and CAD/CAM. 


Any standard IGES file (or NASA-IGES file and NINO) from a CAD/ 
CAM system can be transformed into NURBS and then output the desired 
geometric entities to discrete grid data set, such as the PLOT3D format, for 
the utilization in grid generation tools. Or from the other direction, any geom- 
etry defined by those geometry generation functions in CAGI can be output to 
a standard IGES for the use of any CAD/CAM design tools. The algorithms for 
transforming the various geometric entities into the NURBS representation 
are presented in chapters 2, 3 and 5. Table 6.1 lists those IGES geometric en- 
tities which can be transformed to NURBS definition in CAGI as well as those 
defined in NASA-IGES and NINO. 

After an IGES file is read in, the geometric entities will be shown on the 
image display window, and the associated “name list” of each entities will be 
listed on the “Entity Name List” browser located on the left side of the main 
panel as shown in Figure 6.5. If the name list is highlighted, then the 
associated entity is defined as “active” and will be displayed in graphic win- 
dow, whilp those without highlighting are defined as “inactive” and will be 
hided from the graphic window. The user can interactively turn on/off of cer- 
tain entity by selecting the corresponding name list in this browser. However, 
no matter the entities are highlighted or not, they are all in the CAGFa data- 
base unless the user “delete” them by selecting the proper function. 
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Figure 6.3 CAGI screen layout (1). 
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Figure 6.4 CAGI screen layout (2). 


Figure 6.3 CAGI screen layout (2). 
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Figure 6.5 Illustration of IGES Transformer Module in CAGI. 


Figure 6.5. Illustration of IGES Transformer Module in CAGI. 








Table 6.1 List of IGES Entities Supported by Different Packages. 
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Entity 

Number 

Entity Name 

NASA 

IGES 

NINO 

CAGI 

0 

Null Entity 

V* 

V* 

y* 

100 

Circular Arc Entity 

y 


y 0 

102 

Composite Curve Entity 

V 0 


y 

104 

Conic Arc Entity 

y 


y 0 

106 

Copious Data Entity 

y \ 


y 0 

108 

Plane Entity 




110 

Line Entity 



y 0 

112 

Parametric Spline Curve Entity 



y 

114 

Parametric Spline Surface Entity 



y 

116 

Point Entity 

y 


y 

118 

Ruled Surface Entity 



y 

120 

Surface of Revolution Entity 



y 

122 

Tabulated Cylinder Entity 


V* 

y 

124 

Transformation Matrix Entity 

V 

V* 

y 

126 

NURBS Curve Entity 

y 

V* 

y 

128 

NURBS Surface Entity 

mm 

y 

mm 

130 

Offset Curve Entity j 



y 

140 

Offset Surface Entity 



y 

141 

Boundary Entity 

y 


y 

142 

Curve on a Parametric Surface 

y 

y 

y 

143 

Bounded Surface Entity 

y 

V* 

y 

144 

Trimmed Surface Entity 



y 


The fields without u y” under NASA-IGES and NINO don’t mean these 
translators can not process those entities, it only indicates those entities must 
be transformed first to any entities allowed in NASA-IGES and NINO format. 
And the “Bounded Surface” and “Trimmed Surface” must be processed interac- 
tively in CAGI. 
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In the process of designing the sculptured geometry in the CAD/CAM 
systems, many geometric entities which are not necessary for the CFS analy- 
sis may be created. For example, a surface may be formed by revolving a 
boundary curve which was defined by set of discrete points, but the final IGES 
file created by the CAD may contain all the points and curves which are not 
needed. Reading all the entities and transforming them to NURBS definition 
requires a lot of CPU time and memory for those un-interested entities. As 
well displaying all of the entities on the graphic windows not only clogs the 
entire screen but also makes it hard for any manipulations. CAGI provides a 
solution to this problem by allowing the user to filter the entities before they 
are processed. For example, turn off entity 106 and 110 will make CAGI to 
skip processing the points and lines contained in an IGES file. And the other 
alternative to avoid the hurdle is to read all entities and store them in data- 
base, but use the “group function” to turn on/off or delete groups of entities. 
For example, the user can turn on/off all the points, curves or surfaces or cer- 
tain specific entity type by selecting the proper buttons: A status report re- 
grading the IGES entities is provided after this procedure. 

Geometry Manipulation Module 

After transforming an IGES file, one may select this module to modify 
the geometry by utilizing the NURBS properties. The NURBS information 
such as the degree, number of control points of the selected curve or surface 
will be shown in the proper position of this panel. The user can change the 
control polygon coordinates or the associated weights and visualize the change 
of the geometry interactively. Nevertheless, the most useful features for the 
NURBS representation in the applications of the numerical grid generation 
are the changing of resolutions and the distribution functions without distort- 
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ing the original geometric definition [Ref 88,90,91]. When the user select this 
module, a distribution panel also pops up at the right-lower comer of the 
main interface. This distribution panel provides eight different packing func- 
tions to define the spacing of the NURBS curves, surfaces and volumes. How- 
ever, as is aforementioned, the largest obstacle which inhibits the engineer to 
utilizing the NURBS definitions is the difficulty in controlling the spacing of 
the NURBS boundary or surface. This is because the NURBS is represented 
in the parametricform, the distribution function defined in parametric space 
may not always be reflected in the physical space. In other words, an even 
spacing in the parametric space may not result in an even distribution in 
physical space (more detail information about this problem is discussed in 
chapter 4). To solve this problem, CAGI provides options for the user to choose 
whether the distribution is defined in parametric or physical space. The sharp 
comers (discontinuous points) on the curves or surfaces can also be main- 
tained during the redistribution procedure. Moreover, the NURBS tool kits, 
such as knot insertion and degree elevation, data reduction are implemented. 
When reading the IGES file, the geometric entities may have different 
orientations which make them difficult for CFS analysis. For example, two 
adjacent surfaces may have complete different I and J directions. For this 
problem, CAGI also provides the functions which allow the user to “reverse” 
the directions, “swamp” the orientations of any selected NURBS entities. Fig- 
ure 6.6 shows the functions implemented in this module. 

Geometry Generation Module 

CAGI not only imports the geometry defined in the CAD/CAM system, 
it also has the capability to create geometry with NURBS definition for 3D 
curves, surfaces and volumes. Since the grid points / lines may not necessarily 
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pass through the NURBS control polygons (or the control nets), it has been a 
difficulty for grid generation software to model the geometry utilizing the 
NURBS information. For example, most of the engineers who don’t have a sol- 
id NURBS background don’t know how to define the NURBS control net, the 
associated weights and the knot vectors to model a sphere (or other geometry). 
However, in CAGFs generation module, this problem has been overcome. The 
user doesn’t even have to know the NURBS definition to generate the bound- 
aries and surfaces (volumes). The corresponding NURBS definitions (control 
points, weights , orders and knot vectors) will be automatically created by the 
user’s simple input. For example, to build a NURBS circular arc, the user only 
need to specify the center, the radius and the desired angles, CAGI automati- 
cally transforms these specified information to the NURBS representation. 
Same for the generation of the surface and volume. The functions implement- 
ed in this module include the generations of point, curve, surface and volume. 
For NURBS generation functions, they are listed as follows. 

NURBS Curve Generation Functions: 

The functions implemented in this module include generating of the 
NURBS “Straight line”, “Circular arc”, “Conic arc”, “Reading a NURBS curve”, 
“Extract a NURBS curve from a NURBS surface”, “Interpolate data set to a 
NURBS curve”, “Translate / Scale / Rotate of a existed NURBS curve”, “Rotate 
a point with any arbitrary axis to create a NURBS arc”, “Split a NURBS curve 
to two NURBS segments”, “Join two NURBS curves to one NURBS curve”, 
“Create a NURBS curve by trimming a existed NURBS curve”, “Creating a 
NURBS curve which lay on top of a NURBS surface”, “Duplicate a NURBS 
curve with respect to any arbitrary axis”, “Reflect a NURBS curve with re- 
spect to any arbitrary plane”. 
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NURBS Surface Generation Functions; 

The functions implemented in this module include generating of the 
“NURBS ruled surface”, “NURBS surface of revolution”, “NURBS Tabulated 
cylinder (Extruded surface)”, “Reading a NURBS surface”, “Translate / Scale / 
Rotate of a existed NURBS surface”, “Transfinite Interpolation of a NURBS 
surface”, “Duplicate a NURBS surface with respect to any arbitrary axis”, “In- 
terpolate data set toform a NURBS surface”, “Split a NURBS surface to sever- 
al NURBS sub patches”, “Join two NURBS surfaces to one single NURBS sur- 
face*, “Create a NURBS surface by trimming a existed NURBS surface”, 
“Reflect a NURBS surface with respect to any arbitrary plane” and “Offset 
NURBS surface”. 

NURBS Volume Generation Functions: 

The functions implemented in this module include generating of the 
“NURBS ruled volume”, “NURBS volume of revolution”, “NURBS Tabulated 
volume”, “Reading a NURBS volume”, “Translate / Scale / Rotate of a existed 
NURBS volume”, “Tr ans finite Interpolation of a NURBS volume”, “Duplicate 
a NURBS volume with respect to any arbitrary axis”, “Split a NURBS volume 
to several NURBS sub volumes”, “Join two NURBS volumes to one single 
NURBS volume”, “Create a NURBS volume by trimming a existed NURBS 
volume” and “Reflect a NURBS surface with respect to any arbitrary plane”. 

Figure 6.7 shows several NURBS surfaces created in this generation 


module. 
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Figure 6.6 Illustration of NURBS Manipulation Module in CAGI. 


Figure 6.6. Illustration of NURBS Manipulation Module in CAGI. 
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Figure 6.7 Illustration of NURBS Generation Module in CAGI. 


Figure 6.7 Illustration of NURBS Generation Module in CAGI. 






CHAPTER VII 

STRUCTURED GRID TOPOLOGY ON TRIMMED SURFACE 

As mentioned earlier, CAD/CAM systems generally utilize the IGES 
format for geometry Input and Output. The trimmed parametric surface (en- 
tity 144) is the most frequently used entity in the IGES format. Unfortunate- 
ly, this entity is also the most difficult entity for structured grid topology. This 
trimmed entity can not exist alone, it must be coupled with another entity — 
entity 142, which is the curve on a parametric surface entity. This entity is 
used as the trimmed curves to define the boundaries of the trimmed surface. 
Hence, to know how a trimmed surface is defined, entity 142 must be pres- 
ented first. 


Curve on a Parametric Surface Entity (Type 142) 

This entity is defined by a parametric surface, called an “untrimmed” 
surface along with the parametric curve, and identifies the curve as the entity 
lying on top of the surface [Ref 35]. 

Let the domain D of the parametric surface be a rectangle in (u,v ) space, 
then the surface S(u,v) can be defined as equation (7.1). 


S(u,v ) = (x(u,v), y(u,v), z(u,v )) 

D = { {.u,v)\u l <1 u <, u 2 i>i ^ v < v<J 


(7.1) 


The domain of the parametric curve, denoted as C(t), must lie in the range of 
the parametric domain of the surface D, and assume the parametric values 
which are used for evaluating C(t) are represented as B(t). Then B(t) must 
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satisfy the equation (7.2) shown as follows. 

B(t ) = ( u(t),v(t )) a ^ t ^ b a €E e [u^Og] (7.2) 

Then the parametric curve C(t) which lies on top of the surface can be repre- 
sented as equation (7.3) 

Of) = S(B(t)) = S(u(*),o(f)) = 

(x(u(t\ v(t)), y(u(t),v(t)\ z{u(t), v(t))) 

Hence, there are three geometric entities included in entity 142. They 
are the untrimmed surface, parametric curve B(t), and the parametric curve 
C(t). The curve B(t) is a 2D plane curve and lies on the parametric domain of 
the untrimmed surface, while the curve C(t) is a 3D curve lies on the un- 
trimmed surface in physical space. These two curves are mutually related — 
C(t) can be generated by the composite mapping of S(B(t)), while B(t) can also 
be obtained by projecting curve C(t) onto the untrimmed surface and finding 
the corresponding parametric values in the parametric domain. Because 
these two curves are both given in the IGES file, there may be conflicts if one 
like to specifies one curve differently from the other (say, obtaining the curve 
B(t) from curve C(t) vs. obtaining curve C(t) from B(t)). This means projecting 
the curve C(t) onto the untrimmed surface may not be able to obtain a curve 
exactly the same as B(t), or on the other hand, evaluating the B(t) on surface 
S(u,v) may not be able to create the same curve as C(t). Obviously, this will 
lead to the accuracy problem. Hence, in an IGES file, another flag is used to 
control which curve should be used to avoid this problem. This flag is denoted 
as “PREF” (preferred flag). If PREF is 0, IGES does not specify which curve to 
use, if PREF is 1, the B(t) should be used, if PREF is 2, then curve C(t) is pre- 
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ferred, if PREF is 3, then both C(t) and B(t) are equally preferred. Figure 7.1 


illustrates the entity 142. 



Figure 7.1 Illustration of IGES entity 142. 


After knowing the format of entity 142, the trimmed surface can be de- 
fined. In the IGES format, the trimmed surface is defined by a surface, called 
the untrimmed surface, along with two types of boundaries — the outer 
boundary and inner boundary which are lying on the untrimmed surface. 
Both types of boundaries are defined with the format of entity 142. Hence, the 
unt rimm ed surface which the boundaries (entity 142) lie in must be identical 
to the one defined in the trimmed surface entity. If the number of inner 
boundaries is greater than or equal to zero, however, there is only one outer 
boundary defined. If there is no outer boundary curve specified, the bound- 
aries of the untrimmed surface will be used as the default outer boundary (one 
closed curve), and if there are no inner boundaries and no outer boundary spe- 
cified, then the trimmed surface will be equal to the untrimmed surface. Fur- 
thermore, all inner / outer boundaries must be simple closed curves and mutu- 
ally disjoint. Each of the inner boundaries must lie in the interior of the outer 
boundary. Hence, the domain of the trimmed surface is defined as the com- 
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mon region of the interior of the outer boundary and the exterior of each of the 
inner boundaries and includes the boundary curves. Figure 7.2 illustrates the 
definition of this entity. 
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Figure 7.2 Illustration of IGES entity 144. 


The four-sided structured patch requires the number of grid points on 
the opposite boundary be the same; it also requires that the total number of 
grid points be equal to the multiplication of the grid points of two adjacent 
boundaries. Apply these two criteria to the trimmed surface shown in Figure 
7.2, one can immediately understand that the trimmed surface is not a real 
structured patch. Indeed, it is shown only for display. As matter of fact, it is 
plotted by the “scan lines” and the “filling polygons” algorithms [Ref 27]. Be- 
cause all the inner and outer boundaries are simple closed curves, any “scan” 
line which parallels to I (or J) direction of the parametric domain will intersect 
all the boundaries with even number of points. Out of all these intersection 
points, those line segments starting with odd number will be plotted. For ex- 
ample, plotting the line segments connecting the points of (1, 2), (3, 4), ... (N— 1, 
N). Same procedure should be applied to the J (or 7) direction for another set 
of line segments. After this algorithm, the trimmed region (bounded by the 
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interior of outer boundary and the exterior of inner boundaries) can then be 
displayed, however, it can not be utilized for structured grid generation. 

To grid the trimmed surface for a “structured” topology is difficult. The 
reasons for this are summarized as follows: 

Complexity 

As is discussed in previous section, the trimmed surface (entity 144) de- 
fined in the IGES file may contain many inner boundaries, also the bound- 
aries could be any shape (convex or concave curve). For the structured grid, it 
requires 4-sided boundaries. However, in many of the trimmed surfaces, 
there are 3-sided, 5-sided, or even more than 6-sided boundaries. This phe- 
nomenon makes it very complicated for structured grid generation. It is al- 
most impossible for the “automatic” structured grid generation for a compli- 
cated trimmed surface. Figure 7.1 demonstrates the complexity of the 
structure grid generation, in which the trimmed surface contains four inner 
boundaries and the outer boundary is not convex. 

Inconsistency 

The trimmed surface shown in Figure 7.2 is a simple geometry. It con- 
tains only one inner boundary and one outer boundary. However, one can not 
create the “ruled” surface for the trimmed region by just connecting the two 
boundaries. The reason of that is the two boundaries have different directions 
and the starting points of the two curves are far away from each other. The 
directions and the location of the starting points of the boundaries create the 
inconsistency problem which increses the difficulty for griding the trimmed 
surface. 
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Accuracy 

Another difficulty for griding the trimmed surface is related to the accu- 
racy problem. Taking Figure 7.2 as an example, the inner boundary should be 
a closed curve, however, the starting point and the ending point of this curve 
are not exactly identical. The gap between these two points may come from 
the tolerance defined in different design systems. This accuracy problem also 
happens to the gap between the common edges of two adjacent trimmed sur- 
faces. Also, all the trimming boundaries are defined with the format of entity 
142 (curve on parametric surface), if improper parametric values are used, the 
sharp comers (discontinuous points) on the boundaries may be lost, and this 
will lead to an inaccurate geometric definition. Because the trimmed surface 
was designed in CAD/CAM systems, the accuracy problem already existed 
when the geometry was imported to grid generation system. This increases 
the difficulty for generating grids for the trimmed surfaces. 

Grid Topology Requirement 

As it has been mentioned it is very difficult to design an automatic grid- 
ing package for a structured trimmed surface. Even though the simple 
trimmed surface may be handled automatically, the returned structured grid 
may not be useful for all the solution algorithms. This is because different 
solution algorithms may require different grid topologies - such as H type to 
avoid the singularity point. Due to the preferred grid topologies required, the 
difficulty for griding the trimmed surface increased. 

All these difficulties indicate the design of an automatic structured grid 
generator for the trimmed surface (entity 144) would be impracticable. In or- 
der to overcome these obstacles and also let the engineer determine the de- 
sired grid topology, GAG/ currently designs several panels for the user to 
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create structured grids interactively. The functions included in these panels 
are described as follows. 


Resolution Control 

The first panel is developed for the purpose of visualizing the trimmed 
entity. As is discussed in previous section, the untrimmed surface and the 
trimming boundaries are presented as NURBS format. Hence, the resolutions 
will control how the geometry be displayed on computer screen. If the resolu- 
tions of the trimming entities are too low, the shape of the geometry may not 
be truly displayed. On the other hand, if the default resolutions are set too 
large, the program will take much CPU time and more memory to evaluate all 
entities. Since the engineer is only interested in the “trimmed” region, setting 
large resolutions for untrimmed surfaces and the trimming boundaries will 
simply waste the computer resources. Therefore, the motivation for designing 
this panel is to provide easy tools to increase the resolutions for a better un- 
derstanding of the geometry, or decrease the resolutions to reduce the memory 
load after the trimmed region is completed. Figure 7.3 shows the layout of this 
panel. 
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Figure 7.3 Resolution control panel. 
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Breaking Point Define 

In order to generate the structured grid on the trimmed surfaces, it is 
necessary to do the domain decomposition, and the first step of the decomposi- 
tion is defining the breaking points on the trimming boundaries. The func- 
tions implemented in this panel allow the user move the slider bar and inter- 
actively see a breaking point moving on the trimming boundary both in the 
parametric and physical space. The parametric value and the coordinate of 
the breaking point will be updated whenever the user moves the slider bar. 
One can press the “Add” (or “Remove”) button to insert (or delete) the breaking 
points from the database. Moving the slider bar may fail to obtain the proper 
parametric values to catch the sharp comers on the trimming boundary. 
Hence, a function is implemented to allow the user obtain all possible sharp 
comers by finding the proper parametric values. These functions are illus- 
trated in Figure 7.4. 



Figure 7.4 Break point define panel. 

Edge Define 

After the breaking points are defined, the user is allowed to use any two 
breaking points to define an edge. The edge could be a straight line in the 
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parametric domain or a curve along the trimming boundary. While defining 
the edges, the user not only can decide the resolutions but also can chose the 
desired distribution by selecting various packing functions. The distribution 
is referred to the physical space which means the user’s desired distribution is 
reflected on the trimmed surface region. The commonly used functions, such 
as “Reverse the edge direction”, “Combine two edges” are also implemented for 
a easy manipulation. Figure 7.5 shows these functions and the corresponding 
edges in the parametric space as well as the physical space. 
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Figure 7.5 Edge define panel. 


Structured Patch Define 

After the edges are proper defined, the user can select four edges which 
form a closed region and then apply the TFI algorithm to these four edges to 
create a structured patch in parametric space. One can obtain one piece of 
trimmed surface in the physical space by utilizing the parametric patch as the 
parametric values and evaluating them with the untrimmed surface. Because 
the TFI is applied to the edges on the parametric space, the resulting patch in 
parametric space may not always result in a smooth surface in physical space. 
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Hence, the elliptic function is implemented so that the user can improve the 
quality of the final trimmed surface. These functions and example are all il- 
lustrated in Figure 7.6. 
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Figure 7.6 TFI patch define panel. 


After the user generates the pieces of trimmed surfaces, the output 
functions designed in CAGI can be selected for outputing the surfaces to dif- 
ferent data formats (for example, the PLOT3D format), the grid generation 
package, such as the GENIE++ [Kef 61], can use there surfaces for the volume 
grid generation. 

These interactive procedures allow the user to generate the structured 
grids with desired topology on the trimmed surfaces. The algorithm which can 
catch the discontinuous points (described in chapter four) has been installed 
in “Break point define” panel. It allows the user obtain the proper parametric 
values for the sharp comer in physical surface, hence, the proper break points 
can be defined. The re-parameterization algorithm (described in chapter 
four) which can generate the distribution accurately on physical space has 
been installed in the “Edge define panel”. It allows the user to specify the de- 
sired point distribution on the edges in trimmed surfaces. The smooth algo- 
rithm based on the elliptic function has been implemented in the “ TFI patch 
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panel” for a satisfactory structured grid. Another structured trimmed surface 
example obtained from these procedures are shown as follows. 



Figure 7.7 Trimmed surface example. 



CHAPTER VIII 

SUMMARY AND CONCLUSIONS 


The geometry modeling techniques used in Computer Aided Geometry 
Design have been extended and applied to numerical grid generation for com- 
plex real world applications. The general transformation algorithms which 
precisely convert the Non-NURBS entities to NURBS representations have 
been developed. These algorithms can be utilized to bridge the gap between 
the grid generation and the CAD/CAM systems. Also, the two-way linkage 
between the CAD/CAM and grid systems has been achieved through the IGES 
format utilizing entities 126 (NURBS curve) and 128 (surface). The formula- 
tion of NURBS has been extended from curves, surfaces to full 3D NURBS 
control volumes to model the CFS configurations. The development of the re- 
parameterization schemes and their influence on the curve / surface / volume 
grid distributions is demonstrated by computational examples. The applica- 
tions of these re-parameterization techniques to precise grid distribution con- 
trol with accurate geometry fidelity have been demonstrated. An efficient 
NURBS evaluation routine has been developed to facilitate the entire grid 
generation process. In addition, the applications of NURBS to dynamic grid 
generation presented in this study have proven the versatility of NURBS in 
the CFS simulation processes. 

The algorithms developed in this study have been utilized for various 
grid generation related research areas. For example, the interpolation rou- 
tines developed under this research have been applied to Yang’s adaptation 
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[Ref 83,84,85] as well as Shih’s Turbomachinery grid generator “TIGER” [Ref 
52,56]. The NURBS transformation and evaluation routines also have been 
applied to Boyalakuntla’s unsteady simulation for temporally deforming ge- 
ometries [Ref 13]. These algorithms and the associated CFS applications have 
been published in the Journal of “ Computer-Aided Design ” and “ Applied 
Mathematics and Computations * [Ref 51,94]. 

The NURBS generation and manipulation algorithms have been incor- 
porated as modules in the CAGI system. The motivation for developing CAGI 
has been to follow the IGES standard, and incorporate the NURBS represen- 
tations for the two way communications between the CAD/CAM and grid sys- 
tems. CAGI provides the capability of processing the IGES files and converts 
various IGES entities to NURBS representations. These NURBS representa- 
tions are then evaluated, thus preparing the boundaries and surfaces for grid 
generation. The generation module allows the generation of surface and vol- 
ume grids based on NURBS generation algorithms. The NURBS database 
adopted in CAGI allows the geometries / grids to be output to the IGES format 
with entities 126 (NURBS curve) and 128 (NURBS surface). Thus, the grid 
output is compatible with the needs of the CAD/CAM system. The intuitive 
graphical user interface provides a user friendly environment for easy and ef- 
ficient manipulation and generation of the desired sculptured geometry. In 
summary, the algorithms developed and incorporated using CAGD techniques 
in CAGI have demonstrated a significant time saving involved in surface prep- 
arations associated with grid generation. Also, CAGI now provides a two-way 
linkage between the CAD and grid systems. 

The CAGD techniques have been utilized in CAD/CAM systems for de- 
cades, many of them are well developed and documented in the literature. 
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Some of these algorithms are adopted in this study. These algorithms include 
the (NURBS) knot insertion, degree elevation, splitting / joining algorithms 
[Ref 42-47,76,77] and the evaluation of the normalized BSpline basis func- 
tions [Ref 18]. These algorithms have been realized and implemented in this 
research. Also, the data reduction programs [Ref 38] and the FORMS Library 
[Ref 39], available in the public domain are implemented. Applications of 
these existing algorithms make the grid generation process easier. For exam- 
ple, applying the perturbations to the NURBS control net mdses the simula- 
tion of temporally deforming problem efficient. The modeling of volume grid 
using NURBS control polygons reduces the memory storage requirement. 

The grid generation process generally has a more restrictive require- 
ment for modeling the geometry than that of CAD/CAM system. For example, 
the surface of revolution option for grid generation is not always for a full rev- 
olution (360® revolution). Also the generation of the conic arc usually requires 
different rotation angle and the knowledge of the semi— major and semi— minor 
axis information. The generation of a circular arc may require the different 
staring and ending angles. Hence, the transformation algorithms for circular 
arc, conic arc, surface of revolution to NURBS (described in Chapter Two) are 
all enhanced and generalized for the grid generation procedures. Also, the 
grid adaptation algorithm [Ref 83-85] is enhanced by replacing the interpola- 
tion algorithm with the new reparameterization algorithm (described in 
Chapter Four) for the precise grid distribution control. When the NASA IGES 
committee defined the NASA IGES specification, the committee accepted vari- 
ous suggestions from this study (such as the exclusion of the IGES plane enti- 
ty, entity number 108, from the specification). 
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The development of many new algorithms related to NURBS is the key 
contribution of this study. These new algorithms include the transformation 
of superelliptic arc to a NURBS curve, the transformation of cubic curve (sur- 
face) to a BSpline curve (surface) and the transformation of a composite curve 
to one NURBS curve (described in Chapter Two). Another contribution is re- 
lated to the generation algorithms. For example, the generation of a cascad- 
ing surface by NURBS and various NURBS volume generating functions (de- 
scribed in Chapter Three) have not been found in related literature. The new 
NURBS evaluation algorithm presented in Chapter Four provides a competi- 
tive alternative to the de Boor algorithm (in terms of memory and computa- 
tional time). The new reparameterization algorithm described in Chapter 
Four already received positive feedback after it was published in the Journal 
of “ Computer Aided Design” [Ref 94] . 

Considerable advantages have been realized in the development of 
these algorithms as demonstrated in this study. However, various research 
issues remain to be addressed. They are described as follows. 

Even though the re-parameterization allows the user to control the grid 
points (lines) distributed along the NURBS entity in the physical space, in 
many CFS simulations, especially for the viscous and turbulent problems, this 
distribution function may not be adequate. The orthogonal grid near the 
boundaries needs to be considered. Since the NURBS is an algebraic method, 
the orthogonality is an important issue for further study. 

The unstructured grid and the hybrid grid techniques have become pop- 
ular and important in today’s CFS simulation. These two techniques provide 
the flexibility to construct a complex configuration easily. As the examples 
shown in chapter three, one can generate the unstructured / hybrid grids in 
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the 2D parametric domain and then utilize NURBS parametric property to 
evaluate the final 3D unstructured / hybrid grid in the physical space. Howev- 
er, the parametrization problem occurs. Even though the behaviors (aspect 
ratio, skewness ...) of the triangles or hybrid grids are satisfactory in paramet- 
ric space, this does not necessarily result in the same grid quality in the physi- 
cal space. The re-parameterization algorithm for structured grids described 
in chapter four will not be applicable to the unstructured (triangulated) grid. 
Hence, sthe NURBS re-parameterization algorithm for unstructured and hy- 
brid grid generation is an important issue for further research. 

Another open question is related to the NURBS deforming geometry. 
As is described in chapters two and five, NURBS is composed of the control 
points, orders, knot vectors and weights. Changing any of this information 
will lead to a new NURBS entity. However, the difficulty lies in the fact that 
for certain geometries, there is not sufficient information to model them by 
NURBS exactly. For example, exact modeling of the sine wave or a bell- 
shaped curve with concise NURBS control polygons still has not been formu- 
lated. For the deforming geometry problems, if the deforming shapes are spe- 
cified within different time steps, for example, a semi-circular change to a 
sine wave within time interval 0 ~0.5 sec, and from the sine wave to a bell- 
shaped curve within the time interval 0.5 ~ 1.0, the NURBS modeling tech- 
nique may fail and a new combination change of the NURBS information 
must be found. 

A final important issue relates to the automatic structured grid gen- 
eration for the trimmed surfaces defined in the IGES file. The strategies for 
handling the structured trimmed surface presented in this study allows the 
engineers to interactively generate the desired structured grid around the 
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trimmed region. However, in some of the industrial applications, an IGES file 
may contain thousands of trimmed surfaces. It is necessary to enhance the 
procedures so that the overall structured grid generation on the trimmed sur- 
face can be reduced. 

Although there is still room for further work, as presented, CAGI has 
already successfully bridged the gap between grid systems and CAD/CAM sys- 
tems, and it has been effectively used in industrial applications. The NURBS 
volume generation and manipulation algorithms will be extended to enhance 
the CAGI to be a solid modeling package in the future work. 
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This file was produced by CAGI (Computer Aided Grid Interface) betal .0 S 1 

„19HERC/NSF Miss. State, 10Hsphere.igs, G 1 

17HMSU/ERC CAGI V1.0.19HCAGIV1.0, 10/10/94,32,38,6, G 2 

308,15,llHUnspecified,l,l,4HINCH,32,32,13H951109.173219,0.000001, G 3 

2.00000000000000,8HRobertYu,llHNSF/ERC/MSU,9,0,13H951109. 173219; G 4 

128 1 0 0 0 0 0 0 00000000D 1 

128 0 0 52 0 0 0 B_SPL_SF ID 2 

128, 4, 8, 2, 2, 0, 0, 0, 0, 0, IP 1 

O.OOOOOOE+OO, O.OOOOOOE+OO, O.OOOOOOE+OO, 5.000000E-01, IP 2 

5.000000E-01, l.OOOOOOE+OO, l.OOOOOOE+OO, l.OOOOOOE+OO, IP 3 

O.OOOOOOE+OO, O.OOOOOOE+OO, O.OOOOOOE+OO, 2.500000E-01, IP 4 

2.500000E-01, 5.000000E-01, 5.000000E-01, 7.500000E-01, IP 5 

7.500000E-01, l.OOOOOOE+OO, l.OOOOOOE+OO, l.OOOOOOE+OO, IP 6 

l.OOOOOOE+OO, 7.071068E-01, l.OOOOOOE+OO, 7.071068E-01, IP 7 

1 .OOOOOOE+OO, 7.071068E-01, 5.000000E-01, 7.071068E-01, IP 8 

5.000000E-01, 7.071068E-01, l.OOOOOOE+OO, 7.071068E-01, IP 9 

l.OOOOOOE+OO, 7.071068E-01, l.OOOOOOE+OO, 7.071068E-01, IP 10 

5.000000E-01, 7.071068E-01, 5.000000E-01, 7.071068E-01, IP 11 

l.OOOOOOE+OO, 7.071068E-01, l.OOOOOOE+OO, 7.071068E-01, IP 12 

l.OOOOOOE+OO, 7.071068E-01,5.000000E-01,7.071068E-01, IP 13 

5.000000E-01, 7.071068E-01, l.OOOOOOE+OO, 7.071068E-01, IP 14 

l.OOOOOOE+OO, 7.071068E-01, l.OOOOOOE+OO, 7.071068E-01, IP 15 

5.000000E-01, 7.071068E-01, 5.000000E-01, 7.071068E-01, IP 16 

l.OOOOOOE+OO, 7.071068E-01, l.OOOOOOE+OO, 7.071068E-01, IP 17 

1 .OOOOOOE+OO, 1 .OOOOOOE+OO, O.OOOOOOE+OO, O.OOOOOOE+OO, IP 18 

1 .OOOOOOE+OO, 1 .OOOOOOE+OO, O.OOOOOOE+OO, -4.37 1 1 39E-08, IP 19 

1.000000E+00, 0.OOOOOOE+OO-l.OOOOOOE+00, l.OOOOOOE+OO, IP 20 

0.000000E+00,-l .000000E+00,-8.742278B-08, O.OOOOOOE+OO, IP 21 

l.OOOOOOE+OO, O.OOOOOOE+OO, O.OOOOOOE+OO, l.OOOOOOE+OO, IP 22 

l.OOOOOOE+OO, 1 .OOOOOOE+O0,-4.371 139E-08, l.OOOOOOE+OO, IP 23 

l.OOOOOOE+OO, -l.OOOOOOE+OO, l.OOOOOOE+OO, l.OOOOOOE+OO, IP 24 

-1.000000E+00.-8.742278E-08, O.OOOOOOE+OO, l.OOOOOOE+OO, IP 25 

O.OOOOOOE+OO, O.OOOOOOE+OO, l.OOOOOOE+OO, -4.37 1139E-08, IP 26 

1 .OOOOOOE+00,-4.37 1 1 39E-08,-4.37 1 1 39E-08, l.OOOOOOE+OO, IP 27 

-1.000000E+00,-4.371139E-08, l.OOOOOOE+OO-l. OOOOOOE+OO, IP 28 

-8.742278E-08, O.OOOOOOE+OO, l.OOOOOOE+OO, O.OOOOOOE+OO, IP 29 

O.OOOOOOE+OO, 1 .OOOOOOE+OO-1 .OOOOOOE+OO, l.OOOOOOE+OO, IP 30 

-4.37 1 1 39E-08.-1 .OOOOOOE+OO, l.OOOOOOE+OO, -l.OOOOOOE+OO, IP 31 

-l.OOOOOOE+OO, 9.999999E-01,-1.000000E+00,-8.742278E-08, IP 32 

O.OOOOOOE+OO, l.OOOOOOE+OO, O.OOOOOOE+OO, O.OOOOOOE+OO, IP 33 

1 .OOOOOOE+O0,-1 .0O000OE+O0,-8.742278E-O8,-4.37 1 1 39E-08, IP 34 

-l.OOOOOOE+OO,- -8.742278E-08.-1.000000E+00-1 .OOOOOOE+OO, IP 35 

-8.742278E-08,-l .00000OE+0O.-8.742278E-O8, O.OOOOOOE+OO, IP 36 

l.OOOOOOE+OO, O.OOOOOOE+OO, O.OOOOOOE+OO, l.OOOOOOE+OO, IP 37 

-1 .000000E+00,-l .000000E+00,-4.371 139E-08-9.999999E-01 , IP 38 

-1 .OOOOOOE+OO ,-l .OOOOOOE+O0,-9.999999E-Ol ,-i .OOOOOOE+OO, IP 39 

-1 . OOOOOOE+OO, -8.742278E-08, O.OOOOOOE+OO, l.OOOOOOE+OO, IP 40 

O.OOOOOOE+OO, O.OOOOOOE+OO, 1 .OOOOOOE+OO, 1 .3 1 1342E-07, IP 41 
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-1 .000000E+O0,-4.37 1 1 39E-08, 1.311342E-07 -l.OOOOOOE+OO, 
-l.OOOOOOE+OO, 1 .3 1 1342E-07.-1 .OOOOOOE+OO.-l .OOOOOOE+OO, 
-8.742278E-08, 0.OOOOOOE+OO, l.OOOOOOE+OO, O.OOOOOOE+OO, 
O.OOOOOOE+OO, l.OOOOOOE+OO, 1.000000E+00-9.999999E-01, 
-4.37 1 139E-08, 1 .000000E+00,-9.999999E-Cl ,-l .OOOOOOE+OO, 
1 .OOOOOOE+O0,-9.999998E-Ol,-1 .OOOOOOE+OO, -8.742278E-08, 
O.OOOOOOE+OO, l.OOOOOOE+OO, O.OOOOOOE+OO, O.OOOOOOE+OO, 
l.OOOOOOE+OO, l.OOOOOOE+OO, 1. 748456E-07, -4.371 139E-08, 
l.OOOOOOE+OO, 1.748456E-07, -l.OOOOOOE+OO, l.OOOOOOE+OO, 

1 .748456E-07-1 ,OOOOOOE+O0,-8.742278E-08, O.OOOOOOE+OO, 
O.OOOOOOE+OO, l.OOOOOOE+OO, O.OOOOOOE+OO, l.OOOOOOE+OO; 

S 1G 4D 2P 52 


IP 

42 

IP 

43 

IP 

44 

IP 

45 

IP 

46 

IP 

47 

IP 

48 

IP 

49 

IP 

50 

IP 

51 

IP 

52 

T 

1 
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#The number of control points of the surface 
52 

#The two degrees of the surface 
21 

#The two knot vectors 
0.000000 
0.000000 
0.000000 
0.500000 
0.500000 

1.000000 

1.000000 

1.000000 
0.000000 
0.000000 

1.000000 

1.000000 

#The coordinates of the control points and weight 
- 1.950000 0.000000 0.000000 1.000000 
- 1.950000 1.200000 0.000000 0.707107 
0.000000 1.200000 0.000000 1.000000 
1.950000 1.200000 0.000000 0.707107 
1.950000 0.000000 0.000000 1.000000 
- 8.000000 - 0.000001 0.000000 1.000000 
- 8.000001 8.000000 0.000000 0.707107 
0.000000 8.000000 0.000000 1.000000 

8.000000 8.000000 0.000000 0.707107 

8.000000 0.000000 0.000000 1.000000 
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Introduction 


Why is CAGI needed? 

The computational Field Simulation (CFS) for physical problems has become more practical due to the 
progress made in computer memory availability and the computer speed. This CFS process generally 
involves the steps: 

1. Pre-processing: includes numerical grid generation, boundary condition set up and definition 
involving sculptured geometry description of all solid components and field boundaries / surface. 

2. Processing: Numerical solving a set of non-linear partial differential equations at the grid point 
developed by a pre-processor. 

3. Post-processing— Numerical visualization of flow field properties simulated by the solution 
processor. 

In this process, the sculptured modeling and the grid generation are the most time consuming and 
laborious work. It is necessary to develop a system which can reduce the geometry contraction and grid 
generation time for the CFS process. 

For a routine application of CFS, in an industrial environment, the overall response time for CFS must 
be reduce considerably. As noted by several scientists: 

" ...The industrial requirement is for reliable one hour grid generation turnaround for 
one-time geometries when ran by designers. The system must include CAD-to-Grid links 
which resolve tolerance issues and produce grids with a quality good enough for the flow 
solver. The designer has to feel that the grid generation processes is under control and is 
predictable." 

— David Ives, Robert Miller, William Siddons and Kevin VanDyke. May 1995 

In many of today’s industrial applications, most of the geometrical configurations of interest to practical 
problems are designed using a CAD/CAM system. Setting up the communication between CAD/CAM 
design tools or other grid generation system would be the ideal situation. And this idea was pointed out 
by the NASA Steering Committee on Surface Modeling and Grid Generation on 1992. Unfortunately, 
the CAD/CAM systems have many different geometry output formats which force the designer to spend 
a gread deal of time manipulating geometrical entitites in order to achieve a useful sculptured 
geometrical description for grid generation. In addition, there is a danger of losing the fidelity of the 
geometry in this process of data transfer between different Input/Output (I/O) formats. The other issure 
related field simulation is the grid quality. The quality of the grid affects the accuracy of the solution and 
the computation time. It may be necessary to reconstruct the grids for a more satisfactory result after 
obtaining the first solution. This reconstruction procedure involves a change of either resolution (the size 
of the grids) or the spacing (the distance between grid points) functins. However, this proess is tedious 
and very time consuming; this is especially true for a complex geometry. 



To bridge the gap between the CAD/CAM systems and the grid generation systems, it is necessary to 
establish communicatin paths so that the geometries and grids defined within these two systems can be 
linked with each other. For most of the CAD/CAM systems, the Initial Graphics Exchange Specificatin 
(IGES) is a widely accepted standard for geometry exchange. Most CAD/CAM systems support the 
IGES format as an Input/Output of resulting geometries. And for the grid generation part, there are many 
approaches for representing sculptured geometry/grid, such as the Bezier curve/surface, parametric 
cubic curve/surface. Hermit representation, Transfinite Interpolation (also referred as TFI), ... and so on. 
Among these representations, the Non Uniform Rational BSpline (NURBS) has been widely utilized in 
many applications. NURBS is getting popular because it has a powerful features, such as the local 
control property, variation diminishing and convex hull, .. etc. Also, the geometry tool kits, like the 
curve/surface interpolation, data reduction, degree elevation, knot insertion and splitting are all 
well-developed. These features make the NURBS very useful not only in the CAD/CAM packages but 
also in Grid Generation systems. And what important is - the IGES file already included the NURBS 
curve (entity 126) and surface (entity 128). Therefore, if a software package can read in an IGES file, 
and convert any Non NURBS entities to NURBS curve/surface, then the geometry defined in IGES can 
be utilized by grid system. Or in the other way, if the geometry/grid defined in grid system is 
represented as NURBS, then it is possible to output the grids to an IGES format with entity 126 and 
128. If this can be done, then the geometry and grid are communitable, and the construction time for the 
pre-posser can be reduced. CAGI is aimed at this goal. It is a package with NURBS database. It can read 
in the IGES file and transforms the geometrical definitions to NURBS, or define the grids with NURBS 
representation and output the grid with IGES data. 


What is CAGI? 

CAGI stand for Computer-Aided Grid Interface. It is a grid genration package with NURBS database. 
The Graphic User Interface (GUI) is made by utilizing the FORM Library, and the SGI Graphics 
Library is utilized for the graphic display. CAGI can either read the standard IGES format or generate 
grids form NURBS definition. The representation of NURBS has been extended from curve (ID), 
surface (2D) to volume (3D) definition. This project is sponsored by the NASA Marshall space flight 
center. Any suggestion or request can be forwarded to the technical monitors, Mr. Ted Benjamin and 
Robert Williams. 

Currently, CAGI contains several modules: 

O Geometry Generation Module 

This module allows the user to generate the grids from NURBS definition. The user can 
create points, curves, surfaces or volume by the different NURBS options. 

O Geometry Manipulation Module 

This Module allows the user to manipulate the selected geometric entities (either curves, 
surfaces or volume). The user is allowed to change the orentation, the resolutions, or even 
the distribution of the selected entity easily by clicking the proper button or slider bars. Since 
CAGI transforms the geometry to the NURBS definition, this module takes advantage of 
NURBS properties/tools such as knot insertion, degree elevation, data reduction and even the 
alternation of the location of control polygon and weight to change the shape of the 
geometry. 



O Volume Grid Module & Geometry Viewing 

These two modules are temperady not fully function. They are under construction. 


Contact Information 


General Contact: 


Dr. Bharat K. Soni 
P. O. Box9627 

Mississippi State, MS 39762-9627 
(601)325-2647 or (601)325-8278 


bsoni @ erc.msstate.edu 


Technical Problem or CAGI Question: 


Dr. Tzu-Yi (Robert) Yu 
P. O. Box 9627 

Mississippi State, MS 39762-9627 
Tel:(60 1)325-2467, Fax:(601)325-7692 



yu @ erc.msstate.edu 


CAGI request: 


Mr. Ted Benjamin 

ED32, NASA Marshall Space Flight Center 

Huntsville, AL 35812 

Tel: (205)544-9402, Fax:(205)544-1215 


tedb @ tyrell.msfc.nasa.gov 


Questions or comments about this document: 



crystal @ erc.msstate.edu 


Organization of the CAGI Manual 

In this manual, you will find descriptions of the commands in CAGI and instructions on how to 
use them. You will also find tutorials that give you step-by-step instructions on some examples. 



The manual is arranged in the following manner: 

O CAGI Basics introduces the CAGI Interface. 

O CAGI Men Buttons describes the CAGI options and commands. 

O Tutorials/Demos help you get started learning the functions of CAGI. 
O Reference allows you to find the information you need more easily. 


Navigating the CAGI Manual 


In addition to hyperlinked texts, there are navigation buttons at the end of each document. The 
navigation buttons are described as follows: 


takes you to the Table of Contents 


takes you to Index 


takes you to the previous file in linear sections of the document; takes you to the 

parent file in non-linear sections 




Hsi takes you to the beginning of the Manual 


takes you to the next file in in linear sections of the document; takes you to the 
first child in non-linear sections 






CAGI Basics 


m Starting CAGI 

O Limition: CAGI is implemented by C language, the memory allocation function is used. Unless 
the user runs a very big IGES file or create several huge grids to consume the memory, otherwise, 
otherwise, it can be executed in a small box like Personal IRIS. However, since the Graphics 
functions and the GUI are implemented by utilizing the SGI Graphics library, this program can 
only be executed on SGI machine which supports the GL library. 

O Execution: To run CAGI, simply type "cagi 


CAGI Main Window 


This image shows an orientation of CAGI Screen Layout. It contains the geometry 
displayed in the graphics window and the entire GUI layout with several main modules. 
Click here for a large image. For each general section of the CAGI is discussed in detail 
as follows: 



• Entity Name List 

The Entity Name List panel is placed at the left side of CAGI. The panle lists the names of all the 
entities in the database (either created from an IGES file or generated from the build function). If a 
geometry (any grid) is created, then it will be displayed in the graphics window with the proper 
named (with a highlight bar with light blue color) showed in this Entity Name List panel. If an 
entity name is highlighted with a light blue bar, then that entity is referred as active and will be 
displayed on the graphichs window. If an entity name is not highlighted, the entity is referred as 
inactive and will not be shown in the graphic window. The user can turn on or off entities by using 
the mouse to highlight or unhighlight an entity’s name. For example, the use can put the mouse to 
any of the entity name listed in this panel, and use the right mouse to click the name to toggle 
on/off of the entity. 

• Message Window 

The Message Window is placed at the bottom left comer of CAGI. It is designed to show any error, 
warning or status messages for the user. 

• Transformation Panel 

The last panel in the main CAGI window is related to the tranformation and graphics window. The 
Transformation Panel provides another options for rotating objects. 





O The 111 buttons indicate rotation with respect to the positive axes. 
O The ill buttons indicate rotation with respect to the negative axes. 


There are three pairs of rotation buttons availabl e, each pair is coupled with one principal 

direction (i.e, the X, Y or Z axis). For example, issllia&l press this rotation button 
indicate the user like to rotate the active objects with respect to positive or negative X axis. 
The rotation buttons for the same direction are mutual excluded (either positive or negative 
direction). The user can select the rotation option for different principal axis at the same 
time. 


O The MiWWM&l button allows the user to Rotate the objects continually once the rotation 
direction is decided. In other words, after the rotation buttons are selected, press this button 
will make the geometry objects continually rotated. 

The sensitivity of the Rotation, Translation and Scale may be adjusted by using the buttons 
found under the Sensitivity portion of the panel: 

O The il buttons decrease the sensitivity of Translation, Rotation, and Scaling. 

O The HH buttons increase the sensitivity of Translation, Rotation, and Scaling. 


The single arrow allows a minor decrease (or increase), while the double arrow button 
allows a major decrease (or increase). 


The Transformation Panel also has functions which allow the user to adjust the center of rotation: 

O The button causes CAGI to calculate the center of all the active entites and bring 

the active entities to the center of the graphic window. This button is used when the entities 
are outside the graphics window after rotation, translation or scaling. CAGI (actually, the 
GL library) uses the origin by default for transformation. This can be very inconvenientif the 
entities are far away from the origin. The rotation would be very hard to control if the 
rotation radius is very huge. The Reset_V button counteracts this problem. 


O The button tells CAGI to use the center of the entites as the rotating center. This 

provides an each way to view the geometry even when it’s far away from the origin. To 
calculate the current geometric center, press the reset_v button. 


O The laMBll button displaces the xyz-coordinate frame. 


This Transformation Panel also allows the user to change the background color of the graphic 
window. 


BackG 


O The fewM BBBB i M l button allows the background color of the graphic window to be changed. 
Once the button is selected, a color panel will appear. The user may change the color by 
clicking on the desired color. The left and right arrow keys found on either side of the 





pop-up screen give more color palettes to from which to choose. 

Other functions available in the Transformation Panel include: 

O The button allows the user to take a "snapshot" or picture of the entities displayed 

in the graphics window. Once the button is pressed, the snapshot icon will appear. The user 
then must: 

1. Place the cursor over the snapshot icon until the cursor turns into a camera. 

Note: If the snapshot icon appears in an area where it is undesirable, the user can 
move the icon by placing the cursor over it and selecting the Alt and the F7 keys on 
the keyboard. The cursor will become crosshairs and the user can move the icon to a 
spot which is more convenient 

2. Hold down the middle mouse button. This keeps the camera cursor even when the 
cursor is moved from the snapshot icon. 

3. Select the right mouse button. 

4. Select Redraw Rubberband by moving through the menu until it is highlighted and 
releasing the right mouse button. 

5. Select and hold the left mouse button. Move the mouse to draw a rubberband around 
the entity you wish to take a snapshot of. Once the rubberband is drawn, it may be 
resized using the horizontal, vertical and comer cursors or moved using the crosshair 
cursor. 

6. Select New file name by using the same procedure outlined in step 4. Once the File: 
menu appears, enter the name and save it with either an .rgb or .gif file extension. 

7. Select Save and Exit using the same procedure outlined in step 4. 

Note: If the user wishes to take more than one picture, the user can choose Save as 
<file> and then go to steps 4 through 7 for each picture. 

Note: For more information on using SnapShot, please refer to the Man Pages 


O The ImBMiI button makes the graphics window resize to fill the entire CAGI screen. The 
user may exit this view by pressing the Esc key on the keyboard. 

O The button allows the user to delete unwanted entities. Generally, the user can 

delete the geometric entity by using the button designed in the Geometry Manipulation 
Panel. Also, the user can turn on/off of an entity by placing the cursor to the Entity Named 
List Panel and clicking the entity. However, if the user like to delete an object without going 
to manipulation panel, one can use this button to reach the goal. To utilize this function, 
press the DelObj icon (and the button will keep on pressed). Next pick the entity name from 
the Entity Name List. To stop deleting objects, click the DelObj button again (and the button 
will be released). After the button is released, the user can use the mouse to click the entities 







The functions of this pop-up menu allow the user to input or output files in various formats. These 
functions are discussed below: 


■ Input /Output 


The Input /Output field 
allows the user to either 
read in or write out a file. 
The user can click the 
button to switch for either 
Input or Output function. 


■ File Type 

The File Type field allows 
the user to select the file 
format. The user can select 
the IGES file format as 
input or create an output 
file for all the active 
entities in either IGES, 
FAST(Plot3D) or NURBS 
format. 


If a Plot3D file is to be read in, 
the user should go to the 
Geometry Generation Module and select either curve or surface for input. For reading the Plot3D 
file, please refer to Interpolate and Plot3d functions in Geometry Generation Panel. 

A Plot3D volume file can only be read after the Volume Grid Module has been utilized. 

■ Data Form 

The Data Form field allows the user to select either Formatted or Unformatted data for 
output. 

■ File Name 

The File Name input field allows the user to type in the file name. Or if the user clicks the 
file name listed in the browser, the selected file name will be displayed in this input field. 

After selecting the proper data formats and keying the file name, the user should press the 
Accept button to process the case, or press the Close button to close this form. 






Examples of using these menu buttons are shown as follows: 

1. Input Example: Select INPUT, IGES data file, type in the name of the IGES file (for 
the IGES data, there is no "Unformatted" file), and then press the Accept button. After 
the CAGI process the IGES file completely, press Close to close this form. 

2. Output Example: Suppose CAGI already creates 2 curves and 3 surfaces (and these 
entities are also actived), if the user would like to output these entities to Plot3d files. 
Then select the OUTPUT, PLOT3D, UNFORMAT and type in the desired file name 
for output, press the Accept for output. If the file name is already existed, CAGI will 
pop up a warming message and ask the user confirm for overwriting the file or quitting 
the output process. 


O The 



menu button cause a pop-up menu to appear: 


All the functions provided in this pop-up menu will change the properties of all active entities. 
These functions are discussed below: 



Shade 

The Shade button will set all 
the active surfaces, i.e. the ones 
which are visisble in the 
graphics window with the entity 
name hilighted in the Entity 
Name List browser, to be 
shaded. 

Wireframe 

The Wireframe button turns all 
active shaded surfaces to 
wireframe (only the grid lines 
are plotted) ones. 

Sur clr 

The Sur clr button causes a 
color panel to pop-up menu to 
appear, from which the user can 
choose the desired color for the 
active surfaces. 

Cur clr 

Similar to Sur clr, this Cur clr 
button causes a color panel to 
pop-up from which the user can 
choose the desired color for the 
active curves. 


Backface 














The Backface button is used to 
turn backfacing polygon 
removal on and off. Please refer 
to the GL man pages for more 
detail information. This 
function is only effected when 
the Transparency is applied to a 
shaded object. 

■ Reshape 

The Reshape button in 
conjunction with the NI, NJ and 
NK input fields is used to 
changed the resolutions of all 
active entities. To utilize this 
function, keyin the desired 
values into die NI, NJ and NK 
input fields. Next, select this 
Reshape button to change the 
resolutions of the active 
entities. 

■ Single/Double 

The Single/Double button 
controls the precision of the 
output file format. If Single is 
active, then the output files will 
be outputted as single precision. 
If Double is active (after the 
user press the Single, the button 
will become Double) then the 
output files will be output as 
double precision. This is 
important if the user like output 
the geometry / grids and like to 
visualize them by FAST, since 
it can not read in Unformatted 
Double precision file, the user 
is advised to use Single 
precision (the default) for 
output. 

■ SurNet/ 

The SurNet/ button controls the 
display of the NURBS control 
nets associated with the active 
surfaces. If this button is 
pressed, all the active surfaces 




will display the associated 
NURBS control nets. The user 
can click the button again to 
release to turn off the function. 

■ CurNet/ 

The CurNet/ button controls the 
display of the NURBS control 
nets associated with the active 
NURBS curves. Press the 
CurNet button to turn on this 
function and release the button 
to turn the function off. 

■ CurNetcolor 

The CurNetcolor button causes 
a color panel to pop-up menu to 
appear, from which the user can 
choose the desired color of the 
NURBS curves control nets. 

■ SurNetcoIor 

The SurNetcoIor button causes 
a color panel to pop-up menu to 
appear, from which the user can 
choose the desired color of the 
NURBS control nets. 

■ adjust knot 

The adjust knot button causes a 
pop-up menu to appear, from 
which the user can adjust the 
knot vector (s) for the NURBS 
entity (curve/surface or 
volume) based on the 
arc-length. This button is used 
when a NURBS entity has a 
"bad" distribution. 

■ Reserved 

The Reserved buttons are for 
future functions. If users have 
any suggestions for adding 
extra functionalities, please 
send e-mail to Dr. Tzu-Yi 
"Robert" Yu at 
yu @ erc.msstate.edu 




O The 


menu button causes a new pop-up menu to appear. 


All the functions provided in this pop-up menu will modify the shading parameters for all active 
and shaded entities. Using the right mouse button to press the Option button on the pop-up menu, 
it will display four choices -Material Editor, Light Editor, Model Editor, and Close. Each choice 
is discussed below: 

1. Material Editor 

■ Ambient 

The Ambient field specifies the ambient reflectance 
of the material. It is followed by three color slider 
bars with floating point values, typically in the range 
0.0 through 1.0. These values specify red, green, and 
blue reflectances. If the values are (1.0, 1.0, 1.0), it 
indicates the Ambient light is pure while. The 
default values are as: 0.40, 0.20, and 0.00. The user 
can modify these values by moving the the desired 
slider bars. 

■ Emission 

The Emission field specifies the color of light 
emitted by the material. It is followed by three 
floating point values, typically in the range 0.0 
through 1.0, which specify red, green, and blue 
emitted light levels. The default values for this field 
are: 0.50, 0.00, and 0.00. The user can modify these 
values by moving the color slider bars. 


■ Diffuse 

The Diffuse field specifies the diffuse relfectance of 
the material. It is followed by three floating point 
values, typically in the range 0.0 through 1.0, which 
specify red, green, and blue diffuse reflectances. The 
default values for this field are set as: 0.90, 0.50, and 
0.00. The user can modify these values by moving 
the slider bars. 


■ Specular 

The Specular field specifies the specular reflectance 
of the material. It is followed by three floating point 
values, typically in the range 0.0 through 1.0, which 
specify red, green, and blue specular reflectances. The default values for this field are: 
0.90 , 0.90 , and 0.00 . The user can modify these values by using the slider bars. 

■ Shininess 

The Shininess field specifies the specular scattering exponent, or the shininess, of the 







material. It is followed by a single floating point value, typically in the range 0.0 
through 128.0, which specifies the shininess. Higher values result in smaller, hence 
more shiny, specular highlights. If this value is set as 0.0, it will effectively disables 
specular reflection. The default value for this field is: 50.00. The user can modify this 
value by using the slider bar. 

■ Transparency 

The Transparency field specifies the transparency of the material. It is followed by a 
single floating point value, typically in the range 0.0 through 1.0, which specifies the 
transparency. The default value for this field is: 1 . 00 . The user can modify this value 
by using the slider bar. 

■ Reset 

The Reset button located in the right upper comer of this form enable the user to 
specify the desired material properties (like how to set the Ambient, Emission, Diffuse 
and Specular values). Using the right button of the mouse to select the desired 
material. Currently, the material of Brass, ShinyBrass, Pewter, Silver, Gold, Shinygold, 
Plaster, Cyanplastic, Greyplastic, Yellowplastic, Redplastic, Greenplastic, Blueplastic, 
Pinkplastic, Lavpolstone, Brownplostone are available. 


2. Light Editor 




■ Ambient 

The Ambient field specifies the ambient light 
associated with the light source. It is followed by 
three floating point values, typically in the range 0.0 
through 1.0, which specify the red, green, and blue 
ambient light levels. The default values for this field 
are: 0.20 , 0.20 , and 0.60 . The user can modify 
these values by using the slider bars. 


■ Lcolor 

The Lcolor field specifies the color and intensity of 
the light that is emitted from the light source. It is 
followed by three floating point values, typically in 
the range 0.0 through 1.0, which specify the levels of 
red green and blue light emitted from the light 
source. The default values for this field are: 1.00, 
1.00, and 1.00. The user can modify these values by 
using the slider bars. 


■ X, Y, Z 

The X, Y, and Z fields specify the position of the 
light source in the graphics window. Each is 
followed by a floating point value which specifies 
the position in terms of object-coordinates. 

■ Local/Infinity View 

The Local/Infinity View button allows the user to 
choose where the light source be treated as locally 
distant or infinitely distant. The default position for 
the button is Local View. 

■ Active Light 

The Active Light field allows the user to choose the 
number of active light sources. Values of 1 through 7 are available. The use can add or 
delete light sources using either the Add or Del buttons. The arrows key allows the user 
the select the light source for manipulating its properties. 



3. Model Editor 






■ Ambient 

The Ambient field specifies an additional ambient 
light lievel that is associated with the entire graphic 
window, rather than with a light source. This light is 
added to the ambient light in the scene. It is followed 
by three floating point values, typically in the range 
0.0 through 1.0, which specify the red, green, and 
blue ambient light levels. The default values for this 
field are: 0.40, 0.20 , and 0.00. The user can modify 
these values by using the slider bars. 

■ Attenuation 

The Attenuation field specifies the constant and 
linear attenuation factors associated with all 
non-infinite tight sources. It is followed by two 
floating point values in the range 0.0 through 1.0. 
The first attenuation factor is used to directly reduce 
the effect of a tight source on objects in the graphics 
window. The second tight source factor specifies 
attenuation that is proportional to the distance of the 
tight source from die object(s) being lighted. The 
default values for this field are 1.00 for constant 
attenuation and 0.00 for linear attenuation. 

■ Twoside 

The Twoside button specifies whether lighting 
calculations are done assuming that only frontfacing 
polygons are visible, or are corrected for each 
polygon based on whether it is frontfacing or 
backfacing. If the button is pushed, a lighting model 
that is correct for both frontfacing and backfacing 
polygons is utilized. If the button is not pressed (the 
default state), a lighting model that is correct only 
for polygons whose visible face is the facet for 
which normals have been provided is utilized. 



Note:: The user can use the man page for "man lmdef ' for more detail information. 


O The 



menu button causes a new pop-up menu (see below) to appear: 


The functions designed in this form, mainly dealing with the geometric database and the IGES 
entities. They allow the user to tun on/off (or delete) groups of geometric entities contained in the 
database, and also allow to " filter " an IGES file before the process. The functions are discussed as 
follows: 




■ Entities Statistic Table 

The Entities Statistic Table displays 
the statistics of the entities contained 
in an IGES file which has been 
process by CAGI. 


■ Inactive all 

The Inactive all field contains three 
buttons. These buttons allows the 
user to inactivate all points, curves 
or surfaces by selecting the Point, 
Curve or Surface buttons, 
respectively. 

■ Delete all 

Similar to the Inactive all, the 
functions allow the user to delete 
(erase from database) all points, 
curves or surfaces by selecting the 
Point, Curve or Surface buttons, 
respectively. 


■ Group active 

The Group active field allows the 
user to activate or inactivate all 
entity groups displayed in the entity 
choice — the violet box between the 
Inactive and Active buttons. The 
user can select the entity groups by 
clicking on the entity choice box 
with the right mouse button. It will 
display a list of the entity groups. To 
select a group, go through the list until the desired group is highlighted and release the right 
mouse button. Until the desired group is selected, the user can press Active button or 
Inactive button to make the select geometric entities to be active or inactive. 

■ Group Del 

The Group Del field allows the user to delete all entity groups selected from the entity 
choice (described in the previous function). The user can select the entity groups by clicking 
on the entity choice(the violet box) with the right mouse button to display a list of the entity 
groups. To select a group, go through the list until the desired group is highlighted and 
release the right mouse button. To delete the selected entity group, push the Delete button. 

■ Reset 

The Reset button allows the user to reset the Entities Statistic Table by calculating the 
current status of all entities in the database. 



Quit 















The Quit button allows the user to exit the pop-up menu. 

■ IGES Entity Filter 

An IGES file generally contains a lot of entities — they could be points, curves or surfaces. 
However, for a CFD application, the points (or some boundary curves) may not be that 
important. Before reading an IGES file, the user has options to " filter " out the undesired 
entities by releasing the button with the entity name on it. After releasing the button, the 
name of the button will be "Closed" to indicate that even an IGES file contains this entity, it 
will not be process by CAGI. The user can select different entities from curve or surface 
defined in the IGES file. 


O The 


button allows the user to access the help menu. Since the current 


documentation is written in the web site, currently, this help button is not function. 


O The! 


button allows the user to terminate the CAGI program. A warming message 


will appear for the user to confirm or cancel the exit. 


•CAGI Modules 

Cagi has been designed with four modules. Please click on the hypertext to read in-depth descriptions of 
the functions available in each module. 

O Geometry Generation Module 

Selecting this module allows the user to generate the NURBS entities, i.e. a point, NURBS 
curves, surfaces or even NURBS volume. 

O Geometry Manipulation Module 

Selecting this module allows the user to manipulate the selected geometric entities (curve, 
surface or volume). Since CAGI transforms the geometry to a NURBS format, this module 
takes advantage of NURBS properties such as knot insertion, degree elevation and the 
alternation of the control polygon and weight function to change the shape of the geometry. 

O Volume Grid Module 

Selecting this modules allosws the user to generate volume grids. 

O Geometry Viewing 

Selecting this module allows the user to visualize the solution or complex geometry. 








Introduction 


Why is CAGI needed? 

The computational Field Simulation (CFS) for physical problems has become more practical due to the 
progress made in computer memory availability and the computer speed. This CFS process generally 
involves the steps: 

1. Pre-processing: includes numerical grid generation, boundary condition set up and definition 
involving sculptured geometry description of all solid components and field boundaries / surface. 

2. Processing: Numerical solving a set of non-linear partial differential equations at the grid point 
developed by a pre-processor. 

3. Post-processing— Numerical visualization of flow field properties simulated by the solution 
processor. 

In this process, the sculptured modeling and the grid generation are the most time consuming and 
laborious work. It is necessary to develop a system which can reduce the geometry contraction and grid 
generation time for the CFS process. 

For a routine application of CFS, in an industrial environment, the overall response time for CFS must 
be reduce considerably. As noted by several scientists: 

" ...The industrial requirement is for reliable one hour grid generation turnaround for 
one-time geometries when ran by designers. The system must include CAD-to-Grid links 
which resolve tolerance issues and produce grids with a quality good enough for the flow 
solver. The designer has to feel that the grid generation processes is under control and is 
predictable." 

— David Ives, Robert Miller, William Siddons and Kevin VanDyke. May 1995 

In many of today’s industrial applications, most of the geometrical configurations of interest to practical 
problems are designed using a CAD/CAM system. Setting up the communication between CAD/CAM 
design tools or other grid generation system would be the ideal situation. And this idea was pointed out 
by the NASA Steering Committee on Surface Modeling and Grid Generation on 1992. Unfortunately, 
the CAD/CAM systems have many different geometry output formats which force the designer to spend 
a gread deal of time manipulating geometrical entitites in order to achieve a useful sculptured 
geometrical description for grid generation. In addition, there is a danger of losing the fidelity of the 
geometry in this process of data transfer between different Input/Output (I/O) formats. The other issure 
related field simulation is the grid quality. The quality of the grid affects the accuracy of the solution and 
the computation time. It may be necessary to reconstruct the grids for a more satisfactory result after 
obtaining the first solution. This reconstruction procedure involves a change of either resolution (the size 
of the grids) or the spacing (the distance between grid points) functins. However, this proess is tedious 
and very time consuming; this is especially true for a complex geometry. 



To bridge the gap between the CAD/CAM systems and the grid generation systems, it is necessary to 
establish communicatin paths so that the geometries and grids defined within these two systems can be 
linked with each other. For most of the CAD/CAM systems, the Initial Graphics Exchange Specificatin 
(IGES) is a widely accepted standard for geometry exchange. Most CAD/CAM systems support the 
IGES format as an Input/Output of resulting geometries. And for the grid generation part, there are many 
approaches for representing sculptured geometry/grid, such as the Bezier curve/surface, parametric 
cubic curve/surface, Hermit representation, Transfinite Interpolation (also referred as TFI), ... and so on. 
Among these representations, the Non Uniform Rational BSpline (NURBS) has been widely utilized in 
many applications. NURBS is getting popular because it has a powerful features, such as the local 
control property, variation diminishing and convex hull, .. etc. Also, the geometry tool kits, like the 
curve/surface interpolation, data reduction, degree elevation, knot insertion and splitting are all 
well-developed. These features make the NURBS very useful not only in the CAD/CAM packages but 
also in Grid Generation systems. And what important is — the IGES file already included the NURBS 
curve (entity 126) and surface (entity 128). Therefore, if a software package can read in an IGES file, 
and convert any Non NURBS entities to NURBS curve/surface, then the geometry defined in IGES can 
be utilized by grid system. Or in the other way, if the geometry/grid defined in grid system is 
represented as NURBS, then it is possible to output the grids to an IGES format with entity 126 and 
128. If this can be done, then the geometry and grid are communitable, and the construction time for the 
pre-posser can be reduced. CAGI is aimed at this goal. It is a package with NURBS database. It can read 
in the IGES file and transforms the geometrical definitions to NURBS, or define the grids with NURBS 
representation and output the grid with IGES data. 


What is CAGI? 

CAGI stand for Computer-Aided Grid Interface. It is a grid genration package with NURBS database. 
The Graphic User Interface (GUI) is made by utilizing the FORM Library, and the SGI Graphics 
Library is utilized for the graphic display. CAGI can either read the standard IGES format or generate 
grids form NURBS definition. The representation of NURBS has been extended from curve (ID), 
surface (2D) to volume (3D) definition. This project is sponsored by the NASA Marshall space flight 
center. Any suggestion or request can be forwarded to the technical monitors, Mr. Ted Benjamin and 
Robert Williams. 

Currently, CAGI contains several modules: 

O Geometry Generation Module 

This module allows the user to generate the grids from NURBS definition. The user can 
create points, curves, surfaces or volume by the different NURBS options. 

O Geometry Manipulation Module 

This Module allows the user to manipulate the selected geometric entities (either curves, 
surfaces or volume). The user is allowed to change the orentation, the resolutions, or even 
the distribution of the selected entity easily by clicking the proper button or slider bars. Since 
CAGI transforms the geometry to the NURBS definition, this module takes advantage of 
NURBS properties/tools such as knot insertion, degree elevation, data reduction and even the 
alternation of the location of control polygon and weight to change the shape of the 
geometry. 







The manual is arranged in the following manner: 

O CAGI Basics introduces the CAGI Interface. 

O CAGI Men Buttons describes the CAGI options and commands. 

O Tutorials/Demos help you get started learning the functions of CAGI. 
O Reference allows you to find the information you need more easily. 





Appendix A. 1 


Genie++ : Brief Description 




GENIE ++ - General Grid Generation System 

INTRODUCTION 

NASA maintains an applications-oriented computational fluid dynamics (CFD) effort complemen- 
tary to and in support of aerodynamic-propulsion design and test activities. This is especially true 
at NASA/MSFC where the goal is to advance and optimize present and future liquid-fueled rocket 
engines. Numerical grid generation plays a significant role in fluid flow simulation utilizing CFD. 
The first step in a CFD simulation is the generation of an appropriate grid.. The geometry of interest 
must be accurately modeled and the points distributed in an efficient and smooth manner. These 
constraints often cause this step to be often one of the most time-consuming. Several grid generation 
codes of increasing capability have developed. During 1983-85 SVTGD2D-3D was developed and 
operated in a batch mode. As graphics workstations progressed in power and capability, IN- 
GRID2D-3D was developed during 1985-88 to use available interactive techniques to speed the 
process. As an outgrowth of these programs and other research activities Genie was developed dur- 
ing 1988-91 as a semi interactive grid generation package. From 1992 to the present Genie ++ has 
been under development as a completely interactive grid generation system. Genie ++ has demon- 
strated the capability to generate grids about very complex configurations of interest to MSFC with 
complete geometric fidelity. Thus, current development efforts concerning Genie ++ deal with tech- 
niques to decrease the labor time required and to enhance the fidelity of the geometry representation. 
The transfer of information directly from a CAD system to the grid generation system has the poten- 
tial to facilitate the fulfillment of these two goals. Therefore the computer aided grid interface 
(CAGI) has been developed and the Initial Graphics Exchange Specification (IGES) translator im- 
plemented in Genie ++ . The IGES translator is fully compatable with the NAS A-IGES standard and 
allows direct transfer of information from a CAD system to Genie ++ . The development of CAGI 
has been supported by NASA MSFC and complex configurations of current NASA interest have 
been used as test cases to validate the grid generation system. 
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Genie ++ PROGRAM DESCRIPTION 


The computer code GENIE ++ is a continuously evolving grid system containing a multitude of prov- 
en geometry/grid techniques. The generation process in GEN1E ++ follows earlier versions. The pro- 
cess uses several techniques either separately or in combination to quickly and economically gener- 
ate sculptured/analytical geometry descriptions and grids for arbitrary geometries. The 
computational mesh is formed by using an appropriate algebraic method. Grid clustering is accom- 
plished with either exponential or hyperbolic tangent routines which allow the user to specify a de- 
sired point distribution. Grid smoothing can be accomplished by using an elliptic solver with proper 
forcing functions. B-spline and Non-Uniform Rational B-splines (NURBS) algorithms are used 
for surface definition and redistribution. The built-in sculptured/analytical geometry definition 
with desired distribution of points, automatic Bezier curve/surface generation for interior bound- 
aries/surfaces, surface redistribution or remapping based on NURBS weighted Lagrange/Hermite 
transfinite interpolation methods, interactive geometry/grid manipulation modules, and on-line 
graphical visualization of the generation process are salient features of this system, which result in 
a significant time savings for a given geometry/grid application. 

One recently developed capability has proven to be very useful in several applications. Even when 
geometry data is obtained directly from a CAD system, it is generally not in a form suitable for gen- 
eration. Often patches defining portions of a surface do not match with their background surface. 
Also surface/surface intersections are not defined, and sometimes do not fully intersect. In this case 
an extrapolation must be performed. This is performed automatically when a surface/surface inter- 
section is performed in Genie ++ . Once this intersection curve is determined the job is still not com- 
plete. The intersection curve must conform to a series of I, J or K constant lines in the computational 
space. Experience has indicated that the typical practice of breaking a surface into a larger number 
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Original Wing Fairing 


} 



Wing Fairing After Intersection 


Figure lb. Fuselage Patch Before & After Blending 



Figure 2a. Original IGES Data for Veins 
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Figure 2b. Selected Surfaces of Veined Elbow Volume Grid 
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